I! "RD-A133  554 ADVANCED  TELEPROCESSING  SYSTEMS(U)  CALIFORNIA  UNIV  LOS  i/2 
1  ANGELES  DEPT  OF  COMPUTER  SCIENCE  L  KLEINROCK  20  SEP  82 

I  UCLA-ENG-83-24  MDA903-82-C-0064 


UNCLASSIFIED 


F/G  i 7/2.  1  NL 


UNIVERSITY  OF  CALIFORNIA,  LOS  ANGELES 


SEMI-ANNUAL  TECHNICAL  REPORTS 
_  Sponsored  by  the 

Defense  Advanced  Research  Projects  Agency 
Contract  Nos.:  DAHC-I5-C-0368,  MDA  903-77-C-0272:  DARPA  Order  No.  2496 


COMPUTER  NETWORK  RESEARCH 


DATES 

DDC  ACCESSION  NUMBER 

August  1969  to 

February  1970: 

AD  705  149 

August  1970: 

AD  711  342 

June  1971: 

AD  727  989 

December  1971: 

AD  739  705 

June  1972: 

AD  746  509 

December  1972: 

AD  756  708 

June  1973: 

AD  769  706 

December  1973: 

AD  A 0041 67 

June  1974: 

AD  A008422 

December  1974: 

AD  A0 16823 

June  1975: 

AD  A02067I 

December  1975: 

AD  A02C9I4 

June  1976: 

AD  A034I7I 
(Final) 

ADV  ANCED  TELEPROCESSING  SYSTEMS 


June  1976  to 
December  1976: 
June  1977: 

June  1978: 
September  1979: 
March  1980: 
September  1981: 


AD  A039018 
AD  A 047 496 
AD  A077404 
AD  A08I938 
AD  A 088839 
to  be  assigned 


ADVANCED  TELEPROCESSING  SYSTEMS 
Annual  Technical  Report 
September  30,  1982 


Contract  Number:  MDA  903-82-C-0064 
DARPA  Order  Number:  2496 
Contract  Period:  October  1, 1981  to  September  30,  1983 
Report  Period:  October  1,  1981  to  September  30,  1982 

Principal  Investigator:  Leonard  Kleinrock 

Co-Principal  Investigator:  Mario  Gerla 

(213)  825-2543 


Computer  Science  Department 
School  of  Engineering  and  Applied  Science 
University  of  California,  Los  Angeles 


Sponsored  by 

DEFENSE  ADVANCED  RESEARCH  PROJECTS  AGENCY 


The  views  and  conclusions  contained  in  this  document  are  those  of  the  authors  and  should  not  be 
interpreted  as  necessarily  representing  the  official  policies,  either  express  or  implied,  of  the  Defense 
Advanced  Research  Projects  Agency  or  the  United  States  Government. 


Slew  MTV  CLASSIFICATION  OF  THIS  SSCI  fSSwi  0«*  towN) 


REPORT  DOCUMENTATION  PAGE 

READ  INSTRUCTIONS 

BEFORE  COMPLETING  FORM 

1.  RKCIFIKN  T*S  Catalog  NUMBS  A 

4.  TITLE  fans  SaMlllal 

ADVANCED  TELEPROCESSING  SYSTEMS: 

ANNUAL  TECHNICAL  REPORT 

*•  TVFK  OF  RERONT  A  RKRIOO  COVKRKO 

Annual  Technical 

1  OCT  81  through  30  SEPT  82 

A.  RKRFORMINO  ORO.  RKROAT  NUMBER 

?.  AUTNORfaJ 

Leonard  Kleinrock 

4.  CONTRACT  OR  GRANT  NUMBSRfaJ 

MDA  903-82-C-0064 

».  RIRFOAMING  OROANIZATION  namk  ano  ADDRESS 

School  of  Engineering  and  Appoied  Science 
University  of  California 

Los  Anueles,  CA  90024 

16.  R NOONAN  ELEMENT,  RROJKSt^  TASK  . 

AREA  A  WORK  UNIT  NUMBSIIS 

DARPA  Order  No.  2496 

II.  CONTNOLLINO  OFFICS  NANS  ANO  AOONSSS 

Defense  Advanced  Research  Projects  Agency 

1400  Wilson  Boulevard 

Arlington,  VA  22209 

IZ.  RERONT  OATS 

September  30,  1982 

151 

14.  NONITONINO  AOCNCV  NANI  A  AODRKSSfff  Bfffaranf  fraai  CmntrmJIlHI  out  mm) 

It.  SECURITY  CLASS,  (ml  (Ufa  NRII) 

ISa.  OEC^ASSIFICATION/OOWNaRAOINO 

1C  DISTRIBUTION  STATENENT  (ml  tflla  Umpmrtj 


Approved  for  Public  Release;  Distribution  unlimited 


IT.  DISTRIBUTION  STATKMKNT  (ml  Ma  NllrH  MmS  Im  IlMk  SB.  If  BlffaraM  tnm  Kmpmrt) 


I*,  surrlkmbntary  notks 


It.  KIY  BOROS  (CmIImn  an  rmm  aMa  It  nmmmmimry  aN  IfHitliy  fer  MaaA  m—  >ar> 

Computer  communication  networks.  Packet  switching.  Multi-hop  networks.  Packet 
radio  networks,  Channel  access  Protocols,  ALOHA  network,  Spatial-TDMA, 
rude-CSMA,  Resource  sharing,  Multi-access  communications.  Computer 
communications,  Broadcast  communications ,  Channel  sharing. 

tO.  ABSTRACT  fCanlftwa  an  ravaraa  «Ma  If  aaaaaaary  aM  IBanilfr  Sr  MaaS  wtar>  " 

;)This  Annual  Technical  Report  covers  research)  carried  out  by  the 
Advanced  Teleprocessing  Systems  Group  at  UCLA  under  DARFA  Contract 
No.  MDA  903-82-C-0064 1 covering  the  period  from  October  1,  1981  to 
September  30,  1982.  _ 

'  y>  r 


DO  I  j2J*T*  1473  COITION  OF  I  NOV  00  IS  OBSOLETE 

%m  oioa-LF  Oisasoi  HMrT.aggTpTvn  r 

IKCURITY  CLASSIFICATION  OF  THIS  RAOE  flNlan  Data  tmimrmt) . \ 


KCUfllTV  CL  AMI  fie  AVION  OF  THIS  PAOCPHim  Of 


xThia  contract  has  three  primary  designated  research  areas:  packet  radio 
systems,  resource  sharing  and  allocation,  and  distributed  processing  and 
control.  _ 

This  report  contains  the  abstracts  of  the  publications  jnhich  suonarize 
our-  research  results  in  those  areasduring  this  annual  period,  followed 
by  the  main  body  of  the  report  which  is  devoted  to  a  study  of  channel 
access  protocols  that  are  executed  by  the  nodes  of  a  network  to  schedule 
their  transmissions  on  multi-access  broadcast  channel.  In  particular, 
the  main  body  consists  of  «e-Eh.D.  dissertation. iay  Randolph  Nelson, 
/"^Channel  Access  Protocols  for  Hilti-Scp  Broadcast  Packet  Radio  Net¬ 
works,  vl^orfwctedunder  the  supervision  of  Professor  Leonard  KLeinrock 
(Principal  Imstl^tor^for  this  research)  .V-This  work  discusses  same 
new  channel  access  protocols  useful  for  mobile  radio  networks.  Included 
is  an  analysis  of  slotted  ALOHA  and  same  tight  bounds  on  the  performance 
of  all  possible  protocols  in  a  mobile  environment.  A  new  protocol  known 
as  Mft  is  also  introduced  and  analyzed;  It  is  collision-free 

protocol  which  allows  one  to  find  the  minimus  possible  delay  foe  mes¬ 
sages  in  a  network  which  is  perfectly  scheduled.  Yet  another  protocol 
known  as  jfcife  is  introduced  and  studied.  Here  a  terminal  is  per¬ 
mitted  to  transmit  even  when  he  detects  another  terminal's  transmission, 
under  certain  circumstances  which  are  discussed  in  detail;  the  cases 
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INTRODUCTION 

This  Annual  Technical  Report  covers  research  carried  out  by  the  Advanced  Teleprocessing  Sys¬ 
tems  Group  at  UCLA  under  DARPA  Contract  No.  MDA  903-82-C-0064  covering  the  period  from 
October  1,  1981  through  September  30,  1982.  Under  this  contract  we  have  three  designated  tasks  as 
follows: 


TASK  I.  PACKET  RADIO  SYSTEMS 

The  extension  of  our  analytic  and  design  techniques  to  modern  multi-hop  packet  radio 
networks  will  be  studied.  The  applications  and  extensions  include  access  methods, 
large  network  control  and  management,  queueing  network  models,  approximation 
methods,  capture  phenomena,  conflict-free  algorithms,  reliability,  routing  procedures, 
topological  studies,  TDMA  in  a  multi-hop  environment,  and  multiplexing  methods. 

TASK  II.  RESOURCE  SHARING  AND  ALLOCATION 

Extended  concepts  of  “power”  in  networks  will  be  studied.  The  extensions  include 
more  complex  topologies  and  configurations,  extended  queueing  disciplines,  general 
distributions,  other  definitions  of  power,  effects  of  varying  the  traffic  matrix,  fairness. 

The  problems  of  large  scale  internetting  with  respect  to  resource  allocation  and  sharing 
will  also  be  studied  further. 

TASK  III.  DISTRIBUTED  PROCESSING  AND  CONTROL 

Overall  principles  of  distributed  processing  and  distributed  control  will  be  studied.  The 
issues  of  sequencing  in  data  base  updates,  distributed  control,  and  distributed  process¬ 
ing  (involving  the  calculation  of  concurrency  of  processing)  are  the  subjects  of  concern 
here. 

A  major  contribution  of  our  research  during  this  period  is  contained  in  Reference  5  listed 
below,  namely  “Channel  Access  Protocols  for  Multi-Hop  Broadcast  Packet  Radio  Networks”,  by  Ran¬ 
dolph  Nelson.  This  Ph.D.  dissertation  (supervisored  by  Professor  Leonard  Kleinrock)  makes  a  number 
of  contributions  to  the  field  of  access  protocols  for  multi-hop  systems.  In  this  report,  slotted  ALOHA 
is  analyzed  in  a  mobile  packet  radio  network.  Furthermore,  a  new  protocol  known  as  Spatial-TDMA  is 
described  and  analyzed;  this  is  a  collision-free  system  for  which  we  find  the  mean  response  time  and 
allows  us  to  determine  the  absolute  minimum  delay  averaged  over  all  messages  in  the  network.  Yet 


another  protocol  referred  to  as  Rude-CSMA  is  discussed  and  analyzed  and  is  an  access  protocol  which 
permits  terminals  to  transmit  even  after  they  sense  a  busy  channel.  The  performance  of  the  system  and 
the  conditions  under  which  it  is  useful  are  discussed.  The  entire  dissertation  is  reproduced  as  the  main 
body  of  this  report  immediately  following  the  list  of  research  publications  below.  This  list  of  publica¬ 
tions  summarizes  the  results  of  this  annual  period  and  the  abstract  of  each  paper  is  given  along  with  the 
reference  itself. 


RESEARCH  PUBLICATIONS 


1.  Nelson,  R.  and  L.  Kleinrock,  “Spatial  TDMA:  A  Collision  Free  Multi-Hop  Channel  Access 
Protocol,”  ICC  ‘82  Conference  Proceedings,  Philadelphia,  Pennsylvania,  June  13-17,  1982,  pp. 
1C. 4.1  -  1C.4.4,  submitted  to  IEEE  Transactions  on  Communications. 

In  this  paper  we  define  a  broadcast  channel  access  protocol  called  Spatial-TDMA 
which  is  designed  specifically  to  operate  in  a  multi-hop  packet  radio  environ¬ 
ment  where  nodes  are  assumed  to  be  stationary.  The  protocol  defined  assigns 
transmission  rights  to  nodes  in  the  network  in  a  local  TDMA  fashion  and  is 
collision-free.  Methods  for  determining  slot  allocations  and  slot  sizes  are 
developed  and  a  solution  is  found  for  determining  the  assignment  of  capacities 
for  the  links  of  the  network  that  minimizes  the  average  delay  of  messages  in 
the  system. 

2.  Levy,  H.,  “Record  Recognition  and  Duplicate  Elimination  of  Personnel  Lists,”  Computer 
Science  Department,  University  of  California,  Los  Angeles,  Master’s  Thesis,  June  1982. 

4 

Personnel  lists  consist  of  records  which  contain  descriptions  of  individuals. 

Each  such  description  can  be  the  person’s  name,  address,  birth  date,  profes¬ 
sion,  etc.  Due  to  misspelling  or  other  reasons  these  descriptions  are  not 
mathematically  exact  in  their  computer  description.  For  this  reason,  it  is  very 
likely  that  one  person  can  be  represented  in  different  lists  by  records  which  are 
not  exactly  identical  to  each  other. 

In  the  presence  of  this  property  it  is  desired  to  recognize  personnel  records. 

More  specifically,  it  is  required  to  identify  if  two  records  ri,r2  represent  the 
same  person. 

This  problem  is  confronted  in  this  work.  While  the  solutions  offered  by  exist¬ 
ing  methods  are  based  on  heuristics  only,  we  propose  a  method  which  com¬ 
bines  heuristics  with  the  statistical  behavior  of  the  lists  and  which  finds  a  sub- 
optimal  decision  criteria  for  this  problem.  Assuming  that  a  similarity  measure  is 
defined,  we  find  an  optimal  selection  method  (which  identifies  if  two  records 
are  the  same),  based  on  this  measure. 


Takagi,  H.  and  L.  Kleinrock,  “Output  Processes  io  Contention  Packet  Broadcasting  Sys¬ 
tems,”  Computer  Science  Department  Report  No.  CSD-821006,  July  1982. 

The  processes  consisting  of  the  packet  interdeparture  times  in  contention-type 
packet  broadcasting  systems  are  studied.  The  channel  access  protocols  con¬ 
sidered  include  slotted  and  unslotted  ALOHA  and  carrier-sense-multiple-access 
(CSMA)  with  collision  detection  or  with  delay  capture  effect.  Through  analysis 
of  the  channel  activity  cycle,  the  distribution,  mean  and  coefficient  of  variation 
of  the  packet  interdeparture  times  are  explicitly  derived.  Taking  the  reciprocal 
of  the  mean  interdeparture  time,  we  obtain  the  channel  throughput.  Cases  with 
dissimilar  users  are  also  considered.  Application  of  the  present  results  to  the 
packet  queueing  processes  is  suggested. 


Kleinrock,  L.,  “A  Decade  of  Network  Development,”  Journal  of  Telecommunications  Network, 

Spring  1982,  pp.  1-11. 

The  network  technologies  developed  over  the  past  decade,  due  largely  to 
incredible  advances  in  integrated  chip  production,  have  set  the  stage  for  an 
enormous  DP  revolution  the  impact  of  which  will  soon  be  felt  in  the  business 
community  and  in  the  home  and  consumer  markets.  We  have  been  able  to 
identify  those  fundamental  principles  which  render  this  new  technology  so  cost 
effective.  End-to-end  communication  services  are  now  coming  into  place  which 
will  pass  data,  voice,  video,  fax,  graphics,  etc.  from  an  individual  terminal  in  an 
office  or  a  living  room  across  vast  distances  to  remote  processing  facilities.  The 
technological  issues  have  been  addressed,  and  in  this  paper  we  describe  some  of 
the  successes  and  some  of  the  remaining  problems  to  be  solved.  As  engineers, 
we  cannot  restrict  our  efforts  to  the  purely  technical  (and  “nice”)  problems, 
but  rather  we  must  accept  the  further  responsibility  of  familiarizing  ourselves 
with  the  overall  environment  in  which  our  “products”  must  perform;  if  we  fail 
in  this  regard,  the  world  of  commerce,  trade  and  business  will  ignore  our  best 
efforts  as  irrelevant  to  its  needs. 


Nelson,  R.,  “Channel  Access  Protocols  for  Multi-Hop  Broadcast  Packet  Radio  Networks,” 
Ph.D.  Dissertation,  Computer  Science  Department,  University  of  California,  Los  Angeles, 
July  1982,  Report  No.  CSD-820731. 
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ABSTRACT 


A  multi-hop  packet  radio  network  consists  of  a  set  of  geographically  distributed  nodes 
(e.g.  computers,  terminals,  etc.)  which  communicate  using  radio  transceivers  over  a  shared, 
broadcast  radio  channel.  In  this  dissertation  we  study  algorithms,  called  channel  access  protocols. 
that  are  executed  by  the  nodes  of  the  network  to  schedule  their  transmissions  on  the  channel. 
In  particular,  we  analyze  the  performance  obtained  in  mobile  networks  where  nodes  use  the 
slotted-ALOHA  protocol  and  determine  network  parameters  that  maximize  throughput.  These 
results  are  extended  when  we  obtain  a  tight  upper  bound  for  the  performance  of  all  protocols  in 
a  mobile  environment.  This  bound  can  be  used  as  a  standard  by  which  we  compare  the  perfor¬ 
mance  of  other  protocols  in  multi-hop  mobile  networks  and  we  show  that  slotted-ALOHA  is 
about  36%  efficient  in  such  environments. 

In  a  related  area,  we  define  a  new  protocol,  spatial-TDMA.  designed  expressly  for  net¬ 
works  in  which  nodes  are  stationary.  This  protocol  is  collision-free  and  we  find  the  mean  delay 
for  messages  passing  through  nodes  using  the  protocol.  This  leads  us  to  formulate  and  solve  a 
capacity  assignment  problem  for  networks  using  spatial-TDMA,  which  minimizes  the  delay 
averaged  over  all  messages  in  the  network. 

Another  protocol  we  define,  rude-CSMA.  uses  information  only  about  a  node's  local 
environment.  This  protocol  is  called  "rude”  because,  under  certain  circumstances,  the 
throughput  of  the  network  is  increased  if  nodes,  after  sensing  a  busy  channel,  transmit  packets 
with  a  non-zero  rate.  We  derive  the  performance  equations  for  this  protocol  and  then  find  sys¬ 
tem  parameters  that  maximize  throughput  for  different  network  topologies. 
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LIST  OF  NOTATION 


The  following  is  a  list  of  notations  used  in  the  thesis.  Notations  are  consistent 
throughout  and  we  have  organized  the  list  according  to  the  chapter  in  which  the  notation  first 
appears. 
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Mean  density  of  terminals  on  the  plane. 

Transmission  radius  of  nodes  in  the  network. 

The  average  number  of  neighbors. 

The  capture  parameter  of  the  network. 

The  number  of  nodes  in  the  network. 

The  probability  of  transmitting. 

The  effective  number  of  neighbors. 

The  expected  number  of  successes  in  the  network. 
The  offered  load  of  packets  from  a  circle  of  radius  R . 
The  effective  offered  load. 


CHAPTER  3 

R0  The  radius  of  the  entire  network. 

R  i  The  random  variable  for  the  euclidian  distance  for  nodes  one  hop  away. 

R 2  The  random  variable  for  the  euclidian  distance  for  nodes  two  hops 

away. 

R 3  The  random  variable  for  the  euclidian  distance  for  nodes  three  hops 

away. 

0  The  random  variable  for  the  angle  between  the  vectors  of  the  first  and 

second  hop. 

/‘(A)  The  maximum  fraction  of  successful  transmissions. 
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The  size  of  the  time  frame. 

Clique  i. 

The  total  time  from  the  time  frame  allocated  to  C . 

The  time  vector. 
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The  set  of  all  cliques.  C  -  {  Cj,  C2.  •  •  •  ,  Ck  1. 

The  average  flow  over  a  time  frame  of  messages  from  node  /  to  node  /. 
The  total  flow  of  messages  in  the  network  over  a  time  frame. 

The  average  flow  of  messages  over  a  time  frame  over  arc  /. 

The  external  arrival  time  allocated  to  a  specified  queue  over  one  time 
frame. 

The  internal  arrival  time  allocated  to  a  specified  queue  over  one  time 
frame. 

The  service  time  allocated  to  a  specified  queue  over  one  time  frame. 

The  average  number  of  external  arrivals  to  a  specified  queue  over  one 
time  frame. 

The  average  number  of  internal  arrivals  to  a  specified  queue  over  one 
time  frame. 

The  rate  used  in  the  fluid  approximation  over  external  arrival  periods 
for  a  specified  queue. 

The  rate  used  in  the  fluid  approximation  over  internal  arrival  periods 
for  a  specified  queue. 

The  rate  used  in  the  fluid  approximation  over  service  periods  for  a 
specified  queue. 

T„  corresponding  to  queue  i. 

T,n  corresponding  to  queue  i. 

T,  corresponding  to  queue  /. 

L,x  corresponding  to  queue  /. 

L  corresponding  to  queue  /. 

The  average  delay  of  messages  passing  through  queue  /  for  time  vector 
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Parameter  to  rude-CSMA  protocol. 

Parameter  to  rude-CSMA  protocol. 

Average  rate  of  message  arrival  to  nodes  in  rude-CSMA. 

Average  transmission  time. 

Equal  to  yo/fi. 

The  state  of  node  /,  s,(S)  =  1  means  node  /  is  transmitting,  otherwise  it 
is  not  transmitting. 

The  state  of  the  system,  S  -  (s^S),  sy(S),  ■  ■  ■  ,s„(S)). 

The  number  of  neighbors  of  node  /',  in  state  S,  that  are  transmitters. 
The  number  of  neighbors  of  node  /,  in  state  S,  that  are  not 
transmitters. 

The  number  of  adjacent  transmitters  in  state  5. 

The  number  of  adjacent  non- transmitters  in  state  S. 

The  number  of  transmitters  in  state  5. 

The  expected  amount  of  successful  transmissions  in  state  S. 

The  rate  at  which  node  i  presents  packets  to  the  channel  in  state  S. 

The  rate  at  which  node  i  ends  its  transmission  in  state  S. 

The  probability  of  being  in  state  S  for  parameter  values  x  and  y. 

The  expected  throughput  for  parameter  values  .v  and  y. 

The  effective  number  of  idle  neighbors. 


CHAPTER  1 
INTRODUCTION 


In  this  section  we  will  briefly  recount  some  of  the  events  that  lead  to  the  development 
of  computer  networks.  We  will  then  define  terms  used  to  describe  packet  radio  networks  and 
delineate  the  problems  which  are  addressed  in  this  dissertation. 

1.1  A  Little  History 

The  basic  problems  of  computer  science  are  those  involved  with  sharing  resources.  In 
the  1960’s,  since  computing  machinery  was  expensive,  designers  of  computer  systems  broke 
away  from  the  one-job  one-machine  system  configuration  and  developed  the  policy  of  multi¬ 
programming.  This  decreased  the  cost  of  the  machine  time  dedicated  to  each  job,  but  required 
considerable  effort  to  create  and  implement  operating  system  software.  As  the  complexity  of 
the  system  environment  increased,  so  did  the  intricacies  of  the  operating  system,  and  issues 
concerning  its  functionality  and  implementation  became  major  efforts.  In  some  sense  the 
operating  system  became  another  component  of  the  resource  sharing  problem  it  was  designed 
to  solve.  The  difficulty  associated  with  managing  the  complexity  of  these  systems,  combined 
with  economics,  served  as  a  catalyst  to  °teer  things  away  from  building  larger  systems.  Mass 
fabrication  of  integrated  circuitry  drastically  reduced  the  cost  of  computers  to  the  point  where  it 
was  no  longer  vital  to  operate  systems  at  high  utilizations  to  be  cost  effective,  and  thus  provid¬ 
ing  more  computers,  instead  of  making  them  larger  and  more  powerful,  became  more  feasible. 
The  resource  sharing  problems  that  were  created  with  this  trend  consisted  of  developing 
efficient  methods  to  allow  distributed  users  to  share  many  smaller  machines  rather  than  provid¬ 
ing  them  sufficient  processing  power  from  one  large,  commonly  shared  machine.  Besides  prob¬ 
lems  in  distributing  processing,  this  approach  created  the  need  to  connect  these  machines  with  a 
communications  network.  The  resource  sharing  problems  of  operating  systems  were  thus 
extended  to  include  those  associated  with  computer  networking. 

Experimental  networks,  such  as  the  ARPANET  of  the  department  of  defense 
[Klei76a],  demonstrated  the  feasibility  of  the  network  approach.  This  network  used  point  to 
point  cable  to  link  the  communication  processors  of  the  network  together,  which,  since  the  net¬ 
work  itself  was  a  shared  resource,  became  the  next  part  of  the  resource  sharing  problem.  The 
inflexibility  of  point  to  point  cable  to  changes  in  the  topology  of  the  network  motivated  network 
designers  to  seek  a  more  flexible  transmission  media  for  the  links  of  the  network.  The  success¬ 
ful  implementation  of  the  ALOHA-NET  [Abra70a]  of  the  University  of  Hawaii  showed  that  it 
was  possible  to  use  broadcast  radio  as  this  new  media.  A  new  resource  then,  the  radio 
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spectrum,  was  added  to  the  set  of  resources  which  had  to  be  efficiently  shared. 

In  such  networks,  if  two  or  more  nodes  transmit  simultaneously,  their  signals  collide 
and  neither  signal  is  successfully  received.  Thus  the  resource  sharing  problem  of  scheduling 
transmissions  on  the  channel  to  minimize  the  number  of  collisions  in  the  network  became  an 
important  research  area.  This  brief  historical  account  leads  us  to  the  present  day  and  to  the 
topic  of  this  dissertation. 

Our  major  purpose  in  this  dissertation  is  to  lend  understanding  to  these  scheduling  algo¬ 
rithms,  called  channel  access  protocols,  for  radio  networks.  Before  delineating  the  set  of  prob¬ 
lems  addressed  in  this  dissertation  and  its  major  contributions,  we  must  make  a  few  definitions 
in  the  next  section.  In  the  sections  following  this,  we  will  then  relate  in  a  more  precise 
manner,  the  basic  principles  of  resource  sharing  with  the  results  of  our  investigations. 

1.2  Preliminary  Definitions 

In  this  section  we  briefly  define  the  essential  elements  of  a  packet  radio  network.  A 
packet  radio  unit  is  a  hardware  device  consisting  of  a  communications  section  which  transmits  and 
receives  signals,  and  a  logic  section  which  constructs,  decodes,  and  processes  messages.  The 
unit  of  information  sent  by  the  transmitter  is  called  a  packet.  We  will  distinguish  two  types  of 
radio  units,  regular  packet  radios  and  repeaters  This  classification  depends  on  the  nature  of  the 
programs  the  unit  executes  and  we  allow  units  to  change  from  one  classification  to  the  other. 
Typically  the  function  of  repeaters  is  to  form  communication  links  between  regular  units  which 
are  separated  by  distances  greater  than  their  transmission  range,  and  thus  relay  messages  from 
one  unit  to  another.  To  prevent  repetitious  language,  we  will  use  terminals  and  nodes  to  mean 
regular  packet  radio  units. 

The  topology  of  the  network  is  the  spatial  placement  and  transmission  range  of  the 
nodes  of  the  network.  The  channel  is  the  communications  media  used  by  the  transceivers  and 
in  our  work  will  be  broadcast  radio.  We  say  the  channel  is  idle  with  respect  to  a  certain  node  at 
a  certain  point  in  time,  if  that  unit  does  not  detect  a  signal  on  the  channel  at  that  time,  and  that 
the  channel  is  busy  otherwise. 

Suppose  a  terminal  t\  is  sending  a  message  to  terminal  a  and  that  a  is  within  range  of 
fi's  transmission.  Also  assume  that  rfs  message  is  sent  over  the  time  period  [O.r].  We  say 
that  a  captures  fi's  signal  at  time  r,  (0  <  r  <  t ),  if  a  can  decipher  the  information  in  rfs  signal 
at  r.  We  say  that  fj’s  transmission  is  successful  if  a  captures  t\'s  signal  over  all  r  in  [O.f]  other¬ 
wise  we  say  ft' s  signal  collided  with  other  signals  (possibly  noise)  on  the  channel.  In  a  typical 
radio  environment  there  are  two  parameters  that  govern  the  capture  of  signals,  power  and  time. 
In  power-capture ,  ft's  signal  will  be  received  by  a  in  the  presence  of  a  transmission  by  f:.  if  the 
ratio  of  the  power  of  the  two  signals  received  by  a  is  greater  than  a  certain  number.  fi~].  called 
the  capture-ratio.  If  /3-0  ( non-capture )  a  will  receive  r{s  transmission  only  if  there  are  no 


other  transmitters  within  a 's  hearing  range,  and  if  /3*  1.0  (perfect-capture)  a  will  receive  if  fj  is 
closer  to  a  than  ty.  With  time-capture,  t\'s  signal  is  received  by  a  if  /i  began  transmitting  at 
least  a  certain  time,  the  capture  time,  prior  to  the  beginning  of  ty's  transmission. 

Packets  originating  in  one  node  may  need  to  be  sent  to  several  intermediary  nodes  to 
reach  their  final  destinations.  We  call  each  transmission  in  this  relay  a  hop  and  distinguish 
between  single-hop  and  multi-hop  environments.  A  single-hop  environment  is  said  to  be  line-of- 
sight  if  each  node  detects  transmissions  from  all  the  nodes  in  the  network,  and  said  to  be  hidden 
if  some  nodes  do  not  hear  all  the  nodes  in  the  network.  In  the  hidden  environment  we  assume 
terminals  can  only  send  packets  to  terminals  they  can  hear.  In  the  multi-hop  environment  ter¬ 
minals  send  packets  to  destinations  that  are  outside  their  transmission  ranges  and  a  routing  algo¬ 
rithm  is  needed  to  specify  the  next  node  that  is  to  be  the  recipient  of  the  packet. 

Nodes  execute  an  algorithm  called  the  channel  access  protocol  to  determine  when  to 
transmit  their  packets  and  the  information  used  as  input  to  this  algorithm  is  called  the  policy 
state  information.  If  terminals  use  data  containing  information  about  a  proper  subset  of  the 
nodes  in  the  network,  the  policy  is  said  to  be  a  local  policy,  otherwise  we  say  it  is  a  global  pol¬ 
icy. 


We  will  judge  the  performance  of  the  network  by  its  delay-throughput  characteristics. 
The  throughput  of  the  network  is  the  average  number  of  messages  reaching  their  final  destina¬ 
tions  over  a  unit  time  interval  over  a  specified  region  of  the  network.  In  single-hop  environ¬ 
ments,  this  region  is  defined  to  be  the  entire  network  and  thus  the  throughput  of  the  network  is 
equal  to  the  average  number  of  successful  messages  on  the  channel  per  unit  time.  In  multi-hop 
environments  however,  several  successes  can  occur  simultaneously  in  different  parts  of  the  net¬ 
work,  and  thus  we  calculate  throughput  over  a  region  of  unit  area.  The  maximum  throughput 
of  a  particular  network  for  a  particular  channel  access  protocol  is  called  the  capacity  of  the  net¬ 
work.  The  delay  of  the  network  is  the  average  time  between  message  arrival  at  the  originating 
node  and  delivery  to  the  final  node.  A  certain  policy  will  be  said  to  be  optimal  if  we  can  show 
that  no  other  policy  using  a  subset  of  the  policy  state  information  achieves  a  greater 
throughput. 
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1.3  Statement  of  the  Problem 


The  problem  we  will  study  in  this  dissertation  consists  of  creating  efficient  transmission 
policies  for  multi-hop  packet  radio  broadcast  networks.  Thus  the  shared  resource  is  the  broad¬ 
cast  radio  channel  and  the  objective  is  to  maximize  the  capacity  of  the  channel  over  given  sets 
of  policy  state  information.  To  elucidate  some  of  the  characteristics  of  this  problem,  let  us  first 
discuss  problems  in  developing  protocols  for  the  single-hop  line-of-sight  environment.  In  net¬ 
works  of  this  type  all  nodes  share  identical  information  about  the  status  of  the  channel.  This 
information  is  obtained  by  sensing  the  channel  and  is  exploited,  for  example,  in  the  CSMA 
family  of  protocols  [Toba74a].  In  CSMA,  a  terminal  wishing  to  transmit  a  packet  on  the  chan¬ 
nel,  first  senses  the  channel  to  determine  if  it  is  idle  or  busy.  If  idle,  the  terminal  transmits  its 
packet  immediately,  otherwise  it  refrains  until  some  future  time  and  again  senses  the  channel. 
Collisions  occur  only  if  two  or  more  terminals  sense  the  channel  within  a  propagation  delay 
(defined  to  be  the  average  delay  between  the  time  a  message  is  transmitted  on  the  channel  and 
the  recognition  that  this  terminal  is  transmitting  by  all  the  nodes  in  the  network)  of  each  other, 
and  hearing  it  idle,  all  transmit  their  packets.  If  the  propagation  delay  between  terminals  is 
small,  the  probability  of  this  event  is  also  small  and  collisions  are  infrequent. 

What  other  information  would  be  valuable  in  determining  efficient  transmission  policies 
for  single-hop  systems?  Suppose  the  identities  of  the  busy  terminals  were  known  by  every  ter¬ 
minal  in  the  network.  An  optimal  policy  which  both  maximizes  throughput  and  minimizes 
delay  could  then  be  created  by  giving  busy  terminals  permission  to  transmit  in  some  specified 
order,  say  first-come  first-served.  It  is  well  known  that  the  performance  of  such  classical 
queueing  systems  are  a  lower  bound  for  all  possible  access  schemes.  Such  state  information 
must  be  sent  over  the  channel,  though,  at  a  cost  in  the  delay-throughput  performance  of  the 
system.  Indeed  this  communication  and  the  optimal  use  of  specific  control  information  lie  at 
the  root  of  the  difficulty  of  determining  optimal  policies  for  single-hop  line-of-sight  environ¬ 
ments  [Moll81al.  Thus  we  have  that  to  efficiently  share  the  channel,  one  needs  ways  to  distri¬ 
bute  control  information  and  since  the  channel  is  the  only  means  to  perform  this  communica¬ 
tions,  the  control  information  requires  efficient  methods  to  use  the  channel.  We  are  reminded 
of  such  a  cycle  mentioned  in  the  introduction  in  the  case  of  operating  systems,  where  part  of 
the  structure  it  needed  to  manage  to  effectively  share  computing  resources,  was  the  operating 
system  itself. 

The  complexity  of  the  problem  is  greatly  increased  when  the  hidden  environment  is 
considered  since  each  node  hears  only  a  proper  subset  of  the  nodes  in  the  network,  and  neigh¬ 
boring  nodes  do  not  necessarily  agree  on  the  state  of  the  channel.  If  we  try  transmitting  to 
another  node  in  the  network  after  sensing  the  channel  idle,  as  in  CSMA.  there  is  no  guarantee 
that  the  receiver,  who  hears  a  different  subset  of  nodes  than  the  transmitter,  will  also  sense  an 
idle  channel.  The  set  of  subnetworks  created  by  the  local  information,  in  essence  one  for  each 
terminal  in  the  network,  have  a  high  degree  of  interdependence.  For  any  two  nodes,  there  is  a 
set  of  terminals  that  can  hear  both  of  them,  and  simultaneous  transmission  to  this  interference 
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set  will  result  in  a  collision. 


FIGURE  1.1 
A  Sample  Network. 

For  example,  in  Figure  1.1  transmitters  fj  and  t2  have  interference  set  N,  to  which  they  both 
transmit,  and  sets  Nl{  and  N,2  which  hear  only  their  respective  transmitters.  There  are  three 
distinguishable  cases  for  successful  transmission  on  the  channel.  There  are  no  successes  if  both 
r i  and  t2  send  to  Nr  or  if  both  do  not  transmit.  One  success  will  occur  if  exactly  one  of  the  ter¬ 
minals  ti  transmits  or  if  both  transmit  and  one  terminal  sends  to  Af,  while  the  other  sends  to 
Nt.  Finally  there  are  two  successes  if  both  terminals  send  to  N,x  and  N,r  Coordinating  their 
transmissions  to  maximize  the  number  of  successes,  is  called  the  interference  problem  and  will 
be  a  major  area  of  concern  in  this  dissertation. 

In  the  hidden  terminal  environment  knowing  the  identities  of  busy  terminals  is  not 
sufficient,  as  in  one-hop  networks,  to  generate  optimal  transmission  policies.  One  must  now 
also  know  the  destinations  for  packets  in  each  of  the  queues,  and  the  traffic  characteristics  for 
future  arrivals  of  messages. 

There  are  some  theoretical  results  we  can  prove  for  optimal  transmission  policies.  In 
figure  1.1  let  qt[  and  q,2  be  the  conditional  probabilities  that  when  the  transmitters  transmit  they 
send  packets  to  and  N,2  respectively  and  suppose  transmitters  send  packets  with  probabilities 
P,x  and  P,r  Using  the  above  definitions  we  can  derive  an  expression  for  the  expected  number 
of  successes.  There  will  be  one  success  if  exactly  one  of  t\,  and  t2  transmits,  or  if  both 
transmit,  if  one  of  them  transmits  to  the  interference  set  (and  its  transmission  suffers  a 
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collision)  and  the  other  to  either  set  /Vf]  or  ,V,v  Writing  this  out  using  the  above  notation  we 
have: 


P[  One  Success!  —  P,,(l  —  PJ  +  P,,<  1  -  P,,)  +  Pf|P.;((l- ?,,)?,,) 

Two  successes  can  only  occur  if  both  transmitters  send  packets  to  sets  ,V,t  and  V,,.  The  proba¬ 
bility  for  this  is: 

P(  Two  Successes]  -  Pl(  P,2 qt]  q ,, 

The  expected  number  of  successes  is  then  given  by: 

£  -  P(  One  Success!  +  2  Pi  Two  Successes] 


E  “  P(,  +  Pr2  -  Pf) 


We  wish  to  maximize  this  over  choices  of  (P,t,Pfj).  The  solution  will  occur  in  the  interior  of 
the  feasible  set  0<P,,P,  <  1  only  if  the  Hessian  matrix  (the  matrix  of  second  partial  deriva¬ 


tives  i.e. 


&E 

dp,  dp, 


)  is  negative  semi-definite  lAvri76a).  The  Hessian  however  is 


H 


0 

0 


which  is  clearly  not  negative  semi-definite  since  the  determinant  of  the  second  principle  minor 
is  negative.  We  thus  conclude  that  the  optimal  solution  occurs  at  one  of  the  vertices  of  the 
feasible  set.  There  are  four  such  vertices  and  we  can  eliminate  (0,0)  immediately  since  it  has  a 
throughput  value  of  zero.  This  leaves  three  pairs  to  consider  namely  (0,1),  (1.0),  and  (1,1). 
The  expected  number  of  successes  for  vertex  (1,1)  is  equal  to  2-  iq,  +  q,2)  and  equals  1  for 
(1,0)  or  (0,1).  We  see  then  if  q,x  +  q,2  4 1  the  optimal  policy  is  to  be  rude  and  transmit  packets 
regardless  of  the  state  of  the  channel,  and  if  qt]+qh  >1  the  optimal  policy  lies  on  either  the 
(1,0)  or  (0.1)  vertice.  We  should  note  that  to  make  the  policy  fair,  in  the  sense  that  each  user 
obtains  an  equal  share  of  the  capacity  of  the  channel,  we  would  alternate,  in  some  fashion, 
between  (1.0)  and  (0.1).  This  result  generalizes  to  the  case  of  n  interfering  terminals  where 
the  optimal  policy  lies  on  one  of  the  vertices  of  a  (0. 1 ) "  hypercube.  To  achieve  maximal 
throughput,  it  does  not  make  sense  to  assign  transmission  permission  to  a  node  if  its  transmis¬ 
sion  will  cause  a  collision,  and  thus  the  above  example  implies  that  an  optimal  assignment  of 
transmission  permissions  does  not  allow  collisions  in  the  network.  This  motivates  why  we 
define  a  collision-free  protocol  for  multi-hop  networks  called  spatial-TDMA  in  Chapter  4. 
Another  observation  we  can  make  from  the  above  example,  which  motivates  another  protocol. 
rude-CSMA.  which  we  define  in  Chapter  5,  arises  from  the  fact  that  under  certain  cir¬ 
cumstances  nodes  in  the  above  example  should  transmit  even  if  they  sense  a  busy  channel.  In 
rude-CSMA  nodes  adjust  the  rate  at  which  they  present  packets  to  the  channel  according  to  the 
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activity  of  the  channel  in  their  local  environments  and  do  not  necessarily  refrain  from  transmit¬ 
ting  if  they  sense  a  busy  channel.  The  idea,  as  in  the  example,  is  that  a  message  may  be  des¬ 
tined  towards  a  node  outside  the  range  of  the  other  transmitters  in  the  local  environment  and 
thus  transmitting  will  tend  to  increase  the  throughput  of  the  network. 

We  can  refine  the  above  model  by  first  noting  in  the  above  strategy  that  collisions  do 
sometimes  occur,  and  we  have  assumed  that  terminals  have  packets  to  send  at  ail  times.  If  we 
relax  these  assumptions  we  can  create  a  policy  that  maximizes  the  expected  throughput  over  all 
time  by  analyzing  a  more  complex  model.  Suppose  now  we  know  the  messages  and  their  desti¬ 
nations  for  all  the  terminal’s  queues  as  well  as  the  topology  of  the  network.  Also  suppose  we 
know  the  arrival  rate  of  new  messages  to  each  terminal  and  the  distribution  of  message  destina¬ 
tions.  Then  by  modeling  the  system  as  a  Markov  decision  process  [Mine70a]  we  can  assign 
permissions  to  transmit  to  nodes  in  the  network  and  using  a  dynamic  programming  technique 
called  policy-iteration  lHowa60a]  can  determine  a  global  transmission  policy  that  will  maximize 
the  expected  number  of  successes  over  all  time.  This  decision  policy  is  collision-free  and  is 
easily  implemented  if  information  about  the  queues  of  the  network  were  known  to  all  nodes. 
In  general  however,  dissemination  of  this  information  would  require  considerable  channel  capa¬ 
city  and  thus  tends  to  counteract  the  optimality  of  the  proposed  policy. 


1.4  Previous  Work  in  the  Area 

We  briefly  give  in  this  section  an  overview  of  related  research  in  the  area  of  channel 
access  protocols  and  will  concentrate  on  broadcast  protocols  for  the  single-hop  and  multi-hop 
environments. 

Single-hop  -  There  is  a  vast  literature  for  the  single-hop  line-of-sight  environment  and 
we  will  only  discuss  results  relevant  to  later  sections.  The  most  elementary  contention  protocol 
for  this  environment  was  developed  by  Abramson  [Abra70a]  and  is  called  the  ALOHA  protocol. 
In  this  protocol  a  terminal  transmits  a  message  as  soon  as  it  arrives,  and  messages  transmitted 
concurrently,  collide  and  are  not  received  by  their  intented  nodes.  Collided  packets  are  re¬ 
transmitted  after  nodes  wait  a  randomly  selected  period  of  time.  It  can  be  shown  that  under 
the  assumptions  of  an  infinite  population  of  terminals  offering  Poisson  traffic  at  a  constant  rate, 
the  maximum  throughput  of  the  system  over  a  finite  period  of  time  is  1/2  e.  or  about  18%  of 
the  bandwidth  of  the  channel,  and  that  the  system  will  eventually  form  an  infinite  backlog  of 
messages  waiting  to  be  sent  over  the  channel  [Lam75al.  The  throughput  of  the  system  can  be 
doubled  [Robe72a]  by  using  the  ingenious  device  of  dividing  the  time  axis  into  slots  equal  to 
the  length  of  packets.  This  protocol,  called  slotted-ALOHA .  operates  in  a  similar  fashion  to 
ALOHA  except  that  nodes  are  allowed  to  begin  transmission  only  at  the  beginning  of  time 


slots.  This  system  is  also  unstable  [Lam75a],  but  several  methods  for  stablizing  it  have  been 
proposed. 

Observe  that  in  the  ALOHA  system  radio  units  transmit  even  if  the  channel  just  prior 
to  their  transmission  is  busy,  hence  guaranteeing  a  collision.  Tobagi  [Toba74a]  analyzed  a  class 
of  protocols  called  Carrier  Sense  Multiple  Access  (CSMA)  which  attempt  to  avoid  this  futile  colli¬ 
sion.  In  CSMA  when  a  node  wishes  to  transmit  a  packet,  it  first  senses  the  channel  to  deter¬ 
mine  if  the  channel  is  busy.  If  the  channel  is  idle,  the  node  transmits  the  packet  immediately. 
Collisions  only  occur  if  two  or  more  nodes  sense  the  channel  idle  and  begin  transmitting  within 
one  propagation  time,  r,  of  each  other.  If  the  channel  is  sensed  busy,  the  node  refrains  from 
transmitting  the  message  and  executes  one  of  the  following  variants  of  the  protocol. 

1-persistent  —  The  node  waits  until  the  channel  becomes  idle  and  transmits  with  probability  1. 

p-persistent  —  The  node  waits  until  the  channel  becomes  idle  and  transmits  with  probability  p. 

If  the  channel  remains  idle  for  a  propagation  delay  then  the  node  again  transmits 
its  packet  with  probability  p.  If  the  channel  becomes  busy,  the  node  repeats  the 
protocol  after  waiting  a  randomly  selected  period  of  time. 

non-persistent  —  The  node  reschedules  sensing  the  channel  to  a  random  time  in  the  future. 

The  magnitude  of  the  propagation  delay  r  is  critical  to  the  performance  of  the  protocol. 
If  r  is  small,  the  probability  of  a  collision  is  also  small  since  a  transmitting  terminal  is  heard  by 
all  the  nodes  in  the  network  in  a  very  short  period  of  time.  On  a  satellite  channel  however, 
where  the  propagation  delay  is  large  (on  the  order  of  1/4  second),  a  terminal  transmitting  on 
the  channel  after  hearing  it  idle,  will  not  be  heard  until  1/4  of  a  second  later.  During  this  time 
there  is  a  large  probability  that  other  terminals  will  sense  the  channel,  hear  it  idle,  and  ?,1sn 
transmit  their  packets  on  the  channel  thus  resulting  in  a  collision.  The  throughput  of  .tun 
persistent  CSMA  for  propagation  delays  typically  found  in  land  based  packet  radio  networks  is 
about  87%  of  the  bandwidth  of  the  channel. 

Tobagi  also  studied  the  single-hop  hidden  environment  and  devised  a  protocol.  Busy- 
Tone-Multiple-Access  (Toba75a]  in  which  a  central  station,  assumed  to  be  within  line-of-sight  of 
all  terminals  in  the  network,  emits  a  signal  indicating  if  the  channel  is  busy  or  idle.  In  this  way 
all  terminals  by  listening  to  the  busy  tone  share  the  same  channel  status  information  and  can 
use  this  information  to  control  their  transmissions.  This  is  actually  a  simplification  of  Tobagi's 
model  since  he  considered  problems  associated  with  correctly  detecting  the  busy  tone  in  the 
presence  of  noise,  and  therefore  nodes  do  not  always  have  the  same  channel  status  information. 
It  should  be  noted,  that  operating  the  CSMA  protocol  in  the  hidden  environment  without  the 
busy  tone  concept,  results  is  severely  decreased  throughput  since  nodes  sensing  an  idle  channel 
will  transmit,  and  have  the  possibility  of  colliding  another  transmission  already  in  progress.  W’e 


will  derive  equations  for  the  maximum  throughput  of  CSV1A  in  the  multi-hop  environment  in 
Chapter  3  of  this  dissertation. 

Multi-hop  -  Silvester  (Silv80a]  considered  a  packet  radio  network  consisting  of  ran¬ 
domly  placed  terminals  on  the  plane,  with  finite  transmission  radii,  using  slotted-ALOHA.  His 
derivation  assumed  the  existence  of  a  routing  algorithm  which  would  send,  for  each  transmis¬ 
sion,  packets  to  terminals  that  were  closest  to  the  final  destination.  His  results  showed  that  to 
maximize  the  throughput  of  the  system,  terminals  should  transmit  with  a  radius  such  that  their 
transmissions  would,  on  the  average,  be  heard  by  approximately  six  other  nodes.  Because  of 
this  finite  range,  there  could  be  several  simultaneous  successful  transmissions  in  the  network 
during  any  given  slot,  and  this  phenomena  was  named  spatial  re-use  of  the  channel.  In  Chapter 
2  we  generalize  this  model  to  include  capture  and  find  optimal  operating  points  for  the  network. 

Akavia  [Akav78a]  studied  a  multi-hop  model  in  which  a  continuum  of  terminals  used 
slotted  ALOHA  to  send  packets  to  neighboring  terminals.  Transmitters  in  his  model  are 
assumed  to  have  perfect  control  over  the  range  of  their  transmission  and  receivers  are  assumed 
to  be  equipped  with  perfect  capture.  Let  us  classify  the  traffic  into  a  node  to  help  us  explain 
some  of  his  results.  Define  N  to  be  the  average  number  of  packets  in  the  queue  of  a  terminal. 
T  to  be  the  average  time  packets  spent  in  the  queue,  and  X  to  be  the  average  arrival  rate  of 
packets  to  the  queue.  By  Little’s  result  [Litt61a]  we  know  that  N  ~  k  T.  We  say  that  the 
traffic  into  a  terminal  is  bursty  if  N  <  <  1  and  is  steady  otherwise.  By  using  an  ad-hoc  formula 
for  the  delay  of  ALOHA  networks,  Akavia  showed  that  if  traffic  into  a  terminal  was  bursty,  ter¬ 
minals  should  transmit  with  enough  power  to  reach  their  destinations  in  one  hop.  On  the  other 
hand,  for  steady  traffic  he  found  there  is  an  optimal  transmission  range  that  minimizes  the 
delay  of  the  system,  and  in  this  case  packets  require  several  hops  before  reaching  their  final 
destinations.  In  another  approach  throughput  was  increased  in  the  single-hop  environment  by 
dividing  terminals  into  groups  with  different  transmission  powers  and  assuming  the  transmis¬ 
sions  from  low  power  groups  do  not  interfere  with  those  of  higher  power.  It  can  be  shown  that 
the  throughput  of  the  ALOHA  channel  increases  and  as  the  limit  of  the  number  of  power 
groups  goes  to  infinity,  it  is  shown  the  throughput  of  the  system  goes  to  one. 

Tobagi  considered  the  2-hop  environment  in  a  series  of  papers  [Toba80a.  Toba80b] 
where  he  analyzed  two  topologies  for  different  access  schemes.  In  both  topologies,  a  set  of  N 
repeaters  are  assumed  to  receive  packets  from  N  independent  sets  of  packet  radio  networks 
each  consisting  of  an  infinite  number  of  terminals  offering  traffic  at  a  Poisson  rate.  Repeaters 
send  this  traffic  to  a  central  station  CS,  and  only  terminal-to-CS  traffic  is  considered.  The  two 
topologies  arise  in  the  manner  in  which  repeaters  are  assumed  to  be  connected  to  the  CS.  In 
the  Star  topology,  repeaters  are  assumed  to  have  a  link  to  the  CS  but  are  not  connected  to  any 
other  repeater,  and  in  the  Fully  Connected  (FC)  topology  all  repeaters  and  the  CS  are  assumed 
to  be  mutually  connected.  Repeaters  initially  are  assumed  to  have  a  single  buffer  for  storing 
packets. 


In  the  first  model,  terminals  use  slotted-ALOHA  to  send  packets  to  their  repeaters. 
Upon  successful  reception  of  a  packet,  a  repeater  forwards  the  packet  to  the  CS  according  to 
one  of  two  protocols.  The  first  protocol.  Delay  First  Transmission  (DFT).  operates  by  having 
repeaters  with  packets  to  send,  transmitting  them  with  probability  p  in  each  slot.  In  the  second 
protocol.  Immediate  First  Transmission  <IFT),  nodes  transmit  newly  arrived  packets  with  pro¬ 
bability  1  in  the  next  slot,  and  collided  packets  are  re-transmitted  in  the  following  slots  with 
probability  p.  This  variation,  which  avoids  the  initial  packet  delay,  would  intuitively  decrease 
the  overall  packet  delay  from  the  repeater  to  the  CS  in  low  traffic  conditions.  This  is  verified  in 
Tobagi’s  paper.  Tobagi  observed  that  although  increasing  the  buffer  size  from  one  to  two 
packets  did  improve  the  delay  of  the  system,  further  increases  in  buffer  size  had  only  a  minor 
effect.  The  system  is  said  to  be  channel  bound  rather  than  storage  bound. 

In  the  second  model,  the  FC  topology  is  studied  where  terminals  use  non-persistent 
CSMA  and  the  repeaters  use  the  IFT  protocol.  Because  of  the  increased  throughput  of  CSMA 
over  slotted-ALOHA,  the  throughput  of  the  two-hop  system  is  also  increased.  Since  successful 
reception  of  a  packet  by  a  repeater,  say  a,  implies  in  the  FC  case  that  no  other  repeater  is 
transmitting,  the  channel  in  the  next  slot  will  be  idle  if  a  does  not  transmit  its  received  packet. 
With  IFT  however,  a  does  transmit  thus  taking  advantage  of  this  information,  and  a  collision 
will  occur  only  if  at  least  one  other  repeater  received  a  packet  simultaneous  to  a.  The  probabil¬ 
ity  of  this  is  low  and  it  is  shown  that  the  delay  suffered  by  a  packet  in  going  from  a  repeater  to 
the  CS  is  about  one  packet  transmission  time.  Once  again  the  system  is  not  storage  bound  and 
increasing  the  buffer  size  has  little  effect  on  the  performance  of  the  system. 

Capture  -  In  the  paper  by  Roberts  [Robe72bl,  power-capture  is  defined  and  the 
throughput  for  a  single-hop  slotted  ALOHA  line-of  sight  system  is  derived.  It  is  shown  that 
capture  increases  the  usable  bandwidth  of  the  channel  and  delay  equations  are  derived  for  satel¬ 
lite  and  ground  packet  radio  systems.  It  is  also  shown  that  as  the  capture  ratio  varies  from  no¬ 
capture  to  perfect-capture,  the  throughput  of  the  system  increases. 

Fratta  and  Sant  [Frat30a]  considered  a  packet  radio  network  with  repeaters  and  regular 
packet  radio  units  distributed  randomly  on  the  plane.  Using  a  power-capture  model  they 
derived  lower  and  upper  bounds  on  the  throughput  per  unit  area  of  the  system  as  a  function  of 
the  traffic  offered  per  unit  area.  It  is  shown  that  the  bounds  for  this  throughput  are  a  function 
of  the  traffic  per  area  and  the  density  of  the  repeaters.  Traffic  is  assumed  to  originate  at  the 
nodes  and  after  being  transmitted  to  the  repeaters,  is  sent,  by  some  unspecified  means  to  a  cen¬ 
tral  station.  In  a  second  model,  a  network  of  randomly  placed  nodes  is  studied  for  throughput 
characteristics  in  the  single-hop  line-of-sight  environment.  Under  the  assumption  that  each 
node  transmits  with  the  same  probability  in  each  slot,  and  that  they  transmit  to  all  other  nodes 
with  equal  probability,  it  is  shown  that  the  system  has  an  asymptotic  throughput  of  1  packet  per 
slot. 


Davis  and  Gronemeyer  [Davi80al  studied  a  spread  spectrum  system  with  time-capture. 
They  studied  a  slotted  ALOHA  line-of-sight  system,  and  developed  a  clever  way  to  prevent  ter¬ 
minals  close  to  a  given  receiver  from  “stealing”  the  channel  from  terminals  further  away.  This 
stealing  would  occur  because  their  proximity  to  the  receiver  would  give  them  an  unfair  advan¬ 
tage  of  being  received  first,  and  thus  be  more  likely  to  be  captured  over  more  distant  nodes. 
By  causing  terminals  to  delay  transmitting  for  a  period  of  time  linearly  dependent  upon  their 
distance  to  the  receiver,  and  by  adding  to  this  a  randomly  selected  interval  of  time,  it  is  shown 
that  all  nodes  will  have  transmission  times  uniformly  distributed  over  some  interval.  Thus  each 
node,  regardless  of  its  distance  to  the  receiver,  has  an  equal  chance  of  being  the  first  to  be 
received. 


1.5  Contributions  of  this  Dissertation 

In  this  section  we  outline  the  main  contributions  of  this  dissertation.  In  Chapter  2  we 
analyze  the  set  of  packet  radio  networks  in  which  nodes  are  randomly  distributed  on  the  plane 
to  determine  throughput  characteristics  of  the  network  as  a  function  of  transmission  probability 
and  range  of  transmission.  This  network  topology  corresponds  to  a  mobile  packet  radio  net¬ 
work  and  the  shared  resource  we  consider  in  this  chapter  is  the  spatial-capacity  of  the  radio 
channel.  In  our  model  of  the  network,  nodes  use  slotted-ALOHA  to  access  the  channel  and 
have  the  ability  to  capture  signals.  We  derive  the  following  equation  for  the  throughput  of  the 
network: 

y(fi,N,p,)  -  (1  —  /?)(1  —  e~sl1)  p  e~Sp  Q 
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where  Q  ^  ^ y  (1  —  0%/0),  n  is  the  number  of  nodes  in  the  network,  N  is 

the  average  number  of  neighboring  nodes  that  are  within  the  transmission  range  of  a  randomly 
selected  node,  p  is  the  probability  of  transmitting  in  any  random  slot,  and  0  is  the  capture 
parameter  of  the  receivers  of  the  network  (recall  a  signal  is  captured  if  the  ratio  of  its  power  to 
that  of  the  other  signals  on  the  channel  is  greater  than  0-1  and  that  0—0  is  non-capture  and 
0-1  is  perfect-capture).  We  can  see  that  the  throughput  of  the  network  increases  as  \fn  and 
this  factor  represents  the  spatial-capacity  of  the  channel.  This  dependency  results  from  the  fact 
that,  because  a  transmission  is  heard  by  only  a  small  subset  of  the  nodes  in  the  network, 
several  transmissions  can  be  simultaneously  successful  in  any  given  slot.  Typical  receivers  in 
such  networks  have  a  capture  parameter  of  0-0.7.  For  this  value,  the  point  on  the  LV.p) 
plane  that  optimizes  throughput  can  be  found  to  be  .V-4.99  and  p=. 21  at  which  point  the 
throughput  of  the  network  is  given  by  .OliVn  .  If  we  define  the  offered  load.  G,  of  the  system 


to  be  the  average  number  of  packets  presented  to  the  channel  from  an  area  equal  to  that  of  a 
transmission  radius,  then  our  results  show  that  (7s  1  optimizes  performance.  This  is  a  generali¬ 
zation  of  the  same  result  for  the  single-hop  environment.  In  essence  this  model  corresponds  to 
finding  optimal  solutions  to  the  interference  problem  using  a  minimal  set  of  information. 

We  extend  these  results  in  the  third  chapter  by  determining  an  upper  bound  for  the 
spatial-capacity  of  random  networks  over  all  protocols  in  a  mobile  environment.  The  goal  is  to 
free  the  model  of  the  second  chapter  from  the  slotted-ALOHA  assumption  to  determine  which 
insights  obtained  from  that  model  generalize.  In  particular,  we  find  that  the  maximum 
expected  fraction  of  successful  transmissions  in  a  connected  random  planar  packet  radio  net¬ 
work,  over  all  protocols,  is  upper  bounded  by: 

f'(N)  -  .9278//V 

where,  again,  jV  is  the  average  number  of  neighboring  nodes  that  are  within  transmission  range 
of  a  randomly  selected  node.  This  is  a  good  rule  of  thumb  to  use  in  judging  the  performance 
of  other  protocols  with  these  network  assumptions.  When  we  evaluate  the  performance  of  our 
results  for  optimal  slotted-ALOHA  against  this  upper  bound,  we  find  that  this  protocol  achieves 
about  Me  of  the  performance  given  by  the  above  equation.  Recall  that  this  is  equal  to  the 
efficiency  of  slotted-ALOHA  in  a  single-hop  environment,  and  thus  we  have  the  intuitively 
pleasing  result  that  the  maximal  capacity  of  slotted-ALOHA  is  the  same  in  single-hop  and 
multi-hop  environments. 

In  the  fourth  chapter  we  seek  to  find  the  average  delay  messages  encounter  in  passing 
through  a  network  using  a  collision-free  protocol.  This  causes  us  to  define  a  new  channel 
access  protocol,  called  spatial-TDMA,  for  multi-hop  environments  in  which  nodes  are  assumed 
to  be  stationary.  This  corresponds  to  solving  the  interference  problem  with  full  global  informa¬ 
tion.  One  advantage  of  using  broadcast  radio  in  such  an  environment  lies  in  the  fact  that  the 
capacity  of  the  links  in  such  a  network  can  be  dynamically  allocated  in  an  adaptive  manner. 
Allocating  this  capacity  is  a  resource  sharing  problem  which  we  also  address  in  chapter  four.  In 
particular  we  determine  an  approximation  to  the  average  delay  that  messages  experience  in 
passing  through  the  nodes  of  a  network  using  such  a  protocol.  We  determine  that  the  delay  is 
highly  dependent  upon  the  total  time  in  a  cyclic  time  frame  during  which  specified  nodes  in  the 
network  are  allowed  to  transmit,  and  also  upon  the  ordering  of  such  transmissions.  This  leads 
us  to  define  a  capacity  assignment  problem  which  attempts  to  minimize  the  average  message 
delay  for  all  messages  in  the  network  over  all  possible  orderings  and  lengths  of  transmission 
periods.  Although  this  problem  is  mathematically  intractable,  we  do  find  such  optimal  operat¬ 
ing  points  for  networks  in  which  the  ordering  of  the  intervals  has  been  randomly  assigned. 
This  is  a  upper  bound  on  the  delay  that  would  be  obtained  for  the  best  possible  ordering. 
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In  the  fifth  chapter  we  seek  to  determine  optimal  use  of  the  information  that  nodes 
obtain  when  they  sense  the  channel  in  their  local  environments.  This  causes  us  to  define  a  new 
class  of  protocols  called  rude-CSMA  which  use  this  information.  We  formulate  the  protocol, 
derive  equations  for  its  performance,  and  optimize  it  over  different  classes  of  networks.  This 
investigation  corresponds  to  solving  the  interference  problem  using  only  local  information.  In 
these  protocols  we  assume  that  nodes,  after  sensing  the  channel,  can  estimate  the  number  of 
transmitting  nodes  that  are  within  their  hearing  distance.  Intuitively,  if  this  number  is  large, 
there  is  high  probability,  if  the  sensing  node  transmits  a  packet,  that  its  message  would  undergo 
a  collision.  If  however,  there  was  only  one  transmitter,  say,  in  its  local  neighborhood,  then 
there  exists  the  possibility  that  the  sensing  node  would  be  sending  a  packet  to  a  radio  unit  that 
was  outside  the  transmitting  node’s  range,  and  it  would  be  successfully  received.  Thus 
transmitting  with  a  non-zero  rate  might  increase  the  number  of  simultaneous  successful 
transmissions  in  the  network.  This  motivates  the  protocol  and  explains  why  we  have  called  it 
“rude-CSMA”,  since  nodes  sometimes  transmit  even  when  they  sense  the  channel  busy  If  we 
define  the  state  S  —  (s|(5),  s2(S),  •  •  •  ,s„(S))  of  the  system  to  be  a  binary  vector  in  which 
s,(S)— 1  if  node  /  is  transmitting  a  packet  and  equal  to  0  otherwise,  then  the  rate  at  which  node 
/  presents  packets  to  the  channel  (assuming  that  s,(S)-0)  in  state  S  is  given  by: 


r'o(S) 


yox 


Vt  (S> 

y 


where  y0  is  the  arrival  rate  of  messages  to  each  node  of  the  network,  Nq(S)  is  the  number  of 
neighboring  nodes  of  node  i  that  are  not  transmitting,  \\  ( S )  is  the  number  of  transmitting 
neighboring  nodes,  and  x  and  y  are  parameters  of  the  protocol.  On  the  other  hand  if  node  i  is 
a  transmitter,  s,(S)-l,  then  the  rate  at  which  its  transmission  stops  is  given  by. 

r\  ( S )  —  n 


which  corresponds  to  exponential  packet  lengths.  One  sees  that  the  protocols  defined  by  these 
rates  have  ALOHA  Or— 1,  >-l)  and  CSMA  Or-1,  .y-0)  as  special  cases  and  thus  should  have 
an  optimal  performance  as  good  as  either  of  them.  It  can  be  shown  that  the  steady  state  proba¬ 
bility  of  state  S  is  given  by: 

IKS)  -  C/ISI 

where  C  is  a  normalization  constant,  M(S)  is  the  number  of  transmitters.  Bo(S)  is  the  number 
of  neighboring  nodes  that  are  not  transmitters  and  S|(S)  are  the  number  of  neighboring 
transmitters  for  state  S.  Using  these  equations  we  analyze  specific  networks  to  determine  the 
parameter  values  Or,  y)  that  achieve  maximal  throughput.  In  particular,  we  find  that  over  all 
random  topologies  that  we  studied,  y— 0  had  maximal  performance  and  thus  nodes  should 
engage  in  a  type  of  CSMA  protocol  (never  transmit  when  the  channel  is  sensed  busy).  For  lat¬ 
tice  type  networks  we  found  that,  for  certain  input  values,  throughput  was  increased  if  nodes 
after  sensing  a  busy  channel,  transmit  with  a  non-zero  rate  (i.e.  y^O).  It  thus  turns  out  that, 
in  real  networks  (which  would  have  topologies  more  random-like),  the  binary  information 
obtained  from  sensing  the  channel  is  sufficient  to  determine  if  one  should  transmit  or  not. 


13 


We  have  adopted  what  we  hope  will  be  a  clear  writing  style  in  this  dissertation 
Chapters  have  been  made  as  self-sufficient  as  possible  to  allow  the  reader  to  selectively  read 
portions  they  find  most  interesting  without  the  necessity  of  reading  all  the  previous  material. 
This  of  course,  implies  that  there  is  a  certain  amount  of  redundancy  and  repeated  definitions  in 
the  text,  but  it  is  hoped  that  these  repetitions  will  serve  to  clarify  the  material  rather  than  to 
bore  the  reader.  A  list  of  symbols  and  notations  which  are  used  throughout  the  dissertation  is 
included  in  the  beginning  pages  and  all  notations  used  in  a  chapter  are  defined  in  the  introduc¬ 
tory  sections  of  the  chapter. 


1.6  Applications  of  Packet  Radio  Networks 

The  main  advantages  of  packet  radio  networks  over  conventional  networks  is  that  they 
are  not  dependent  on  fixed  topologies,  are  easy  to  establish,  and  can  operate  unattended.  These 
characteristics  allow  terminals  to  be  mobile  and  to  be  attached  to  a  diversity  of  computing  and 
sensor  devices.  Many  interesting  new  applications  of  this  technology  can  be  found  in  [Lick78a], 
In  one  such  application,  sensors  are  carried  by  elderly  persons  which  monitor  their  physical 
orientation  while  walking.  If  the  sensor  detects  a  horizontal  position  (implying  they  fell  down) 
it  notifies  the  client’s  home  computer  which  telephones  a  hospital  with  a  preprogrammed  mes¬ 
sage  and  the  location  (as  determined  by  the  local  radio  network)  of  the  client.  When  medics 
arrive,  they  can  solicit  the  help  of  a  physician  located,  miles  away,  by  attaching  sensors  to  the 
client’s  body  that  transmit  data  the  doctor  can  use  to  control  and  monitor  the  treatment  given. 
Any  number  of  other  applications  using  unattended  sensors  that  communicate  with  computers 
to  control  and  monitor  specific  activities  such  as  alarm  systems  for  security,  fire,  and  gas  or 
water  leakage,  are  natural  applications  of  such  networks. 

The  fact  that  it  is  not  necessary  to  use  wires  to  connect  devices  together  allows  easier 
implementation  of  communication  systems.  One  can  imagine  a  hierarchical  network  in  which 
the  lowest  level  is  a  local  community  radio  network  that  transmits  packets  to  a  central  subnet¬ 
work  station.  These  packets  would  then  be  relayed  along  intra-city  high  bandwidth  lines  to  a 
main  switching  station  that  connects  with  a  long  haul  inter-state  network  having  links  to  inter¬ 
national  satellite  facilities.  Since  radio  is  used  on  the  local  level,  there  is  no  necessity  tor  erect¬ 
ing  poles  and  stringing  cable  directly  to  subscriber  s  equipment.  Besides  beautifying  the 
landscape  of  our  cities,  this  also  has  the  advantage  that  providing  service  to  a  new  subscriber  is 
as  easy  as  outfitting  them  with  proper  radio  equipment.  Both  voice  and  data  would  pass  digi¬ 
tally  over  this  network  thus  supporting  a  mobile  telephone  system  and  providing  access  to  a 
wealth  of  data  and  computing  resources  resident  on  the  net. 


Since  packet  radios  will  be  inexpensive  and  about  the  size  of  a  pocket  calculator 
[Kahn78a],  we  would  imagine  their  use  will  be  widespread  and  thus  revolutionize  the  manner 
in  which  the  society  negotiates  its  information  transactions.  No  longer  will  we  have  to  physi¬ 
cally  wait  in  lines  to  negotiate  our  banking  business,  for  instance,  when  such  transactions  could 
be  performed  from  the  convenience  of  our  homes  (or  even  our  cars!)  and  some  information  in 
the  present  form  of  newspapers,  magazines,  books,  stereo  records,  photographs  and  video  could 
be  sent  digitally  over  the  network,  stored  on  disk,  and  recalled  at  our  convenience. 

In  emergency  situations,  say  as  in  an  earthquake  disaster,  the  flexible  topology  of 
packet  radio  networks  is  necessary  to  establish  communications,  since  in  such  situations  land 
based  lines  will  likely  be  destroyed.  The  immediate  adaptation  of  such  networks  to  changing 
topologies  is  critical  to  the  performance  of  such  an  operation  in  which  new  rescue  teams  con¬ 
tinually  arrive  with  more  packet  radio  units;  and  the  fact  that  it  is  easy  to  change  the  responsi¬ 
bilities  of  a  radio  unit  from  a  regular  packet  radio  unit  to  an  unattended  repeater,  allows  a 
multi-hop  network  to  be  configured  rapidly  thus  increasing  the  useful  range  of  these  networks. 
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CHAPTER  2 

OPTIMAL  THROUGHPUT  FOR  SLOTTED-ALOHA  IN  A  MOBILE  NETWORK 


In  this  chapter  we  derive  the  throughput  equation  for  a  mobile  packet  radio  network  in 
which  nodes  use  siotted-ALOHA  to  access  the  channel.  We  then  optimize  this  equation  over 
certain  parameter  values  to  find  the  optimal  throughput  for  these  networks.  This  corresponds 
to  solving  the  interference  problem  using  a  minimal  amount  of  information. 


2.1  Introduction 

In  this  chapter  we  analyze  a  basic  trade-off  in  multi-hop  radio  networks  which  concerns 
the  throughput  of  the  network.  Recall  that  throughput  is  defined  to  be  the  number  of  packets 
reaching  their  final  destinations  per  unit  time.  To  illustrate  the  problem  discussed  in  this 
chapter  we  will  consider  a  network  under  two  different  cases.  First  suppose  that  all  nodes  of 
the  network  transmit  with  a  very  small  power.  This  implies  that  packets  must  be  relayed  over 
several  hops  before  reaching  their  final  destination  and  thus  tends  to  decrease  the  throughput  of 
the  network.  We  cannot  however  conclude  that  increasing  the  range  of  the  transmitters  of  the 
network  will  increase  the  throughput  of  the  network.  For  suppose  nodes  do  have  an  increased 
range  of  transmissions.  Then  messages  must  be  relayed  fewer  times,  but  at  each  hop  there  is  a 
higher  probability  that  the  transmission  is  not  successfully  received  due  to  the  increased 
amount  of  interference  caused  by  each  transmission.  This  then  increases  the  number  of  colli¬ 
sions  in  the  network  and  tends  to  decrease  the  throughput. 

The  tradeoff  can  thus  be  summarized  as:  increasing  the  power  at  which  a  node 
transmits  decreases  the  number  of  hops  needed  to  reach  the  final  destination  but  causes  more 
interference  at  each  hop,  while  decreasing  the  power  has  the  corresponding  opposite  eflfects. 
Finding  the  operating  point  that  maximizes  the  throughput  for  a  random  network  is  the  prob¬ 
lem  that  is  discussed  in  this  chapter  and  corresponds  to  solving  the  interference  problem  in 
which  nodes  use  a  minimal  set  of  information  (namely  the  density  of  terminals  on  the  plane). 


2.2  The  Models 


Throughout  this  chapter  we  will  make  the  following  assumptions  about  the  network: 

a.  Topology  -  We  assume  that  packet  radios  are  distributed  according  to  a  Poisson 
point  process  on  the  plane  with  a  mean  density  of  A.  packet  radio  units  (also  called 
terminals)  per  unit  area.  We  are  interested  in  finding  the  throughput  for  an  area 
containing  a  large  number,  n ,  of  packet  radios  and  will  ignore  edge  effects.  This 
topology  represents  an  instantaneous  snapshot  of  a  large  mobile  packet  radio  net¬ 
work. 

b.  Stations  -  We  assume  that  each  packet  radio  transmits  with  fixed  power,  and  that 
all  n  stations  (i.e.  terminals)  in  the  network  transmit  with  the  same  power  on  the 
same  frequency  band.  Receivers  are  assumed  to  be  able  to  receive  a  signal  from 
another  station  if  that  station  is  within  a  radius  R  of  the  transmitter,  and  under 
certain  circumstances,  can  successfully  capture  one  of  several  simultaneous 
transmissions  within  its  hearing  range.  Let  us  refine  the  definition  of  capture 
given  in  the  introduction  by  considering  a  receiver  a  which  is  within  range  of  two 
transmitters  t \  and  ti  and  assuming  that  t{  has  a  packet  destined  for  a.  Let  P,  and 
P2  be  the  powers  of  the  signals  received  by  a,  and  r\  and  r2  be  the  distances 
between  a  and  the  two  transmitters.  Whenever  both  r,  and  h_  transmit  their  pack¬ 
ets  simultaneously,  their  signals  interfere  with  each  other.  In  the  absence  of  cap¬ 
ture,  station  r,  will  not  be  received  correctly.  With  capture  however,  station  a  can 
successfully  receive  t\'s  transmission  if  P\/ Pi  >  0~'  (0  <  0  <  1),  where  0  is 
called  the  capture-ratio.  Assuming  omni-directional  antennas  on  the  plane  and 
equal  transmitting  power  for  all  stations,  this  ratio  of  powers  can  be  converted  into 
a  ratio  of  distances  since  the  power  of  a  received  signal  decreases  as  the  inverse 
square  of  the  distance.  Thus,  using  this  distance  measure,  a  will  capture  t\  if 
ri!  r\  >  0-l/2.  From  Figure  2.1  we  see  that  ti  will  be  successful  if  t2  lies  outside 
the  circle  of  radius  r0~1/2,  called  the  capture  radius.  Observe  that  0  —  0  implies 
that  simultaneous  transmission  will  always  cause  a  collision  (non-capture),  and 
0-1  implies  f|  will  be  received  if  it  is  simply  closer  to  a  than  t2  (perfect-capture). 
Well  designed  FM  receivers  have  a  capture  ratio  approximately  equal  to  0.7 
[Robe72b],  Although  the  appearance  of  the  exponent  (-1/2)  of  0  in  the  above 
equation  appears  awkward,  this  selection  simplifies  later  equations.  The  last 

assumption  we  make  about  the  stations  is  that  each  radio  is  always  busy,  and  thus 
we  study  the  heavy  traffic  case. 


FIGURE  2.1. 
Definition  of  Capture. 


Channel  Access  Method  -  We  assume  that  the  time  axis  is  slotted  and  that  the 
probability  a  station  will  transmit  a  packet  on  the  channel  in  any  given  slot  is  p. 

Traffic  Matrix  -  Since  nodes  are  Poissonly  distributed  on  the  plane  with  mean 
density  A,  and  since  any  station  can  send  and  receive  packets  within  a  radius  R. 
every  station  has  on  the  average  N-\nR2  neighbors  (  terminals  within  its  hear¬ 
ing  and  transmitting  range).  We  assume  the  global  traffic  matrix  for  all  the  n 
nodes  in  the  network  is  uniform,  and  thus  the  probability  of  sending  to  any  partic¬ 
ular  node  in  the  network  is  1  In. 

Routing  -  We  choose  to  study  the  case  where  packets  destined  towards  a  particu¬ 
lar  node  F  in  the  network  are  routed  with  equal  probability  towards  one  immediate 
neighboring  node  that  lies  in  the  general  direction  of  F.  For  example  in  Figure  2.3 
there  are  k  terminals  within  transmitting  range  of  t  that  are  closer  to  F  than  to  t. 
Transmitter  t  will  pick  one  terminal  from  these  k  neighbors  with  probability  1  !k. 
Suppose  t  transmits  to  node  i.  We  call  the  difference  of  the  distances  between  t 
and  F ,  and  i  and  F,  the  forward  progress  of  the  transmitted  message. 


Let  us  justify  the  random  routing  assumption  by  comparing  it  to  an  optimal  routing 
model.  In  such  a  routing  policy  packets  are  assumed  to  be  relayed  to  a  neighboring  terminal 
that  is  closest  (furthest  along  the  path)  to  that  packet's  final  destination.  Such  an  optimal  rout¬ 
ing  policy  is  not  realizable  for  a  mobile  packet  radio  network  since  it  requires  the  exact  location 
of  all  terminals  (which  are  assumed  to  be  moving)  as  well  as  that  of  the  final  destination,  and 
would  provide  an  upper  bound  for  network  performance.  This  upper  bound  can  be  easily  calcu¬ 
lated  in  networks  without  capture  because  the  probability  of  being  successfully  received  is 
independent  of  the  distance  between  the  transmitter  and  receiver.  In  the  capture  environment 
however,  nodes  closer  to  the  transmitter  have  a  greater  probability  of  receiving  a  transmitted 
signal  than  those  further  away.  This  non-uniformity  makes  the  calculation  of  the  distance 
covered  in  one  transmission  for  an  optimal  routing  policy  difficult. 


FIGURE  2.2. 

The  Routing  Assumption. 

To  be  specific,  suppose  a  transmitter  /  has  a  message  to  send  to  a  particular  final  node  F  (Figure 
2.2).  Suppose  t  has  k  neighbors  lying  in  the  half  circle  of  his  transmission  radius,  towards  F. 
and  that  their  distances  from  /  are  (rt,  r2 . rk).  Let  (r(,  :2 . rk)  be  the  vector  of  pro¬ 

jected  distances  towards  F,  hence  if  /  transmits  to  node  i  at  r ,  the  progress  towards  F  will  be 
Let  P(r,)  be  the  probability  that  node  /  successfully  receives  t's  transmission.  A  locally  optimal 
routing  algorithm  for  this  system  is  defined  as  one  that  sends  all  packets  tow;,*  r  to  the  node 
j  that  has  the  maximum  value  of  *  Pir).  To  determine  this  value,  one  must  calculate  the 
joint  probability  for  (ri,  r 2 . rk  )  and  (ri.c; . r^l  for  all  k.  to  determine  the  density  for 


the  maximum  projected  distance.  This  would  then  be  unconditioned  on  k  to  determine  the 
density  for  the  maximal  forward  progress.  In  the  non-capture  environment  P\r,)  =  P{r,)  for 
all  i  and  j,  and  thus  maximizing  :,P(r,)  implies  picking  the  maximum  Over  all  sets  of  k 
nodes,  the  probability  that  the  maximum  projected  distance  is  equal  to  a  certain  value,  say  r,  is 
seen  to  be  the  probability  that  there  are  no  terminals  in  the  half  circle  from  t  that  are  closer  to 
F(the  shaded  region  A  in  Figure  2.3). 


F 


FIGURE  2.3. 

Probability  of  Successful  Reception. 

Since  terminals  are  Poissonly  distributed  on  the  plane,  this  probability  equals  In  the  cap¬ 
ture  environment  however,  P(r,)  —  P{r,)  only  if  r,  -  r,,  and  the  above  calculation  is  no  longer 
valid.  The  optimal  routing  algorithm  now  will  no  longer  always  pick  the  node  with  the  max¬ 
imum  r,,  because  the  product  :.P(r,)  may  not  be  maximal.  To  avoid  these  computational 
complexities  and  in  endeavoring  to  create  a  practical  bound  for  packet  radio  networks,  we  have 
chosen  to  assume  a  random  routing  policy  (assumption  e.  above).  Random  policies  have  been 
proposed  for  networks  of  this  kind  [K!ei64a.  Liu80al  and  our  calculations  will  be  a  lower  bound 
on  the  performance  for  algorithms  that  always  send  packets  in  the  direction  of  their  final 
destination. 
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With  this  random  routing  assumption  we  will  analyze  two  models  that  differ  in  the  way 
capture  is  defined  in  relation  to  the  maximum  transmission  distance  R.  Continuing  with  the 
scenario  of  the  previous  section,  it  is  certainly  true  that  in  the  presence  of  multiple  simultane¬ 
ous  transmissions  receiver  a  can  only  successfully  receive  packets  from  its  closest  transmitting 
neighbor.  Thus  we  can  say  that  a  receives  from  t\  if  t]  is  a's  nearest  transmitter,  and  there  are 
no  other  transmitters  within  the  capture  radius.  Let  us  therefore  suppose  in  our  scenario  that  t { 
is  a's  closest  transmitter  and  is  located  a  distance  r  away.  Our  two  models  differ  in  the  way 
they  define  the  capture  radius  in  relation  to  the  maximum  transmitting  distance  R.  Suppose 
rf5~uz  >  R.  or  in  words,  that  the  capture  radius  is  greater  than  a's  maximum  hearing  distance. 
Certainly  a  transmitter  located  further  than  R,  say  at  r’  such  that  R  <  r'  <r  (i~12  will  have  no 
effect  on  a's  reception  since  his  signal  will  be  too  weak  to  be  received.  Thus  in  Model  1  we 
define  the  capture  radius  to  be  equal  to  the  minimum  of  rf3~112  and  R.  The  area  that  must 
contain  no  other  transmitters  for  a  to  successfully  receive  f/s  transmission,  the  dean  area,  for 
Model  1  is  the  annulus  of  inner  radius  equal  to  r  and  outer  radius  equal  to  the  minimum  of 
rp~V2  and  R.  This  definition  for  the  capture  radius  for  Model  1  however  gives  weak  signals 
coming  from  a  transmitter  located  at  a  distance  slightly  less  than  R ,  say  at  R-e.  a  greater  pro¬ 
bability  of  being  successfully  received  than  a  transmitter  with  a  smaller  value  of  r.  since  the 
clean  area,  the  annulus  of  inner  radius  R-e  and  width  e,  is  infinitesimal.  In  practice  however, 
if  there  was  a  transmitter  at  a  distance  slightly  greater  than  R,  say  at  R  +  €,  it  would  disrupt 
reception  since  the  ratio  of  the  powers  of  the  two  transmitters  would  be  close  to  1  even  though 
the  second  transmitter's  signal  was  very  weak.  Model  2  attempts  to  account  for  this 
discrepancy  by  defining  the  clean  area  to  be  r/3~U2  regardless  of  the  relationship  between  r/3~|/J 
and  R .  We  observe  since  the  clean  area  for  perfect  capture  (j3  =  1 )  is  identical  in  both  models, 
we  would  expect  our  equations  to  the  same  for  this  case. 

We  must  comment  that  both  models  make  two  simplifying  assumptions  about  the  cap¬ 
ture  phenomenon.  In  actual  practice,  one  particular  transmitter,  say  t\,  will  be  captured  by  a 
certain  receiver  if  the  ratio  of  its  received  power,  to  the  sum  of  the  received  powers  of  all  other 
signals  simultaneously  heard  by  the  receiver,  is  greater  than  the  specified  capture  ratio.  Letting 
P,  be  the  receiver  power  for  the  i’th  transmitter,  k  be  the  number  of  transmitters  the  receiver 
hears,  and  assuming  that  the  powers  are  sorted  into  decreasing  order 

k 

( P >  Pf1>  •  •  ■  >  P,k),  we  have  that  t]  is  captured  if  P,(/]T  P,  >  /3_l.  In  our  models  we 

approximate  the  sum  of  all  the  powers  of  terminals  r2  through  tk  by  the  power  of  the  next 
strongest  signal  t2.  This  assumption  however,  is  not  crucial  as  we  will  see  later,  since  the 
optimal  choice  of  system  parameters  tends  to  separate  transmitting  stations  so  that  on  the  aver¬ 
age  if  a  receiver  hears  more  than  one  terminal  transmitting,  with  high  probability  it  hears 
exactly  two  transmitters  and  therefore  the  above  sum  and  P.,  are  identical.  The  second  simpli¬ 
fying  assumption  we  make  is  that  capture  is  a  deterministic  phenomenon  such  that  if  the  ratio 
of  the  received  powers  is  greater  than  f3  then  the  signal  is  captured  with  probability  1.  In 


actuality  however,  capture  is  probabilistic  and  has  a  density  that  is  a  function  of  the  ratio  of  the 
received  powers  and  of  the  capture  parameter.  The  results  of  our  deterministic  model  can  be 
applied  to  this  more  realistic  model  however,  without  too  much  error,  by  using  a  value  of  /3  so 

k 

that  if  P,,/X  >  then  the  actual  probability  of  being  captured  is  greater  than  some 

/~2 

specified  confidence  probability  (say  0.95). 


2.3  Analysis  of  Model  1 


2.3.1  Expected  Number  of  Successful  Receptions 

We  first  calculate  the  probability  of  successful  reception  for  a  randomly  selected  termi¬ 
nal  in  the  network.  Let  us  assume  that  terminal  a  captures  the  transmission  of  its  closest 
transmitting  neighbor  t.  Conditioned  on  this,  a  will  successfully  receive  t's  packet  if  the  packet 
was  addressed  to  a  and  if  a  did  not  transmit  in  the  current  slot.  This  occurs  with  probability 

P[E,\no  interference ]  «=  — — ^  ^  e - — 

where  E,  is  the  event  of  a  successful  reception. 

We  can  see  this  by  first  defining  the  following  events 

E,  -  The  event  that  t  sends  to  a. 

Ed  -  The  event  that  t  sends  in  the  direction  of  a. 

N(i)  -  The  probability  that  there  are  /  other  terminals  besides  a  in  the  half  circle  of  radius  R 
from  t. 

In  Figure  2.2  for  example  their  are  k  terminals  in  the  half  circle  from  r  towards  F.  and 
transmitter  t  is  sending  in  the  direction  of  all  the  labeled  terminals  in  the  figure. 
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We  know  that 


P[E,\  -  P[E,\E(i]  P[Ed]  +  P[E,\E%  Plfil 

but  P[E,\Ed\  —  0  since  we  do  not  allow  packets  to  go  away  from  their  destination.  Since  i's 
destination  is  uniformly  distributed  over  the  plane,  a  randomly  chosen  terminal  a  within  a 
radius  R  of  r  is  in  the  direction  of  F.  is  equal  to  the  probability  that  a  lies  in  the  half  circle  of 
radius  R  directed  from  rto  F.  Since  a  is  equally  likely  to  be  in  either  half  we  have  P[Ed\  =  1/2. 
Calculating  the  remaining  term,  P[E,\Ed],  we  have 


Pl£,|£rf]  -  £  P[E,\Ed.Mi)]  PWi)\Ed] 

-o 


Each  of  these  terms  is  known,  for  />[/V(/')|£rf]  —  £(Af(/)]  is  Poisson  with  parameter  KttR2/2 
and  P[E,\£d,  ,V(/)]  =  l/(/  +  l),  since  if  there  are  /  other  terminals  besides  a.  making  a  total  of 
/  + 1  terminals,  t  will  select  one  of  them  with  equal  probability.  Recalling  that  the  average 
number  of  neighbors  is  S  —  \ir/?2,  we  have: 


P[E,\Ed]  -  £ 


/  +  1 


( -V  /  2) ' 


Combining  with  the  previous  calculations  we  obtain: 

P[£;3  -  |  P[E,\Ed\  - 

Knowing  that  the  packet  was  addressed  to  a  in  the  absence  of  interference  we  know  that  a  will 
successfully  receive  t’s  packet  if  a  does  not  transmit,  thus  giving  the  (1  -  p)  term  and  estab¬ 
lishing  the  above  expression. 


We  must  now  calculate  the  probability  that  there  is  no  interfering  traffic.  We  do  this  by 
first  conditioning  on  the  distance  between  a  and  rto  be  r  (r<£),  and  then  analyzing  two  cases: 

Case  1.  f3~h2r<R  (Figure  2.4a) 

In  this  case  a  will  receive  the  packet  if  there  are  no  other  transmitters  in 
(r,r/3~12),  the  clean  area.  This  area  is  equal  to  n(/3~u2r)2-nr2  »  nr2(\/f3-\) 
which  contains  no  transmitters  with  probability 

Case  2.  0-|,2r  >  R  (Figure  2.4b) 

The  area  that  must  now  be  clean  is  seen  to  be  (nR2  -  rrr2)  which  occurs  with  pro¬ 
bability 
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FIGURE  2.4a. 
Capture:  Case  1. 


FIGURE  2.4b. 
Capture:  Case  2. 


The  density  for  the  distance  between  a  randomly  selected  terminal  and  its  closest 
transmitting  neighbor,  can  be  easily  calculated.  Letting  X  be  the  random  variable  for  this  dis¬ 
tance,  and  knowing  that  busy  terminals  are  Poissonly  distributed  on  the  plane  with  parameter 
kp  we  have 

P[X  <  r]  -  1  -  P[no  busy  terminal  in  (O.r)l  -  1  -  e 


Hence  letting  fir)  be  the  density  for  X 


fir)  -  -P[X,  -  r]  -2  nkPre-Kp*rl 
dr 

Using  this  in  the  above,  after  unconditioning  we  obtain: 


PIE, 


■1  (1  ~P)(\  -e-*'2) 

—  nr 


N 


R8'n 

J  e-Kpi,r'-( i/fl-i)  i-rrkpre-^'''  dr 


+  J  2nkpre~Kpi,r''  dr 

R8'n 

Performing  the  integration  we  have: 

/>[£,]  -  °  [/80  -  e~*)  +  (1  -  fi)Npe~Sp\ 

For  future  equations  let  Y  k  >3f  1  —  e~Sp)  +  (1  —  fi)Npe~Sp . 


We  will  discuss  two  special  cases  to  see  that  they  are  intuitively  plausible. 

Case  1 .  /3  -  0  the  non-capture  case  for  which 

/>[£,]  -  (1  -/>)(!  -  e~wl2)  p  e~Wl> 

Intuitive  explanation  -  For  a  to  receive,  it  must  not  transmit  and  this  occurs  with  proba¬ 
bility  (1  -  p).  Receiver  a  also  cannot  be  isolated  from  other  nodes  in  the  network.  In 
particular  one  half  circle  of  radius  R  must  contain  at  least  one  other  terminal  and  this 
occurs  with  probability  (1— e-  V  2>.  Out  of  the  neighbors  in  this  region  only  1  can 
transmit  (occurring  with  probability  Npe~Sp),  and  on  the  average  that  transmitter  is 
surrounded  by  N  neighbors  and  thus  transmits  to  receiver  a  with  probability  1  IN. 
Combining  all  the  above  probabilities  we  obtain  the  same  expression  as  above. 


Case  2.  >3—1  the  perfect  capture  case  where 


/»[£,! 


(1  -/>)(!  -e-v/:) 

,v 


(1  -  e~w") 


Intuitive  explanation  -  Again  a  must  be  silent  which  occurs  with  probability  (1  -  p), 
and  must  not  be  isolated  from  nodes  in  one  half  circle  of  radius  R  which  occurs  with 
probability  (l-e“v/2).  There  must  be  at  least  one  transmitting  station  in  its  neighbor¬ 
hood  (occurring  with  probability  1  -e~y'>)  and  since  only  a's  nearest  neighbor  can  be 
successfully  received  by  a.  the  probability  it  transmits  a  packet  to  a  on  the  average  is 
I/.V. 
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To  calculate  the  expected  number  of  successes  in  the  network  (denoted  by  £)  we 
merely  have  to  multiply  the  previous  probability  by  the  number  of  nodes  n  in  the  network  to 
obtain: 

£-/»£[£,]-  -B-  (1  -p)(l  —  e~v/2)  Y  (2.1) 

/V 

We  can  check  this  equation  against  the  well  know  single-hop  slotted- ALOHA  results  for 
the  infinite  population  model  with  Poisson  traffic  statistics.  We  do  this  by  setting,  in  equation 
(2.1)  and  0,  and  then  performing  two  limit  operations.  By  letting  p— 0  as  n—°°  in 
such  a  way  as  to  preserve  the  product  G  -  np  to  be  a  constant,  we  obtain  Poisson  traffic 
characteristics  with  parameter  G.  We  then  obtain 

lim  -77(1  -  p)  i\-e~Nn)  N  p 

n-N— 00  N 

p— 0 

G—np— constant 

This  reaches  its  maximum  at  <7-1  giving  the  familiar  maximum  throughput  for  the  slotted- 
ALOHA  channel  [Abra77a]  of  1/e. 

It  can  be  seen  that  the  function  £  is  increasing  in  /3,  the  capture-parameter,  by  writing 
£  as  a  linear  function  of  /3: 

£  -  HU 3)  -  K  (0  (1  -  e'Nf  (1  -  Np))  +  N  p  e~Sp) 

where  K  -  (1  -  p)(l  -  e~‘wn).  This  describes  a  straight  line  with  slope 

m  -  1  -e~yp  (1  -Np).  But  this  slope  is  always  positive  since  m<  0  implies  eSp<\-Np 
which  is  seen  to  be  false  by  expanding  eNp  in  a  Maclaurin  series.  We  conclude  that  increasing 
the  receiver’s  ability  to  capture  signals  increases  the  expected  number  of  successes  in  the  net¬ 
work. 

Recall  that  Y ,  as  defined  above,  is  the  probability  that  there  is  no  interfering  traffic  to 
a’s  successful  reception.  We  see  that  K-p(  1  -  e~'Wp)  +  (1  -  / 0)Npe~Sp  can  be  viewed  as  a 
convex  combination  of  0  since  0  <  /3  <  1.  We  cannot  resist  interpreting  /3  then  as  a  probability 
(mathematicians  beware  !)  to  see  if  this  interpretation  provides  some  intuitive  insight  into  our 
problem.  If  we  interpret  /3  as  the  probability  of  capturing  signals,  then  for  the  non-capture 
case,  p— 0,  we  have  that  the  probability  that  a  has  no  interference  is  given  by: 

r3_ 0  -  Npe~Sp 

This  however,  is  the  probability  that  there  is  exactly  one  transmitter  within  a's  local  environ¬ 
ment  which  is  a  necessary  condition  for  non-interference  with  non-capture.  Likewise  when 
jS—1,  the  perfect-capture  case,  we  have: 


which  is  the  probability  that  there  exists  at  least  one  transmitter  within  a's  local  environment. 
This  of  course,  with  perfect-capture,  is  a  necessary  condition  for  successful  reception  since 
non-interference  is  certain  (a  always  captures  its  nearest  transmitter).  We  thus  have  the  intui¬ 
tively  pleasing  result  that  we  can  view  /3  as  being  the  probability  of  capturing  signals. 


We  now  seek  to  determine  the  maximum  number  of  expected  successes  in  the  net¬ 
work.  Certainly  this  number  must  be  less  than  n/2  since  every  successful  receiver  is  associated 
with  exactly  one  successful  transmitter.  It  is  easy  to  verify  this  analytically.  Numerically  calcu¬ 
lating  the  maximum  of  E  for  various  values  of  /3  demonstrates  that  the  maximum  expected 
number  of  terminals  in  the  network  that  could  engage  in  successful  communication  at  any 
given  slot  is  about  21%  for  perfect  capture  and  about  14%  for  the  non-capture  environment. 
The  exact  values  are  reported  in  Table  2.1.  We  note  here  that  the  values  of  N  and  p  that  max¬ 
imize  the  probability  of  success,  do  not  also  maximize  the  throughput  of  the  network  (see 
Table  2.2).  This  is  a  result  of  the  dependency  of  the  probability  of  successful  transmission  and 
the  maximum  transmission  range  R  as  seen  in  the  equation  N^knR2,  and  will  be  discussed 
in  greater  length  in  the  next  section.  Observe  that  the  results  of  Table  2.1  would  only  be  appli¬ 
cable  to  networks  in  which  all  packets  went  exactly  one- hop  to  reach  their  final  destination. 


2.3.2  Expected  Number  of  Successful  Transmissions 

Since  every  successful  packet  in  the  network  has  exactly  one  receiver-transmitter  pair 
associated  with  it,  the  number  of  successful  receivers  in  the  network  is  equal  to  the  number  of 
successful  transmitters.  Thus  our  calculation  in  this  section  should  yield  (2.1)  of  the  last  sec¬ 
tion.  The  approach  in  this  section  however,  allows  us  to  calculate  the  density  for  the  distance 
covered  in  a  successful  transmission,  and  serves  to  support  our  assumption  of  random  routing. 
In  analyzing  an  optimal  routing  algorithm,  the  spatial  information  used  in  determining  the  next 
terminal  in  the  route  must  also  be  used  in  calculating  the  probability  of  successfully  transmit¬ 
ting  a  packet  to  that  terminal.  Such  calculations  in  [Klei78a]  where  the  authors  assumed  an 
optimal  routing  policy,  were  performed  separately  and  produced  an  inconsistency  in  their 
results.  Since  a  determination  of  the  maximal  transmitted  distance,  say  z,  implies  in  their  cal¬ 
culation  that  there  are  no  terminals  in  area  A  of  Figure  2.3,  the  calculation  of  the  probability 
of  success  is  the  probability  there  are  no  transmitters  within  a  radius  of  R  of  the  receiver.  This 
in  Figure  2.3  is  equal  to  where  A'  is  the  intersection  of  circle  C  and  area  A.  In 

[Klei78a]  knowledge  that  A'  contains  no  transmitters  was  not  used  in  their  calculation  of  the 
probability  of  successful  reception.  With  this  inconsistency  the  expected  number  of  successful 
receivers  in  the  network  is  not  equal  to  the  expected  number  of  successful  transmitters.  The 
two  ways  we  calculate  the  expected  number  of  successes  in  the  network  do  however  agree  and 
thus  assure  us  that  our  assumption  of  random  routing  is  consistent  with  both  derivations. 

To  calculate  the  expected  number  of  successful  transmissions  in  the  network  over  one 
slot,  we  first  calculate  the  probability  that  a  randomly  chosen  station  will  successfully  transmit  a 
packet  to  a  neighboring  station,  and  then  multiply  this  probability  by  the  number  of  stations  in 
the  net,  n.  Let  us  define: 

X  -  a  random  variable  denoting  the  distance  between  the  transmitter  and  its  intended  receiver. 

E,  -  The  event  that  the  randomly  chosen  terminal  transmits  a  packet  in  the  given  slot. 

E,  -  The  event  that  there  exists  at  least  1  other  station  in  the  half  circle  of  radius  R  in  the 

direction  of  the  intended  final  destination. 

£s  -  The  event  of  successful  reception  of  the  transmission  by  its  intended  neighboring  termi¬ 
nal. 


Since  we  assume  that  routing  is  uniform  over  all  the  nodes  in  the  right-hand  region  of 
Figure  2.5  we  have: 


Differentiating  we  obtain: 


w  <  -  fim -  f 


P[r  <  AT  +dr\E,,E,)  -  tfr 


Our  definition  of  the  capture  phenomena  for  Model  1  implies: 

P[r  <  X^r  +  dr,Es\E„E,]  —  (1  -  p)  dr  (2.2) 

A 

The  probability  of  successful  transmission,  P[ES),  is  found  after  unconditioning  and  integrating: 


/>(£,]-  p  (1  -  p)  (1— <TV/2) 


f  2  r 

\  R2  6 


minir0-'n.R)\ 


dr 


Breaking  this  integral  up  into  regions  [0,/?/3l/Jl  and  (RfiU2,R}  we  obtain: 


(2.3) 
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Thus  the  expected  number  of  successful  transmissions  is  found  to  be 


£  -  (-£)  (1  -  P)  (l-e-v/2)  Y 
N 

which  is  the  same  as  (2.1). 


2.3.3  Expected  Forward  Progress 


We  are  now  in  a  position  to  derive  the  density  for  the  distance  between  transmitter  and 
receiver  for  a  successful  transmission.  Referring  back  to  equation  (2.2)  and  using  (2.3)  we 
have. 


P[r  <  X  <  r+dr\Es] 


Plr  <  X  <  r+dr ,  FJ 
P[ES) 


p  (1  —  p)(l  —  e-v/J)  2  r  —  >| 
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or  defining  g(r)  io  be  the  density  for  r  and  —  •  we  may  rewrite  equation  (2.4)  in  regions 

S  I 

to  obtain 


g(r) 


-  ^4  e~Kpwrl'B  0  <  /•  £/31/2 

s  /r 

1  ir  ^  r  <  R 

s  /r 


(2.5) 


It  can  be  easily  verified  that  this  integrates  to  1  and  thus  is  a  proper  density. 


Suppose  as  shown  in  Figure  2.6,  transmitter  t  is  sending  a  packet  to  final  destination  F 
through  intermediate  node  a.  We  wish  to  calculate  the  progress  of  the  packet  towards  its  final 
destination.  To  simplify  the  calculation  we  assume  that  forward  progress  will  be  the  same  for 
any  node  on  the  line  perpendicular  to  the  direction  of  the  destination,  line  L  in  the  figure.  This 
assumption  is  reasonable  if  the  distance  D  is  much  greater  than  R.  Because  terminals  are  ran¬ 
domly  distributed  on  the  plane,  for  a  given  distance  r  and  a  given  destination  F  the  angle  0 
will  be  uniformly  distributed  over  (-n/2,nl2).  Define  Z  to  be  the  random  variable  denoting 
the  forward  distance.  We  see  that  for  a  given  r  the  probability  that  Z  is  less  than  some  value  r 
is  the  same  as  the  probability  of  |0|  being  larger  than  cos-l(r/7)  or  letting  F< r)  be  the 


FIGURE  2.6. 

Expected  Forward  Progress. 


distribution  of  Z  we  have 

F(z\r) 


1  r  <  z 

2  cos~l(z/r) 
v 


1  - 


0  ^  2  <  r  =$  R 


Differentiating  this  with  respect  to  z  we  derive  the  conditional  density 

0 


f(z\r) 


n  r  Vl-z2/  r2 
We  can  now  calculate  the  expected  progress  given  r 

£(Z|r]  -  fz  f(z\r )  dz  -  —  f 


r  <  z 

0  <  z  <  r  ^  R 


■dz  — 


We  can  uncondition  this  by  using  the  density  of  (2.5)  to  obtain 


ElZ]  g(r) 

*  o 


S  R2  IT 


Raj'2 

0 


v*  r 

J  2  r 2  e-*'"2'*  dr  +  e~Kpi,R 2  j  2  rdr 

I '2 


(2.6) 


The  second  integral  within  the  brackets  is  equal  to  y  R}  (1  -  0V0)  e  KpvR  and  letting  the 
first  integral  be  denoted  by  INT ,  we  can  write  equation  (2.6)  as 

7-  +  JJL  <,-*">  (1  -fijji)  (2.7) 

s  R1  n  3  s  ir 

In  Appendix  A  it  is  shown  that 

Jo  2  t  ft  (2/  +  l>! 

and  thus  using  this  for  INT  'm  (2.7)  and  reducing  we  finally  obtain: 

=  -7Te""R[^|,i^L  +  T(,-H  ,2-8) 

For  future  equations,  let  Q  k  £&±<m£  + 1  h  . ejfr. 

fyp  /w|  \2j t 1J!  J 


2.3.4  Expected  Throughput 

We  can  calculate  the  expected  throughput  for  the  network  for  each  slot.  For  any  ran¬ 
domly  selected  terminal,  the  expected  path  length  between  it  and  another  randomly  selected 

1/2 

1 28  n 

terminal  is  given  in  [Kend63a]  as  d  —  — - .  Since  :,  as  calculated  in  the  previous  sec- 

4577  \ir 

tion  is  known,  the  number  of  hops  h  a  randomly  selected  packet  will  take  is  given  by  h  -  d/z. 
Therefore  the  average  number  of  messages  delivered  to  their  final  destinations  per  slot,  the 
throughput,  is  given  by  y  -  nP[E,]/h  —  nP[E,)  zld.  Using  the  previously  derived  equations 
for  the  quantities  above  we  obtain  the  final  result: 


yifrN.pJ  -  (1  -  p)i\  -  e~s'2)  p  e~Sp  Q 

The  increase  in  the  throughput  7  with  the  square  root  of  the  number  of  terminals  in  the  net¬ 
work  is  a  result  of  the  spatial  reuse  of  the  channel.  Observe  that  the  equation  obeys  our  intui¬ 
tion  for  p  —  0  or  p  -  1  when  the  throughput  is  zero,  and  that  the  U-?_v'2)  term  is  the  proba¬ 
bility  that  the  network  is  connected  over  one  hop.  Once  again  we  can  show  that  the  throughput 
of  the  system  is  an  increasing  function  of  0  since  if  the  function  is  increasing  in  a  -  0v'0  then 


it  is  also  increasing  in  0.  Thus  7  is  linear  in  a  with  slope  m  -  Y 


4  (.V/?)'-'  /'! 


-  -7  which 


is  clearly  positive  since  the  minimum  of  m  is  2/3.  Thus  as  we  have  seen  previously,  increasing 
the  capture-parameter  can  only  increase  performance. 


2.3.S  Discussion  of  Results 

In  all  the  following  graphs  and  Tables  we  use  normalized  throughput 
y'(/V,p, /8)  «  y (N,p,p)/yfn ,  hence  eliminating  the  dependency  of  the  size  of  the  network  from 
our  equations.  We  must  note  that  the  square  root  of  n  dependency  on  the  throughput  of  the 
network  is  an  important  factor  that  lets  us  achieve,  by  reducing  the  strength  of  transmitted  sig¬ 
nals,  throughputs  greater  than  that  obtainable  by  running  the  network  as  a  one-hop  ALOHA 
network.  For  example,  in  Table  2.2  for  f3  -  0.7,  we  see  that  y'  -  0.0749282.  To  determine  the 
number  of  terminals  needed  in  the  network  to  achieve  throughput  greater  than  Me  we  set 
yfn  y'  >  Me  which  implies  n  =  24.  Thus  in  a  network  with  more  than  24  terminals,  it  pays  to 
voluntarily  limit  the  transmission  ranges  so  that  on  the  average  only  /V*4.99725  other  termi¬ 
nals  are  within  hearing  range.  Figure  2.7  demonstrates  graphically  the  result  we  saw  in  the  pre¬ 
vious  pages  that  increasing  the  capture-parameter  improves  system  performance.  Here  we  plot 
y  as  a  function  of  p  for  a  fixed  N  and  various  values  of  13.  In  Table  2.2  we  have  listed  the 
maximum  y'  over  all  possible  /V  and  p  values  for  a  fixed  j8  value.  We  note  again  that  the  y'  is 
increasing  in  /3.  Observe  that  the  spread  of  the  optimal  values  of  N  and  p  over  all  /3  values  is 
not  wide.  Since  these  values  do  not  change  substantially  with  0,  the  capture-parameter  of 
packet  radios  in  the  network  does  not  need  to  be  known  to  a  high  degree  of  accuracy  to  deter¬ 
mine  the  network’s  W  and  p  values  that  achieve  optimal  performance.  In  Table  2.2  we  have 
also  listed  the  probability  of  a  successful  transmission  as  well  as  the  expected  forward  progress 
for  the  same  N  and  p  values.  By  comparing  Tables  2.1  and  2.2,  we  observe  that  values  of  .V 
and  p  which  maximize  y'  do  not  also  maximize  /*[£’,].  Although  it  might  seem  intuitive  that 
maximizing  the  number  of  successes  in  the  network  by  picking  an  optimal  transmission  range 
R  and  hence  by  picking  <V  -  KnR2  would  increase  the  throughput  of  the  system,  a  little 
thought  shows  that  this  is  not  necessarily  true.  We  can  see  this  from  Table  2.1  where  the  .V 
values  that  maximize  P[E,]  are  seen  to  be  small,  approximately  2.6  for  13  -  .7.  The  network  is 
in  this  case  divided  into  many  receiver-transmitter  pairs  in  an  attempt  to  take  full  advantage  of 
the  spatial  reuse  of  the  channel  and  although  this  increases  the  probability  of  successful 
transmission,  packets  in  such  an  environment  must  pass  over  many  hops  before  reaching  their 
final  destinations.  This  tends  to  decrease  the  number  of  packets  reaching  their  final  destina¬ 
tions  in  any  one  slot,  and  thus  reduces  the  throughput  of  the  system.  This  tradeoff  between 
the  probability  of  success  and  the  throughput  of  the  system  as  governed  by  the  number  of  hops 
between  source  and  destination  is  a  fundamental  issue  of  multi-hop  systems  and  occurs  in 
several  guises.  For  example,  we  have  already  shown  that  increasing  (3  will  increase  y'.  This 
increase  in  y  can  result  from  an  increased  /Mf,!,  an  increased  r  (thus  increasing  the  average 


0  N 


4.33261 
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4.41670 
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4.58656 
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4.83818 

4.99725 

5.17892 

5.38063 

5.59807 


.18012 

0.18193 

0.18531 

0.18978 

0.19523 

0.20155 

0.20867 

0.21647 

0.22474 

0.23324 

0.24164 


0.0584586 
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0.0605900 

0.0624732 

0.0648259 

0.0676626 

0.0710160 

0.0749282 

0.0794432 

0.0845999 

0.0904239 


0.05991 

0.06285 

0.06577 

0.06879 

0.07190 

0.07526 

0.07875 

0.08242 

0.08624 

0.09022 

0.09433 


0.42441 

0.40834 

0.39687 

0.38797 

0.38094 

0.37541 

0.37120 

0.36823 

0.36649 

0.36601 

0.36682 


TABLE  2.2 

number  of  hops  a  packet  takes  from  source  to  destination),  or  a  combination  of  both.  We  see 
in  Table  2.2  that  as  (3  increases  from  0  to  .9,  PffJ  increases  and  z  decreases.  Thus  for  optimal 
throughput  packets  must  travel  over  more  hops  but  they  “hop”  more  frequently,  once  again 
showing  the  tradeoff  between  P[ES\  and  z. 


In  Figure  2.8  we  show  the  relationship  of  y'  as  a  function  of  p  for  fixed  /V  and  £.  We 
notice  that  for  any  N,  optimal  performance  is  degraded  for  small  changes  of  p  from  its  optimal 
value,  p\  but  that  as  N  increases,  the  curves  around  this  p‘  become  narrow.  This  variation  of 
y  for  large  N  results  from  the  fact  that  the  transmission  of  any  packet  radio  interferes  with  a 
larger  number  of  other  terminals.  This  increase  in  the  number  of  collisions  increases  the  sensi¬ 
tivity  of  the  throughput  for  perturbations  of  the  transmission  probability  from  its  optimal  value. 
We  also  observe,  for  distinct  values  N !  and  ;V2,  such  that  /V,  >  jV2,  there  exists  a  transmission 
probability,  p"  such  that  for  all  P\  <  p‘ <  Pz  then  y'(Nz,Pu(3)  >  y'(jVi,pt,/3)  and 
y'(Ni,P2,t 3)  <  y'(N\,pz,fi).  This  can  best  be  seen  by  looking  at  Figure  2.8  where  we  see  that 
for  the  two  model  1  curves  for  W-2  and  iV—5,  there  exist  p*~.45  such  that,  for  all  p>p‘  the 
/V-2  curve  dominates  (has  a  greater  throughput)  the  .V—5  curve,  and  the  opposite  holds  for 
p>p‘ .  Intuitively,  we  expect  that  for  large  p,  systems  with  small  N  would  have  better  perfor¬ 
mance  characteristics  since  being  less  densely  populated  they  are  inherently  more  immune  to 
interference  than  systems  with  higher  N  values.  On  the  other  hand  for  small  p,  since  the  time 
between  transmissions  grows  larger  in  the  small  /V  case,  one  would  expect  that  systems  with 
higher  /V  would  have  better  performance. 


FIGURE  2.8 

Normalized  Throughput  as  a  Function  of  p. 
for  f3  -  0.7  and  ;V  -  2. 5.  and  8. 


We  can  unify  our  discussion  of  Figure  2.8  by  defining  the  offered  load  per  unit  area  to 
be  G  «  Sp.  From  previous  results  [Abra73a]  for  finite  population  slotted-ALOHA  networks, 
we  know  that  G  -  1  optimizes  network  throughput.  In  the  multi-hop  environment  however, 
connectivity  of  the  network  must  be  preserved.  This  consideration,  as  noted  before,  manifests 
itself  in  the  c  -  (1— e~y/2)  term  appearing  in  the  equation  for  y.  If  p  is  large,  G  -  1  implies 
that  N  —  \/p  =  1  which  tends  to  disconnect  the  network  since  c  =  .39.  Obviously  optimal 
throughput  for  this  case  would  have  S  >  1  and  thus  G  >  1.  We  would  thus  expect  G  —  1 
only  in  cases  where  p  is  small  enough  to  make  S  sufficiently  large  to  assure  connectivity.  We 
must  however  take  account  of  capture  in  discussing  the  offered  load.  In  the  non-capture 
environment  we  would  expect  the  offered  load  that  maximizes  throughput  to  be  less  than  that 
for  the  capture  environment  because  the  probability  a  transmission  suffers  a  collision  is  greater 
for  0-0  than  for  (3-1.  Thus  increasing  G  -  Np  has  a  greater  effect  in  increasing  the 
number  of  expected  collisions  in  environments  with  non-capture  than  for  those  with  perfect 
capture.  To  check  this  intuition,  we  numerically  calculated  the  N  value  that  maximized  y  for 
fixed  0  and  p  and  plotted  G  -  Np  against  p  in  Figure  2.9.  We  see  that  curves  for  high  0 
values  dominate  those  for  lesser  values,  justifying  our  belief  that  the  offered  load  that  maxim¬ 
izes  throughput  can  be  greater  for  larger  0,  and  that  as  p  increases  so  does  G.  illustrating  the 
relationship  of  the  connectivity  factor  c  has  in  sparse  environments.  We  can  lend  some 
mathematical  insight  into  these  graphs  by  defining  the  effective  number  of  neighbors,  S',  to 
satisfy  S  —  S'/ (  1  -  e~s  n) .  One  then  interprets  N  as  the  the  average  number  of  terminals  per 
unit  area  that  results  when  we  randomly  distribute  terminals  on  the  plane  with  an  average  den¬ 
sity  of  S'  terminals  per  unit  area,  and  condition  upon  having  a  connected  network.  The 
(1  -  e~sir)  term  is  thus  seen  to  be  the  conditional  probability  of  hop  connectivity.  Using  S' 
instead  of  N  in  calculating  the  offered  load,  G'  -  p  S'  =  1  implies  that  S'  =  \/p  and  using 
this  in  the  above  equation  that  defined  S',  we  would  have  S  =  S  =  \/p(l  -  e~U2p).  To  check 
this  intuition,  we  used  the  values  obtained  in  generating  the  offered  load  curves  of  Figure  2.9 
for  0  -  1  where  we  have  seen  that  for  low  p  values.  G  ~  1.  In  Table  2.3  we  produce  the  S 
and  p  values  that  optimize  y',  as  well  as  their  product,  the  offered  load  G  in  the  left  part  of  the 
Table.  In  the  right  we  tabulate  the  hypothesized  values  using  the  effective  number  of  neigh¬ 
bors  S'.  We  see  that  S  and  S  are  approximately  equal  and  that  the  effective  number  of  neigh¬ 
bors,  S'  is  strictly  less  that  S.  Comparing  G  and  G'  we  see  that  G'  is  much  closer  to  1 
throughout  the  range  of  p  lending  support  to  our  previous  intuitive  arguments. 

In  our  last  plot  for  this  section.  Figure  2.10,  we  graph  y  as  a  function  of  S  for  0-7 
and  various  values  of  p.  Observe  that  for  the  near  optimal  S  for  p  =  .2.  (namely  S  -  5)  the 
curve  is  very  fiat.  This  implies  it  is  not  necessary  to  determine  S  to  a  high  degree  of  accuracy 
to  achieve  near  optimal  performance. 
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FIGURE  2.9 

G  -  Sp  as  a  Function  of  p 
for  optimal  N  and  p,  and  0  -  0.7  and  1.0 
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Recall  in  Model  2  we  assume  that  any  other  transmitter  within  r/3— 1/2  of  a  packet  radio 
receiving  a  packet  from  another  transmitter  a  distance  /•  away,  will  cause  a  collision  on  the 
channel.  We  thus  do  not  need  to  divide  the  clean  area  into  two  regions.  We  will  be  brief  in 
describing  the  results  of  this  section,  since  most  derivations  follow  lines  similar  to  those  in 
model  1. 


2.4.1  Expected  Number  of  Successful  Receptions 


The  clean  area  for  a  transmitter  at  a  distance  of  r  from  the  receiver  is  now  r) 3  1  2.  and 
thus  the  probability  of  this  area  has  no  other  transmitters  is  e~Kpirr'l&.  we  thus  have: 


p[Es)  -  p  (1  —  />)(!  -  <TV/2)  j  ~y  e-Kpw'2,a  dr 


Integrating  this  over  ft),/?]  yields: 


/»(£,]  -  m ^ (l-e-W*) 

/V 


and  hence  the  number  of  successful  receptions  is 


E  -  jj  /?  (1  -  p)(  1  -  e~v/2)  (l-e-Vp/a) 


Once  again,  the  expected  number  of  successes  in  the  network  is  an  increasing  function  in  13. 
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FIGURE  2.10 

Normalized  Throughput  as  a  Function  of  /V 
for  0  —  0.7  and  p  —  0.01.0.2,  and  0.9. 


This  can  best  be  seen  by  writing  £  as  a  function  of  0  and  showing  that  the  slope  is  positive. 
Let  H(0)  -  K  /3(l-e_v,,/fl)  where  K  -  w(l  -/>)(1  -  e_V  2)/,V.  Differentiating  H((3)  with 
respect  to  0  yields  H'(0)  -  1  -  e~ Sl>ni(Npl  0  +  1)  which  must  be  greater  than  zero  since 
gXpis  >  l  +  Npl  0.  It  makes  no  sense  to  try  to  obtain  the  limit  of  the  above  expression  for 
the  infinite  population,  Poisson  traffic  slotted- ALOHA  model  as  we  did  in  model  1.  The  reason 
for  this  concerns  our  definition  of  capture  for  model  2.  If  we  let  0  go  to  zero,  the  capture 
radius  r  0~’h  goes  to  infinity.  Since  all  our  derivations  assume  terminals  to  be  Poissonly  distri¬ 
buted  over  the  plane  with  parameter  X,  this  infinite  capture  radius  will  contain  an  infinity  of  ter¬ 
minals  which  for  any  p  >  0  will  have  an  infinite  number  of  transmitters  with  probability  1, 
thus  guaranteeing  a  certain  collision.  The  equation  above  rightly  indicates  that  for  0  —  0  the 
expected  number  of  successes  is  zero.  In  model  1,  since  the  capture  radius  was  limited  to  a 
maximum  of  R,  it  did  make  sense  to  let  0— 0  since  the  capture  radius  was  bounded.  Once 
again  we  can  determine  analytically  that  £  <  n/ 2  and  produce  Table  2.4  which  contains  the 
probability  of  success  for  various  0  values.  Observe  that  for  /3  =  1  these  results  agree  with 
those  of  model  1 . 
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0.20379 

0.02737 

0.2 

1.5243 

0.24990 

0.04467 
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1.8104 

0.27811 

0.05794 

0.4 

2.0412 

0.29824 

0.%775 

0.5 

2.2373 

0.31373 

0.U7789 

0.6 

2.4089 

0.32622 

0.08578 

0.7 

2.5621 

0.33660 

0.09272 

0.8 

2.7011 

0.34543 

0.09889 

0.9 

2.8284 

0.35307 

0.10443 

1.0 

2.9462 

0.35977 

0.10946 

TABLE  2.4 


2.4.2  Expected  Forward  Progress 

We  first  calculate  the  density  for  the  distance  between  a  successful  transmitter  and  its 
receiver.  Using  the  same  definitions  as  in  the  previous  section,  we  have 


we  can  write  g(r)  as 


g(r)  ^  —  e-Krw*a 

s  R 1 

We  thus  calculate  the  expected  forward  progress  for  each  transmission  to  be 


,  K 

-4  f 

s  n  R2ja 


r 2  e~krl  dr 


where  k  —  Xnp/p.  Using  the  result  from  Appendix  A  to  expand  the  integral  and  then  simpli¬ 
fying,  we  finally  obtain 


2.4.3  Expected  Throughput 


Using  the  previous  results,  we  can  now  write 


y  —  — 'Jn/N  0  (1  —  p)(l  —  e  ,v/2)  e~Np/B  Y.(^Nplp)1 


(2j+\)\ 


(2.10) 


Again  this  is  an  increasing  function  of  0,  has  the  square  root  dependency  on  the  number  of 
nodes  in  the  network,  and  explicitly  accounts  for  p  *  0,  p  -  1  cases  and  network  connectivity. 


2.4.4  Discussion  of  Results 


To  compare  the  performance  of  the  two  models  we  produced  the  same  Tables  and 
graphs  for  Model  2  as  we  did  for  Model  1 .  Earlier  we  observed  that  the  two  models  should  be 
identical  for  0  —  1  since  the  clean  area  in  both  models  for  this  case  are  identical.  Algebraic 
comparison  between  similar  formulas  from  both  models  shows  that  they  are  equal  for  0  -  1 
and  the  results  in  this  section  show  that  as  0  grows  larger,  the  results  of  the  two  models  are 
approximately  equal.  We  therefore  restrict  our  discussion  of  the  curves  and  Tables  in  this  sec¬ 
tion  to  relevant  differences  between  the  two  models. 
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TABLE  2.5 


We  observe  in  Table  2.5  that  for  model  2  there  is  a  wider  spread  for  optimum  N  and  p 
values  for  the  range  of  fi  values.  This  is  due  to  the  larger  capture  radius  for  small  fi.  One 
curious  difference  between  this  and  Table  2.2  is  that  here  as  fi  increases  so  does  r,  whereas  in 
model  1  we  observed  a  decrease  in  the  z  values.  The  increase  for  model  2  can  be  explained  by 
the  fact  that  for  low  values  of  fi,  small  values  of  r.  the  distance  between  transmitter  and 
receiver,  have  higher  probability  than  greater  r  values  because  the  capture  radius  for  these 
values  is  large.  As  /3  increases,  the  capture  radius  for  a  fixed  r,  decreases  in  model  2  and  thus 
larger  r  values  are  more  heavily  weighted  and  thus  increase  z.  This  explains  the  increase  in  : 
for  model  2.  In  model  1  however,  the  capture  radius  was  bounded  to  be  less  than  R .  Thus  the 
clean  area  is  small  for  terminals  that  lie  close  to  R  and  this  increases  the  probability  that  the 
projected  distance  will  be  large.  This  is  the  reason  why  the  throughput  for  model  1  is  always 
larger  than  that  for  model  2. 


In  Figure  2.7  we  plot  y  as  a  function  of  p  for  <Y  -  5  and  various  fi  values.  Once  again 
we  see  the  dominance  of  the  higher  /3  curves.  We  can  also  observe  that  curves  for  model  1 
dominate  those  for  model  2.  Comparison  of  the  curves  for  the  two  models  on  this  plot  shows 
the  dominance  of  throughput  values  of  model  1  over  those  of  model  2.  This  dominance  is  a 
result  that  the  probability  of  collision  is  much  larger  in  the  second  model  for  any  value  of  V 
and  p  and  hence  decreases  the  expected  throughput  of  the  system.  Figure  2.9  shows  again  the 
fluctuations  of  the  offered  load  for  increasing  p  values  and  we  note  that  the  curves  for  (3  -  0.7 
is  almost  identical  to  the  same  curve  for  model  1.  Figure  2.10  shows  the  flatness  of  throughput 
curves  for  fixed  fi  and  p ,  over  values  of  N. 


2.5  Conclusions 


In  this  chapter  we  have  analyzed  two  models  of  capture  in  a  random  planar  network 
where  slotted-ALOHA  was  used  to  broadcast  packets  on  the  channel.  The  results  of  the  two 
models  are  similar  for  capture-ratios  achievable  on  good  FM  receivers  and  thus  either  could  be 
used  to  analyze  networks  of  this  kind.  Since  increasing  the  capture  parameter  increases  the 
throughput  of  the  network  we  can  conclude  that  capture  is  a  beneficial  phenomenon  and 
throughout  the  following  chapters  all  of  our  network  models  will  assume  that  radios  can  capture 
signals.  The  tradeoff  between  the  probability  of  a  successful  transmission  and  the  expected 
number  of  hops  taken  by  a  packet  in  the  network  has  been  delineated,  and  we  have  seen  that 
even  in  ideal  conditions  with  perfect  capture  and  one-hop  messages,  no  more  than  21%  of  the 
nodes  in  the  network  on  the  average,  are  engaged  in  productive  communications  in  each  slot. 
The  square  root  of  n  dependency  on  the  throughput  has  been  shown  to  substantially  increase 
the  throughput  of  the  network  over  conventional  one-hop  ALOHA  networks  when  the  number 
of  nodes  in  the  network  is  sufficiently  large.  The  critical  parameter  to  network  optimization  has 
been  shown  to  be  the  value  of  p ,  the  probability  of  transmitting  in  any  given  slot  and  we  have 
seen  that  an  offered  load  of  (7=1  optimizes  network  performance. 

In  the  following  chapter  we  eliminate  the  slotted-ALOHA  assumption  from  this  model 
and  address  a  more  basic  question.  In  that  chapter  we  seek  to  determine  the  throughput  of  the 
best  protocol  in  the  mobile  multi-hop  environment.  We  will  derive  a  tight  upper  bound  for  the 
performance  of  an  optimal  protocol  in  this  environment  and  then  use  that  as  a  standard  to 
determine  how  well  slotted-ALOHA  performs  in  comparison  to  this  approximation.  Some  of 
the  observations  we  have  made  in  this  chapter,  we  will  see,  generalize  to  this  optimal  protocol. 


CHAPTER  3 

AN  UPPER  BOUND  FOR  THROUGHPUT  FOR  ALL  PROTOCOLS  IN  A  MOBILE  NETWORK 

From  Chapter  2  we  know  how  well  slotted-ALOHA  performs  in  a  multi-hop  mobile 
packet  radio  network.  In  this  chapter  we  remove  the  slotted-ALOHA  assumption  from  the  net¬ 
work  model  and  determine  an  upper  bound  on  the  best  performance  for  any  protocol  in  a 
mobile  multi-hop  network. 


3.1  Introduction 

In  the  previous  chapter  we  analyzed  a  packet  radio  network  consisting  of  randomly  dis¬ 
tributed  nodes  on  the  plane  and  determined  optimal  range  and  transmission  probabilities  for  the 
slotted-ALOHA  protocol.  One  of  the  parameters  of  the  system,  N ,  the  average  number  of 
neighboring  nodes,  contained  the  range  factor  (since  .V-Xn-/?2.  where  A  was  the  density  of 
nodes  on  the  plane  and  R  was  the  fixed  range  of  transmission)  and  was  varied  to  maximize  the 
throughput  of  the  network.  In  this  chapter  we  assume  we  are  given  a  connected,  random,  net¬ 
work  with  a  fixed  N  (it  is  not  a  parameter  that  can  vary)  and  seek  to  determine,  under  these 
assumptions,  the  maximum  number  of  simultaneous  transmissions  in  the  network  that  can  be 
successful.  This  number,  divided  by  the  number  of  nodes  in  the  network,  will  give  the  max¬ 
imum  fraction  of  successful  transmissions  that  are  possible  in  networks  of  this  type  and  will  be 
an  upper  bound  for  any  random  access  protocol.  This  can  be  used  then  as  a  standard  to  evaluate 
the  performance  of  other  protocols.  In  Chapter  4  we  will  refine  the  definition  of  an  optimal 
protocol  of  this  type  and  will  derive  its  delay  characteristics. 


3.2  Model  and  Analysis 


We  will  assume  that  nodes  of  a  packet  radio  network  are  distributed  on  the  plane 
according  to  a  Poisson  point  process  with  a  mean  density  of  A  radio  units  per  unit  area.  With 
no  loss  of  generality  we  will  assume  that  radios  transmit  with  a  range  of  one  unit  (R-l).  The 
average  number  of  neighbors  they  have,  N.  therefore,  is  given  by  .V-a n.  From  any  given 
node,  a,  another  node,  b ,  is  said  to  be  /-hops  away  if  there  exists  a  path  from  a  to  b  that  con¬ 
tains  /—I  other  nodes  and  no  other  path  exists  between  a  and  b  that  contains  fewer  nodes.  In 
our  derivation  we  will  focus  on  a  section  of  the  network  containing  n  nodes,  and  thus  will  be 


concerned  with  a  portion  of  the  network  of  radius  R o.  where  n-\ir/?o  Our  final  results  are 


independent  of  n.  We  will  assume  the  network  formed  by  these  n  nodes  is  connected.  We 
should  mention  here  that  our  network  model,  as  in  Chapter  2,  is  that  of  a  mobile  packet  radio 
network  (thus  the  random  distribution  of  nodes  on  the  plane).  In  such  networks  it  is 
extremely  difficult  for  nodes  to  ascertain  the  locations  of  other  nodes  in  the  network,  and,  in 
particular,  nodes  do  not  know  the  direction  of  the  recipients  of  their  transmitted  packets.  Our 
results  here  are  applicable  only  to  protocols  that  do  not  utilize  information  concerning  the  loca¬ 
tion  or  direction  of  their  neighboring  nodes. 


To  motivate  how  we  propose  to  calculate  our  upper  bound  suppose  a  node,  a,  is 
transmitting  a  packet  to  one  of  its  neighbors.  Let  S,  be  the  set  of  nodes  that  are  /-hops  away 
from  a,  and  define  a  k -order  independent  set  to  be  a  set  of  nodes  that  are  all  mutually  k  or  more 
hops  away  from  each  other.  A  maximal  k -order  independent  set  is  a  /c -order  independent  set  to 
which  no  other  node  of  the  network  can  be  added.  It  is  clear  that  if  all  nodes  of  a  maximal  3- 
order  independent  set  transmit,  and  these  are  the  only  transmitting  nodes  in  the  network,  all  of 
their  transmissions  will  be  successfully  received,  where  by  this  we  mean  that  nodes  in  the  net¬ 
work  that  are  possible  recipients  of  any  of  these  messages,  hear  exactly  one  transmitter. 


FIGURE  3.1 

A  maximal  3 -order  independent  set. 


An  example  of  such  a  set  is  shown  in  figure  3.1.  In  this  figure,  nodes  of  the  3-order  indepen¬ 
dent  set  are  indicated  by  the  larger  circles.  We  can  easily  show  the  following  property  of  maxi¬ 
mal  3-order  independent  sets: 


Lemma:  For  any  maximal  3-order  independent  set  M ,  no  node  not  in  the  set  can  transmit 
without  causing  interference  with  a  transmission  of  at  least  one  node  in  M. 


Proof:  Suppose  a  was  such  a  node.  Using  the  previous  definitions,  it  is  clear  that  S]  P|  Af-0 
(where  0  is  the  null  set)  otherwise  a  would  be  interfering  with  its  own  reception  from  nodes  in 
M.  If  S2  Pi  M  0  then  this  implies  there  is  a  node  that  simultaneously  receives  signals  from 
a  and  also  from  another  node  in  M.  Thus  S2PI  A/  — 0,  but  this  implies  that  M  is  not  maxi¬ 
mal  since  a  is  at  least  3  or  more  hops  away  from  every  node  in  M. 


We  should  observe  that  this  lemma  is  also  true  for  maximal  A: -order  independent  sets 
where  fr-4,5,  and  thus  we  cannot  immediately  conclude  that  a  maximal  3-order  independent 
set  corresponds  to  the  greatest  number  of  transmissions  in  the  network  that  are  guaranteed  to 
cause  no  collisions.  Intuitively,  however,  to  achieve  maximal  throughput  we  would  want 
transmitters  to  be  as  close  to  each  other  as  possible,  without  having  collisions  detract  from  the 
throughput  of  the  channel.  To  make  this  more  precise,  let  Sk  be  the  set  of  all  maximal  k -order 
independent  sets,  and  let  Lk  be  the  cardinality  of  the  largest  set  in  Sk  (we  are  assuming  here  a 
finite  but  arbitrarily  large  graph).  Then,  since  a  k+ 1-order  independent  set  is  also  a  k- order 
independent  set,  we  have  that  Lk  ^  Lk+\.  Thus,  since  collisions  occur  for  k  —  1 , 2.  we  can  con¬ 
clude  that  the  largest  number  of  successful  transmissions,  without  allowing  collisions,  is  given 
by  £3. 


Since  nodes  in  this  scenario  that  are  allowed  to  transmit  are  mutually  at  least  three  hops 
away  from  each  other  and  at  best  exactly  three  hops,  in  the  ideal  case  we  can  imagine  the  plane 
being  tessellated  with  equilateral  triangles  having  sides  equal  to  the  average  distance  between 
nodes  three  hops  away.  This  tessellation  is  motivated  in  figure  3.2  where  we  have  connected 
the  nodes  of  the  maximal  3-order  set  of  figure  3.1  (not  all  such  configurations  will  result  in  a 
hexagonal  shaped  figure).  For  such  a  tessellation  each  vertex  corresponds  to  a  transmitting 
node  of  the  network.  The  number  of  such  triangles,  for  a  given  section  of  the  network,  will 
correspond  to  twice  the  number  of  vertices.  This  can  be  seen  in  figure  3.3  where  we  have 
mapped  each  vertex  to  the  triangle  lying  directly  above  it  (shown  by  the  arrows  in  the  figure). 
The  shaded  triangles  have  no  corresponding  vertices  and  can  then  be  seen  to  be  equal  in 
number  to  those  that  are  mapped  to  vertices.  If  we  let  X  be  the  average  distance  between 
nodes  three  hops  away,  then  the  area  of  each  triangle  is  given  by  v3 X'!A.  Letting  T  be  the 
number  of  triangles  in  network  under  consideration  whose  total  area  equals  nRo ,  we  can  write: 

T  -  /(V3.?2/4) 
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/(/V)  -  .9068/ ,V 


The  data  that  lies  above  the  /(.V)  curve  is  a  result  of  the  edge-effects  obtained  from  generating 
finite  graphs.  In  such  graphs,  nodes  at  the  edge  of  the  graph  are  more  likely  to  belong  to  maxi¬ 
mal  independent  sets  since  they  have  neighbors  only  on  one  side.  This  tends  to  increase  the 


size  of  the  maximal  3-order  independent  sets. 

It  is  interesting  to  graphically  see  what  X  =  2  means  in  terms  of  transmission  areas. 


FIGURE  3.5 

X  in  terms  of  transmission  areas. 

We  see  in  figure  3.5  that  this  average  third  hop  distance  implies  that  the  plane  is  covered  with 
unit  circles  from  transmitters  located  at  the  vertices  of  the  equilateral  triangles.  Observe  that 
there  is  no  overlap  of  the  circles  and  that  very  little  of  the  plane  is  not  covered.  In  fact  it  is 
well  known  that  this  arrangement  of  circles  maximizes  the  density  of  the  area  of  the  plane 
covered  by  exactly  one  circle  provided  that  no  area  is  covered  twice  [Will79a].  Viewing  the 
problem  in  this  manner  leads  us  to  inquire  about  the  arrangement,  allowing  possible  overlap  of 
circles,  that  yields  the  highest  density  of  singly  covered  area.  This  density  will  yield  an  upper 
bound  to  the  fraction  of  the  expected  number  of  successful  transmissions  for  a  mobile  packet 
radio  network  since,  in  the  absence  of  capture,  the  probability  that  a  transmission  is  successfully 
received  is  equal  to  the  fraction  of  area  that  is  covered  by  only  its  transmission.  We  can  use 
figure  3.6  to  this  end.  This  figure  is  a  reproduction  of  one  of  the  triangles  of  figure  3.5.  where 
we  have  drawn  the  triangles  with  a  length  less  than  2  units.  In  this  figure  let  I(.x)  be  the  area 
of  overlap  within  the  triangle,  and  let  E(x)  be  the  area  which  is  not  covered  in  the  circle  tessel¬ 
lation.  First  let  us  derive  equations  for  these  two  functions.  We  can  use  the  equation  derived 
in  Appendix  B  for  A  (•)  (equation  (B.7>)  to  write: 


Section  of  Plane  with  Overlapped  Tessellation. 


which  after  some  manipulation  becomes: 

/Or)  -  cos“'(x/2)  -  y  J  1-y- 

The  amount  of  area  within  the  triangle  that  is  covered  is  C(x)  «ir/2~3/(x)  and  subtracting 
this  from  the  area  of  the  triangle  gives: 

£(x)  —  ~~  x2  —  y-  +  3cos~'(x/2)  — 


Now  suppose  we  wanted  to  “partially”  tessellate  the  plane  with  unit  circles  in  such  a 
way  that  the  sum  of  the  overlapped  and  the  uncovered  areas  was  minimized  (thus  maximizing 
the  amount  of  singly  covered  area).  It  is  clear  from  figure  3.6,  that  lix)  is  increasing  in  x  and 
£(x)  is  decreasing,  and  thus  seeking  a  minimum  is  a  well  formed  problem.  Letting  the  objec¬ 
tive  function  be  Fix)  we  then  have: 

Fix)  -  -y-- x2-  ~  +  6cos_1(x/2)  -3x-^  1  -  — 

This  function  is  minimized  for  x*  — 1.9215  which  gives  the  tessellation  shown  in  figure  3.7  (the 
points  of  intersection  form  a  12-gon).  We  can  also  derive  the  fraction  of  the  plane  that  is 
covered  by  this  type  of  tessellation  by  forming  the  ratio  of  the  covered  area  to  the  area  of  the 
triangle.  When  this  is  done,  the  fraction  is  determined  to  be  .9278.  (as  compared  to  .9068  of 
equation  (3.1))  and  thus  about  92%  of  the  plane  is  singly  covered.  We  can  thus  write  the 
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FIGURE  3.7 

Allowing  Collisions  to  Increase  the  Bound. 


upper  bound  as: 

f'(N)  -  .9278/ A/-  (3.2) 

This  arrangement  of  circles  on  the  plane  was  conjectured  to  have  the  highest  density  of  singly 
covered  area  in  1964  [Toth64a]  and  has  recently  been  proved  by  Ernst  Strass  at  UCLA  (unpub¬ 
lished  manuscript).  We  should  observe  that  for  reasonably  large  A',  equations  (3.1)  and  (3.2) 
agree  for  practical  purposes  and  thus  all  the  approximations  and  conclusions  made  using  equa¬ 
tion  (3.1)  are  also  applicable  to  our  new  bound. 

If  we  let  <7  be  the  offered  load  (the  average  number  of  packets  presented  to  the  chan¬ 
nel  per  unit  time)  of  traffic  coming  from  the  area  of  one  circle,  we  clearly  see  that  (7  s  1.  This 
corroborates  our  results  of  the  previous  chapter  where  we  found  that  (7—1  for  the  slotted- 
ALOHA  network  analyzed  there.  As  mentioned  in  Chapter  2,  this  result  is  a  generalization  of 
that  for  a  single-hop  network  where  it  can  be  shown  that  an  offered  load  of  one  packet  per  unit 
time  optimizes  performance  [Abra70a],  We  can  derive  a  rule  of  thumb  for  multi-hop  packet 
radio  networks  by  approximating  equation  (3.2)  as: 


/(AO  =  1  /N 

It  can  best  be  interpreted  by  looking  at  figure  3.7.  Here  we  see  that  the  plane  can  also  be 
viewed  as  being  almost  tessellated  by  circles  of  unit  radius,  the  center  of  which  contains  a 


transmitter.  Since  each  transmitter  sends  to  an  average  of  N  other  nodes,  it  is  clear  that  only 
one  of  N  nodes  will  be  successful,  thus  giving  the  approximate  1//V  ratio. 


3.3  Comparison  to  slotted-ALOHA 

We  can  compare  this  optimal  fraction  of  transmissions  to  that  of  the  slotted-ALOHA 
radio  network  studied  in  the  previous  chapter.  In  that  chapter,  the  equation  for  the  probability 
of  successful  transmission  in  the  non-capture  case,  (which  corresponds  to  ourfi.VI)  was  given 
by: 

P(N,p)  -  (\-p)p(\-e-Nn)e~N"  (3.3) 

where  p  is  the  probability  of  transmitting  in  any  randomly  selected  slot.  Recall  the  (1— e~v/2) 
term  in  this  equation  is  a  factor  that  represents  the  probability  that  the  network  is  connected. 
Since  our  derivation  of/(/V)assumes  a  connected  graph,  we  should  eliminate  the  connectivity 
term  from  equation  (3.3)  before  comparing  the  performance  of  this  system  to  that  of  the 
optimal  protocol. 

For  a  given  N ,  P(N,p)  achieves  a  maximal  value  for: 

_•  \  N+2-VW+4 

p  &q(N) - — - 

We  can  thus  write  the  maximum  fraction  of  successful  transmissions  for  the  slotted-ALOHA 
network  as: 


P'(N)  -  q(N)(\-q(N))e~w 


(3.4) 


The  efficiency  of  a  protocol  is  defined  to  be  the  ratio  of  its  performance  to  that  of  the 
optimal  protocol  and  thus  we  define: 


e(N)  ±P'(N)/f’(N) 


FIGURE  3.8 

The  Efficiency  of  slotted-ALOHA. 


In  figure  3.8  we  have  plotted  P'(N)  and  e(N)  for  I^<V<11.  In  this  figure  we  observe 

that  e{N)  increases  in  N.  To  find  its  limit,  we  observe: 

lim  q(N)  -  1  /N 

V— oo 


hence  we  have: 


lim  e(N)  -  1/1.9278  e)  -  .3% 

,V— oo 

Thus  the  capacity  of  slotted-ALOHA  at  best  is  about  40®/o  of  the  channel  capacity  as  given  by 
our  approximation.  Observe  that  the  optimal  capacity  of  slotted-ALOHA  occurs  when  the  den¬ 
sity  of  terminals  approaches  infinity.  This  of  course  corresponds  to  an  unrealistic  network 
configuration.  The  increase  in  efficiency  as  grows  larger  is.  however,  very  slow.  For  realistic 
size  networks,  looking  at  figure  3.8.  we  see  that  we  can  approximate  the  efficiency  to  be  about 


Me.  In  a  single-hop  environment  slotted-ALOHA  achieves  a  maximum  capacity  of  Me  and 
thus  we  have  the  intuitively  pleasing  result  that  the  maximum  capacity  of  the  slotted-ALOHA 
protocol  is  about  Me  in  both  single  and  multi-hop  environments. 


3.4  Comparison  with  CSMA 

In  this  section  we  will  use  the  method  defined  in  the  previous  sections  to  derive  an 
equation  for  the  optimal  performance  for  CSMA  in  a  multi-hop  environment  and  then  will  cal¬ 
culate  the  efficiency  of  this  protocol.  In  a  CSMA  environment,  nodes  that  hear  an  idle  channel 
can  transmit  packets  on  the  channel.  In  a  random  connected  network  the  maximal  number  of 
nodes  that  could  transmit  using  the  CSMA  protocol  is  equal  to  the  size  of  the  maximal  2-order 
independent  set.  We  can  calculate  the  fraction  of  nodes  in  such  a  network  using  methods 
developed  in  previous  sections.  If  we  imagine,  in  figure  3.3,  the  length  of  the  equilateral  trian¬ 
gle  to  be  the  average  euclidean  distance  for  nodes  separated  by  2  hops,  defined  to  be  ?,  then 
g(N ),  the  fraction  of  nodes  in  a  maximal  2-order  independent  set,  can  be  calculated  as: 

g(N)  -  2.214 IN 

where  we  have  used  the  result  derived  in  Appendix  B  that  Y  -  1.2881.  If  we  assume  the  frac¬ 
tion  of  the  nodes  in  a  maximal  2-order  independent  set  that  transmit  is  p,  then  the  fraction  of 
transmissions  in  the  network  is  given  by  g(N)p.  Unlike  the  derivation  of  the  optimal  protocol 
however,  in  this  case  not  all  transmissions  are  successful.  We  can  calculate  the  probability  that 
a  randomly  selected  transmitter  is  successful  by  using  figure  3.9.  In  this  figure  we  have  shown 
the  triangle  tessellation  corresponding  to  a  maximal  2-order  independent  set,  and  the  circular 
tessellation  that  results  when  nodes  on  the  vertices  of  the  triangles  transmit.  If  we  suppose  in 
this  figure  that  P  is  a  transmitter,  then  we  can  calculate  the  probability  that  it  is  successful  by 
determining  the  probability  that  P  is  sending  its  message  to  a  neighboring  node  located  in  one 
of  the  areas  labeled  A,  B,  or  C  of  the  figure.  We  will  assume  a  uniform  probability  in  making 
this  calculation  and  also  assume  that  each  transmission  corresponds  to  a  series  of  messages  hav¬ 
ing  infinitesimal  length.  This  assumption  allows  us  to  calculate  the  expected  fraction  of  mes¬ 
sages  that  are  successfully  received  before  being  collided  by  an  interfering  transmission.  An 
alternative  way  to  model  this  is  to  assume  that  each  receiver  captures  the  first  signal  that  is  sent 
to  it.  Although  such  perfect  time  capture  [Davi80al  is  not  attainable  in  real  networks,  this 
assumption  will  allow  us  to  calculate  an  upper  bound  for  network  performance.  Area  C  +  B 
forms  half  of  the  intersection  of  two  unit  circles  separated  by  a  distance  of  ?,  and  can  be 
shown  to  be  equal  to  n/2  —  A(Y,  l)/2  -  .3845  where  A  (•)  is  derived  in  Appendix  B  (equation 
(B.7) ) .  The  area  of  the  triangle  is  easily  seen  to  be  .7094  and  thus  we  can  write  the  following 
equations: 


55 


FIGURE  3.9 

Tessellation  for  CSMA  Calculation. 

ir/6  —  .3845  —  A  +  B 

.7094  -  3  (A+B)  +  C 

which  can  be  solved  to  yield  A  —  .0467,  B  —  .0923,  and  C  -  .2921.  We  can  convert  these  areas 
into  probabilities  by  dividing  by  ir/6  to  obtain  PA  —  .0891,  />fl~.1762.  and  Pc  =  .5579.  Each 
area  is  influenced  by  only  a  subset  of  the  possible  transmitters  surrounding  it,  and  we  can  write 
the  probability  that  P's  transmission  is  successfully  received  as: 

Q  -  PA  +  2/>s(l-p)  +  />c(l-/>)2 

The  fraction  of  successful  traffic  then  is  given  by  C{N,p)  ~g(N)pQ  which  is  maximized  for 
p'-.4599  at  which  point  Q  -  .4421.  Using  this  in  the  expression  for  C(N.p)  allows  us  to  cal¬ 
culate  the  maximal  fraction  of  successful  transmissions  in  the  CSMA  network  which  is  given 

by: 

C*( AO  -  .4504/ /V 

Thus  we  have  the  efficiency  of  the  CSMA  protocol  in  the  multi-hop  environment  under  ideal 
assumptions  is  about  48.5%.  This  performance  is  in  striking  contrast  to  the  efficiency  of  CSMA 
in  the  single-hop  environment  [K!ei75a]  where  under  commonly  held  assumptions  it  has  a 


throughput  of  about  87%.  We  may  also  note  that  the  performance  of  CSMA  in  the  multi-hop 
environment  under  ideal  assumptions  is  not  much  greater  than  that  of  slotted- ALOHA. 


3.5  Densities  of  0,  R\,Ri  and 

As  a  by-product  of  this  research,  we  derive  in  Appendix  B  the  probability  density  func¬ 
tions  for  R i,  R 2,  R 3,  and  0  which  are  plotted  in  figures  3.10,  3.11,  and  3.12.  In  table  1  we 
give  the  means  and  variances  for  these  random  variables. 


FIGURE  3.10 
The  Density  of  R  2. 
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3.6  Conclusions 


In  this  chapter  we  have  calculated  an  upper  bound  for  the  maximal  expected  fraction  of 
successful  transmissions  obtained  in  a  random  connected  planar  packet  radio  network.  The 
method  utilized  is  notable  in  that  it  produced  an  intuitively  pleasing  result  that  /(AO  —  l/N 
and  was  simple,  although  tedious,  to  calculate.  Using  this  optimal  performance  as  a  standard, 
we  compared  the  slotted-ALOHA  and  CSMA  protocols  to  it,  and  have  shown  that  these  proto¬ 
cols,  for  realistic  networks,  have  an  efficiency  of  about  36%,  and  48.5%  respectively.  The  pro¬ 
bability  density  functions  for  the  euclidean  distance  for  nodes  separated  by  one,  two,  and  three 
hops  was  derived  as  a  by-product  of  this  research.  In  the  next  chapter  we  will  investigate  the 
delay  characteristics  of  the  optimal  protocol  which  was  discussed  in  this  chapter.  This  will  lead 
us  to  formulate  a  capacity  assignment  problem  for  networks  using  this  protocol. 


CHAPTER  4 
SPATIAL— TDMA 


In  the  previous  chapter  we  determined  the  maximal  throughput  for  a  network  which 
used  an  protocol  that  gave  a  maximum  number  of  nodes  in  the  network  permission  to  transmit 
over  all  periods  of  time  in  a  manner  that  did  not  allow  collisions.  In  this  chapter  we  formally 
define  this  protocol,  called  spatial-TDMA,  and  determine  its  delay  characteristics.  This  leads  us 
to  formulate  an  optimal  capacity  assignment  problem  for  networks  of  this  type  which  we 
address  in  a  later  section.  First  we  will  review  the  motivation  for  defining  such  a  protocol  and 
the  definition  of  the  capacity  assignment  problem. 


4.1  Introduction 

Networks  have  traditionally  consisted  of  a  set  of  switches  connected  together  by  some 
form  of  cable.  Because  cable  must  be  strung  point  to  point,  and  often  require  land  acquisition 
and  construction  of  supporting  structures,  the  design  of  such  networks  requires  careful,  long- 
range  planning  and  it  is  not  surprising  that  considerable  research  has  been  devoted  to  the  prob¬ 
lems  of  this  initial  design.  One  problem  that  has  been  extensively  studied,  for  example  arises 
when  the  locations  of  the  nodes  of  the  network  and  their  traffic  characteristics  are  assumed  to 
be  known.  In  this  case,  one  must  determine  a  procedure  for  choosing  the  capacities  of  the 
communication  lines.  Since  the  cost  of  the  network  is  an  increasing  function  of  these  capaci¬ 
ties,  it  is  important  that  the  designer  of  such  a  network  select  the  assignment  of  the  capacity  of 
the  links  of  the  network  to  minimize  network  expense  (under  the  constraints  of  preserving  a 
tolerable  delay  for  messages  in  the  system).  Various  solutions 
[Klei64a,Cant74a,  Fran71a,  Meis71a]  have  been  proposed  for  this,  and  other  problems  of  a 
similar  genre,  which  are  commonly  classified  as  Capacity  Assignment  Problems.  As  one  can 
expect  however,  network  specifications  often  change,  and  an  optimal  design  for  an  initial  net¬ 
work  may  be  far  from  optimal  after  changes  are  made  either  to  the  traffic  characteristics  of  the 
nodes  or  to  the  network’s  topological  structure.  Indeed,  even  the  minor  change  of  optimally 
adding  one  new  node  to  a  wire  network  can  be  a  formidable  problem.  Besides  the  difficulty  of 
connecting  cable  from  the  new  node  to  its  adjacent  neighbors,  one  must  also  solve  the  capacity 
assignment  problem  again  for  the  changed  network.  If  the  traffic  offered  by  the  new  node 
changes  the  loads  on  the  links  of  the  already  existing  network,  then  some  of  them  will  need  to 
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be  upgraded  to  higher  capacities.  Likewise,  other  lines  might  be  able  to  have  their  capacities 
reduced  if  some  of  the  flow  through  them  can  be  routed  over  the  new  lines.  Upgrading  exist¬ 
ing  lines  however,  is  expensive,  and  the  cost  of  preserving  this  optimality  in  the  network  after 
the  addition  of  a  new  node  could  be  so  formidable  that  the  designers  would  have  to  settle  for  a 
less  than  optimal  solution.  As  the  process  of  adding  new  nodes  to  the  net  continued,  it  would 
not  be  surprising  to  see  performance  seriously  deteriorate. 

To  avoid  these  problems,  one  needs  to  use  a  more  flexible  medium  for  inter-connecting 
the  switches  of  the  network.  As  we  saw  in  Chapter  1,  a  broadcast  medium  such  as  radio,  offers 
flexibility  to  topological  changes,  and  has  the  property  that  capacities  can  be  changed  to  reflect 
alterations  in  the  specifications  of  the  network.  Connectivity  in  such  networks  is  determined  by 
the  power  of  the  transmitters  of  the  network  and  thus  is  easily  controlled,  and  the  addition  of 
new  nodes  is  free  of  the  cost  of  stringing  cable.  There  still  remains  of  course,  the  problem  of 
re-adjusting  the  capacities  of  the  lines  of  the  system  in  an  optimal  manner.  Indeed,  determin¬ 
ing  the  capacity  of  a  specific  link  in  a  radio  net,  and  adjusting  the  capacities  of  all  the  nodes  in 
the  network  in  some  optimal  manner,  are  not  straightforward  problems  and  depend  upon  the 
protocol  that  nodes  use  to  determine  when  to  transmit  messages. 

We  address  the  above  issues  in  this  chapter  and  formulate  a  collision-free  channel 
access  protocol  that  nodes  use  to  determine  when  they  have  permission  to  transmit.  With  this 
collision- free  property,  a  multi-hop  packet  radio  network  can  be  made  to  simulate  a  wire  net¬ 
work,  and  the  determination  of  the  capacity  of  an  arc  is  straightforward.  The  basic  idea  of  this 
protocol  is  to  allocate  specific  periods,  from  a  specified  cycle  time,  during  which  nodes  can 
transmit.  This  allocation  has  the  property  that  during  any  period,  no  link  is  denied  transmission 
rights  if  it  would  not  potentially  cause  a  collision  (thus  this  is  similar  to  the  protocol  that  was 
discussed  in  Chapter  3).  The  capacity  of  an  arc  is  then  proportional  to  the  amount  of  time 
from  the  cycle  during  which  it  has  permission  to  transmit.  We  first  will  discuss  the  delay 
characteristics  of  such  a  protocol  and  then  propose  a  method  for  assigning  periods  and  adjusting 
their  durations. 


4.2  Description  of  the  Protocol 

Throughout  our  discussion  we  will  assume,  as  in  wire  networks,  that  the  locations  and 
connectivity  of  the  packet  radio  transceivers  and  the  expected  flow  of  traffic  between  any  two 
nodes  i  and  j,  yLj  are  known.  A  fixed  route  between  all  possible  source-destination  pairs  is 
also  assumed  to  be  given.  This  allows  us  to  calculate  the  flow  on  the  i ,h  arc.  a  of  the  net¬ 
work.  These  are  standard  assumptions  made  before  addressing  the  capacity  assignment  problem 
for  wire  networks.  There  are  however,  other  assumptions  we  will  make  that  are  particular  to 


the  radio  environment.  We  assume  the  channel  is  noise-free  and  that  receivers  have  the  ability 
to  capture  one  of  several  signals  on  the  channel.  Our  results  do  not  depend  upon  the  particular 
definition  of  capture  and  will  only  arise  in  determining  the  compatibility  matrix  discussed  in  the 
next  section.  We  will  also  assume  that  nodes  of  the  network  are  time  synchronized  so  that  all 
nodes  share  common  knowledge  about  period  beginning  and  ending  points  in  the  time  cycle, 
which  we  will  assume  has  a  length  of  T  time  units.  In  general,  since  distances  between  nodes 
of  the  network  vary,  time  synchronization  implies  that  each  slot  has  associated  with  it  a  guard 
band  that  is  used  to  assure  that  messages  lie  wholly  within  the  slot.  This  wasted  capacity  is  not 
taken  into  account  in  our  calculations. 

Knowing  the  locations  of  the  nodes  in  the  network  and  the  capture  parameter  of  the 
system,  allows  us  to  generate  a  compatibility  matrix.  In  such  a  matrix,  a  1  in  the  (i,j)  position 
indicates  that  arc  /  and  arc  j  of  the  network  can  simultaneously  transmit  (specifically,  their 
corresponding  nodes  being  able  to  transmit  in  these  directions)  without  causing  a  collision. 
Such  an  arc  is  said  to  be  enabled.  Using  this  compatibility  matrix  as  the  adjacency  matrix  of  a 
graph,  one  can  generate  a  set  of  cliques  containing  arcs  having  the  property  that  all  arcs  in  the 
same  clique  can  be  simultaneously  enabled  without  causing  a  collision.  If  we  let  C,  denote  the 
i'h  clique,  we  can  form  a  clique  cover  C  -  { C[,  C2,  •  •  •  ,  Ck )  with  the  property  that  every  arc  is 
contained  in  at  least  one  member  of  C.  For  each  clique  in  the  clique  cover  one  can  assign  an 
interval  of  time,  from  a  given  time  cycle  (i.e.  a  frame  whose  structure  repeats),  during  which 
arcs  in  that  clique  are  enabled.  Since  any  particular  arc  can  be  contained  in  more  than  one 
clique,  the  times  during  which  an  arc  is  enabled  depend  upon  the  times  assigned  to  the  cliques 
of  which  it  is  a  member. 

In  Figure  4.1  we  exemplify  this  construction.  In  this  figure,  the  six  nodes  (A  through 
E)  have  arcs  which  are  labeled  1  through  10.  When  node  B  enables  arc  7,  nodes  A ,  C,  and  E 
hear  the  transmission  and  B' s  message  is  addressed  to  node  E.  During  this  transmission,  if  we 
assume  that  receivers  cannot  capture  transmissions,  one  of  arcs  5  or  6  can  also  be  enabled 
without  causing  collisions.  In  the  compatibility  matrix  shown  in  the  figure,  then,  the  7'th  row 
contains  a  1  in  the  5’th,  6’th.  and  7’th  positions.  Each  row  of  this  compatibility  matrix  is  gen¬ 
erated  in  a  like  manner.  Using  this  matrix,  we  have  generated  all  possible  maximal  cliques 
C|,  C2,  •••  ,Cio,  and  from  these  have  selected  a  particular  clique  cover  consisting  of 
c-  (C,,C2,  •  •  •  ,C6). 

If  we  denote  the  cycle  time  of  the  frame  as  T  time  units,  we  can  allocate  time  durations 

& 

i,  (not  necessarily  contiguous),  where  T  -  £  t,  to  each  of  the  cliques  during  which  their  arcs 

1-1 

are  enabled.  Each  frame  then  consists  of  a  set  of  durations  during  which  cliques  are  enabled, 
and  the  sequence  of  the  frame  cycles  every  T  time  units. 
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FIGURE  4.1. 
Creation  of  the  Cliques. 


4.3  Queueing  Approximation 


In  this  section  we  will  describe  an  algorithm  for  approximating  the  average  system  time 
(queueing  plus  transmission)  of  messages  passing  through  a  single  node  using  spatial-TDMA. 
The  queueing  system  at  a  given  node  created  by  this  protocol  consists  of  non-overlapping  inter¬ 
nal  arrival,  service,  and  external  arrival  periods  which  are  enabled  during  specific  periods  of  the 
frame.  For  example  in  Figure  4.2  we  have  shown  a  queue  (contained  in  a  single  node  of  the 
network)  with  its  corresponding  frame.  Figure  4.3.  The  queue  contains  two  switches  si  and  si 
which  are  used  to  control  the  internal  input  and  service  processes.  At  most  one  of  these 
switches  can  be  closed  during  any  part  of  the  frame  since  we  do  not  allow  simultaneous  recep¬ 
tion  and  transmission  by  radios  in  the  network.  During  the  first  20  time  units  of  the  frame  in 
Figure  4.3,  we  see  an  internal  arrival  interval,  during  which  si  is  closed  (s;  will  be  open)  to  allow 
arrivals  to  enter  from  inside  the  network  (arrivals  from  the  attached  host  computer  called 
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FIGURE  4.2 

Model  of  Queue  for  Nodes  in  the  Network. 


TIME 


FIGURE  4.3 
A  Time  Frame. 

external  arrivals  can  also  occur  during  this  time).  These  arrivals  occur  at  rate  /,„  message  units 
(assumed  to  be  packets,  bits,  etc...)  per  time  unit  where  we  have  normalized  transmission  rates 
so  that  each  time  unit  can  accommodate  at  most  one  message  unit.  We  will  model  the  internal 
arrival  process  as  a  truncated  Poisson  process.  The  number  of  internal  arrivals  to  the  queue 
during  the  first  interval  then  is  distributed  as  a  Poisson  process  with  parameter  /,„  subject  to  the 


condition  that  no  more  than  20  message  units  arrive  during  this  interval.  In  general,  for  an 
internal  arrival  interval  of  m  time  units,  if  we  let  P[k  int \  m]  be  the  probability  that  k  internal 
arrivals  occur  in  m  time  units,  we  have: 


P[k  int\  m ) 


(m  lm)k  e~(m 

k\  f  (mlin)' 


k  -  0, 1,  •  ••  ,  m 


During  the  next  phase  of  the  frame,  a  service  interval,  of  10  units  in  which  S2  is  closed  (si  is 
open)  and  message  units  are  served  at  the  rate  of  one  per  time  unit.  During  this  time  a  max¬ 
imum  of  10  message  units  can  be  serviced.  The  next  phase  we  show  is  an  external  interval  dur¬ 
ing  which  both  S|  and  s 2  are  open  and  no  internal  arrivals  or  services  are  allowed.  In  the  figure 
we  have  labeled  the  section  of  the  frame  during  which  both  switches  are  open  as  being  idle  to 
indicate  that  no  new  messages  from  the  network  can  arrive  and  that  the  node  cannot  transmit. 
This  does  not  imply  however  that  the  are  no  active  processes  during  this  time.  Indeed,  over 
the  entire  frame  as  mentioned  before,  external  messages  (from  the  node's  attached  host)  can 
arrive.  As  shown  in  Figure  4.2  these  arrivals  immediately  enter  the  tail  of  the  queue  at  a  rate 
of  /„.  Again  this  arrival  process  is  assumed  to  have  Poisson  statistics  and  thus  the  probability 
that  k  external  arrivals  occur  in  m  time  units  is  given  by: 


P[k  ext\  ml 


(m  lex)k  e 


-  ( m  /_, ) 


/t! 


Switches  S!  and  s2  are  then  turned  on  and  off  according  to  the  time  patterns  depicted  in  Figure 
4.3  and  continue  to  cycle  every  100  time  units.  Observe  that  there  is  no  limit  on  number  of 
messages  that  can  arrive  from  the  attached  host  computer  over  the  course  of  a  frame  as  there  is 
from  messages  arriving  from  inside  the  network. 


An  exact  analysis  for  the  average  system  time  has  not  been  carried  out  as  of  yet  and 
thus  we  seek  an  approximate  formula.  Fortunately  a  fluid  approximation  [Klei76a]  to  this  sys¬ 
tem  gives  very  good  results  and  we  will  illustrate  this  method  using  Figure  4.4.  In  the  fluid 
approximation,  waiting  times  are  calculated  by  assuming  the  actual  backlog  of  packets  in  the 
queue  is  approximated  by  the  average  backlog. 

In  figure  4.4  we  have  plotted  the  growth  of  the  average  backlog  of  message  units  in  the 
system  during  the  course  of  the  frame  shown  in  figure  4.3.  For  reasons  to  be  explained  later, 
we  have  started  this  growth  pattern  at  the  beginning  of  the  last  external  interval  (i.e.  at  r-90). 
During  this  interval,  since  internal  arrivals  are  prohibited,  only  external  arrivals  can  add  to  the 
backlog.  If  we  assume  the  rate  at  which  the  queue  increases  during  an  external  period  to  be 
Mex,  we  see  that  the  expected  backlog  grows  linearly  with  this  slope  during  this  interval.  As 
the  frame  progresses  to  the  first  interval,  the  growth  rate  for  the  backlog  is  given  by  A/,„,  which 
is  the  sum  of  A/„  and  the  additional  rate  offered  by  the  internal  arrivals.  A  service  period  fol¬ 
lows  in  the  next  interval  and  the  backlog  drops  by  a  rate  of  A/,  <  0  message  units  per  time 
unit.  Since  external  messages  can  arrive  during  any  type  of  interval,  ,V/5  is  equal  to  Xfex  less 


FIGURE  4.4 

The  Backlog  for  Figure  4.3. 


the  service  rate  of  1  message  unit  per  unit  time.  As  seen  in  the  figure,  the  backlog  of  message 
units  in  the  queue  drops  to  zero  before  the  service  interval  is  finished.  Internal  arrivals  during 
this  time  are  assumed  to  be  immediately  serviced  and  thus  do  not  contribute  to  the  backlog. 
This  process  continues  in  this  manner  until  the  last  external  period,  at  which  point  it  starts 
from  a  zero  backlog  once  again.  We  will  call  a  point  on  the  frame  a  zero-point  (it  is  a  regenera¬ 
tive  point  for  the  frame)  if,  starting  with  an  empty  queue  at  this  point,  yields  a  backlog  of  zero 
after  exactly  one  frame. 

We  now  calculate  the  rates  described  in  the  previous  section.  First  we  make  the  follow¬ 
ing  definitions: 

Let  T  -  Length  of  the  time  frame 

T„  -  Total  time  of  external  intervals  per  frame 
T,„  —  Total  time  of  internal  arrival  intervals  per  frame 
Ts  -  Total  time  of  service  intervals  per  frame 


From  our  previous  discussion  we  can  write: 


(w  “  LexJ  T 
I  m  ”  L  ij  Tn 
Mex  ™  lex 

Mm  “  lin  "F  Mex 
M,  -  M(X  -  1 

The  area  under  the  backlog  curve  represents  the  number  of  message-time  units  accumulated  by 
messages  arriving  during  the  frame.  Dividing  this  by  the  average  number  of  message  units  that 
arrive  during  the  frame,  Lm  +  L„,  gives,  by  Little's  result,  [Litt61a]  the  average  time  spent  in 
the  queueing  system.  Because  we  began  the  calculation  at  a  zero-point,  message  delays  for  all 
arrivals  to  the  frame  are  counted.  It  is  clear  we  can  always  find  a  zero-point  on  any  frame  satis¬ 
fying  Ts>  L,„+  Lex. 

This  then  describes  the  algorithm  for  the  approximation  which  we  summarize  as: 

1.  Find  a  zero-point. 

2.  Calculate  the  area  under  the  backlog  curve. 

3.  Divide  by  Lin  +  to  arrive  at  the  average  system  time. 

It  is  clear  that  the  ordering  of  the  intervals  and  their  lengths  greatly  influences  the  aver¬ 
age  system  time.  Suppose,  for  example,  that  we  change  the  frame  in  Figure  4.3  by  coalescing 
all  the  service  and  intervals  together  and  placing  them  on  the  frame  as  shown  in  Figure  4.5. 
Although  it  has  exactly  the  same  interval  lengths,  this  translation  increases  the  average  system 
time  (in  fact  it  is  a  worst  case  example).  In  this  system,  any  messages  that  arrive  during  the 
internal  arrival  interval  must  wait  at  least  30  time  units  before  being  serviced.  If  the  service 
and  internal  arrival  intervals  are  interchanged,  it  is  clear  that  the  average  system  time  would 
decrease  by  at  least  this  much.  In  fact,  the  system  that  has  the  minimum  average  system  time 
is  one  that  spreads  arrival  and  service  intervals  in  infinitesimal  units  that  alternate  with  each 
other.  In  this  way,  an  arrival  is  immediately  serviced  in  the  following  interval  after  accruing  lit¬ 
tle  waiting  time.  In  a  practical  implementation  of  spatial-TDMA,  however,  there  are  limits  to 
how  small  one  can  make  interval  sizes  since  radios  have  a  finite  switching  time  between 
transmission  and  reception,  and  choosing  a  frame  pattern  that  minimizes  the  system  time  for 
one  particular  queue  in  the  network  does  not  in  general  decrease  the  total  average  message 
delay  for  messages  in  the  entire  network.  In  fact,  finding  the  alternation  that  does  achieve  the 
minimum  delay  for  all  message  in  the  network  is  a  very  difficult  problem  which  we  will  address 
in  a  later  section. 


FIGURE  4.5 
A  Worse  Case  Cycle. 


4.4  Discussion  of  Delay  Results 


In  this  section  we  compare  simulation  results  with  those  found  using  the  fluid  approxi¬ 
mation  described  in  the  previous  section.  Numerous  frames  were  randomly  generated  that  had 
the  same  input  parameters  (7,  7S,  Tm,  T„,  L„,  and  LJ  and  results  of  the  simulation  checked 
against  those  of  the  approximation.  Three  such  frames  are  shown  in  Figure  4.6  where  +1  steps 
correspond  to  internal  arrival  intervals,  0  to  external  intervals,  and  -1  to  service  intervals.  In 
all  three  frames  7-10000,  7„-6000,  7m-2000,  7,  -  2000,  the  average  service  and  internal 
arrival  intervals  have  length  200,  and  the  average  external  arrival  rate  is  200  message  units  over 
the  frame  time  7.  The  mean  service  time  as  a  function  of  p  for  these  frames  is  shown  in  Fig¬ 
ure  4.7. 


There  are  several  interesting  features  of  these  curves.  We  first  see  the  close  match 
between  the  simulation  and  approximation  thus  assuring  us  that  the  approximation  is  accurate. 
The  variation  between  the  mean  system  time  for  the  three  frames  is  very  large  which  shows  the 
dependency  upon  the  ordering  and  size  of  the  intervals  of  the  frame.  For  example,  for  p-.7, 
set  1  has  a  mean  system  delay  of  about  650  time  units  whereas  set  3  has  a  value  of  2800.  more 
than  four  times  as  much.  The  extreme  delays  of  the  third  frame  arise  from  the  long  periods 
(9000, 10000)  and  (0,2000),  during  which  there  are  no  service  intervals.  All  packets  arriving 
during  these  periods  create  a  backlog  that  cannot  be  depleted  until  much  later  in  the  frame.  On 
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FIGURE  4.6 

Three  Randomly  Generated  Frames. 


FIGURE  4.7 

Curves  for  the  Three  Frames  of  Figure  4.6. 

the  other  hand,  the  fortuitous  placement  of  intervals  in  frame  1  consists  of  groups  of  arrival 
intervals  followed  by  service  intervals  that  allow  an  accumulated  backlog  to  be  serviced  quickly. 

Even  though  there  is  a  large  variance  in  the  curves,  there  is  a  striking  similarity  in  the 
shapes  of  the  curves.  The  curves  are  very  well  approximated  by  a  piece-wise  linear  function 
(shown  as  a  dashed  line  in  the  figure).  We  note  that  they  do  not  approach  infinity  at  p-1. 
Since  other  queueing  systems  have  explosive  growth  as  p— *  1,  arising  from  the  omni-present 
1/(1— p)  term  in  most  queueing  equations,  this  behavior  is  very  unusual.  The  explosive 
growth  for  most  queueing  systems  arises,  however,  from  the  fact  that  there  is  variation  in  the 
arrival  pattern  of  messages  to  the  queue,  and  as  p—1  the  probability  that  a  sequence  of  arrivals 
saturates  the  queue  approaches  1.  In  spatial-TDMA  however,  the  number  of  arrivals  for  each 
internal  arrival  interval  is  constrained  to  be  no  greater  than  the  normalized  interval  length,  and 
thus  the  variance  of  arrival  statistics  is  also  constrained.  Besides  reducing  queueing  delays,  this 
restriction  on  the  number  of  messages  that  can  enter  the  system  during  an  interval  also  acts  as 


a  natural  flow  control  mechanism  for  messages  in  such  a  network. 


The  slope  changes  in  the  piece-wise  linear  approximation  occur  when  the  arrival  rate  is 
so  large  that  a  sufficient  number  of  arrivals  to  an  internal  arrival  interval  cannot  be  serviced  in 
the  next  set  of  service  epochs.  For  example,  the  change  about  the  point  p  =  .6  for  frame  2 
arises  from  the  fact  that  for  p  >  .6,  some  arrivals  over  the  interval  (6000,7800)  must  wait  until 
the  next  service  set  of  intervals  (1200,2300)  to  be  processed.  For  lesser  values  of  p,  p  <  .6, 
most  arrivals  to  (6000,7800)  are  serviced  in  the  interval  (7800,9100)  and  thus  suffer  less 
delay.  Naturally  as  p  increases,  the  proportion  of  messages  that  must  wait  until  (1200,2300)  to 
be  serviced  grows  and  so  does  the  mean  system  time.  Each  of  the  breaks  in  the  piece-wise 
linear  approximation  can  be  explained  in  this  manner. 


FIGURE  4.8 

Reducing  the  Average  Period  Size. 


In  Figure  4.8  we  have  shown  a  frame  for  the  same  input  parameters  but  where  the 
average  size  of  the  service  and  internal  arrival  intervals  is  equal  to  20  time  units  instead  of  200 
as  in  those  of  Figure  4.6.  The  corresponding  mean  system  time  curve  is  shown  in  Figure  4.9. 
We  see  a  marked  decrease  in  the  mean  system  time  for  this  frame  in  comparison  to  the  previ¬ 
ous  set  of  frames.  This  demonstrates  the  dependency  of  the  mean  system  time  upon  the  size 
of  the  intervals.  If  we  adjust  the  frame  to  minimize  the  mean  system  time,  as  shown  in  Figure 
4.10  (where  for  illustrative  clarity  we  have  only  shown  a  portion  of  the  frame),  the  resultant 
delay  is  approximately  equal  to  1  time  unit  throughout  the  entire  range  of  p.  For  such  a  frame 
the  najority  of  the  arrivals  to  the  system  are  immediately  serviced  in  the  following  service 
interval. 
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FIGURE  4.9 

The  Curve  for  Figure  4.8. 


4.5  The  Capacity  Assignment  Problem 


4.5.1  Introduction 

In  the  previous  sections  we  determined,  for  a  given  frame,  a  good  approximation  for 
the  average  delay  messages  experience  in  passing  through  a  network  using  spatial-TDMA.  We 
observed  that  the  ordering,  size,  and  number  of  periods  from  the  Irame  allocated  to  internal 
and  external  arrivals  and  service  periods,  had  a  great  influence  on  the  mean  delay.  One  obser¬ 
vation  we  made  was  that,  in  general,  the  mean  delay  decreased  as  the  average  size  of  the 
periods  decreased.  In  practical  implementations  of  this  protocol  there  is  a  minimum  size  for 
any  period  which  is  determined  by  the  switching  time  for  radio  is  in  the  network  and  by  over¬ 
head  considerations.  Let  us  therefore  suppose  that  all  periods  in  a  frame  have  a  length 
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FIGURE  4.10 
An  Optimal  Frame. 

determined  by  this  minimal  switching  time,  and  let  us  normalize  the  frame  duration.  T.  to  this 
slot  size.  We  will  say  that  one  message  unit  can  be  transmitted  during  this  duration.  With  this 
restriction,  then,  every  frame  consists  of  a  number  of  unit  slots  assigned  to  the  cliques  of  the 
clique  cover  in  some  particular  permutation.  This  restriction  does  not  effect  the  optimality  of 
frame  selection  to  any  great  degree,  since  the  unit  slot  is.  in  the  continuous  frame  case,  a  lower 
bound  to  slot  size  anyway,  and  is  usually  very  small  in  relationship  to  T.  We  will  therefore,  in 
a  later  section,  take  the  liberty  to  refer  to  the  size  of  the  intervals  Tex,  T,„,  and  T,,  as  being  .on- 
tinuous  parameters  when  we  pose  a  convex  programming  problem. 

Before  beginning  to  develop  the  capacity  assignment  problem,  let  us  first  establish  our 
notation.  Suppose  we  are  given  a  clique  cover,  C  -  {C1.C2,  •  •  ,  Ql,  and  a  time  vector. 
f2i  •  •  •  ,4),  where  t,  is  the  time  from  the  frame  (having  duration  T)  allocated  to  C . 

We  will  denote  1 1  f  1 1  -  X  r,  and  will  say  a  time  vector  j_  is  feasible  if  _r>0.  |  U  1 1  <  7\  and  cer- 

,-i 

tain  flow  constraints,  which  are  defined  later,  are  satisfied.  Let  the  total  flow  of  messages  into 
the  network  per  frame  be  given  by  y  and  let  a,  be  the  flow  of  messages  into  queue  i  over  that 
time.  Let  q  be  the  number  of  queues  in  the  network.  For  a  given  feasible  time  vector,  K  and 
queue  /',  let  D,U)  be  the  average  delay  of  messages  passing  through  queue  /,  and  let  T'n.  T'„ , 
and  T’  be  the  total  amount  of  time  from  the  cycle  allocated  for  the  internal,  external,  and  ser¬ 
vice  periods  respectively.  Also  for  queue  /  let  L',„  and  L‘ex  be  the  average  traffic  flow  from 
internal  and  external  sources  over  a  time  frame.  It  is  clear  that  a .**(£'„  +  £,.,.)/ T  and  we  will 
define  the  (k.vl)  binary  vectors  M',„,  Mm-  anc*  Ms  t0  satisfy  T'„~M'exl.  and 

T\~MsL-  Observe  that  we  can  form  a  matrix  M'  from  these  vectors  that  will  satislv 

t;„,  T[). 


With  these  definitions  we  can  define  the  capacity  assignment  problem  formally  as: 
Minimize  V  —  D,(t) 

l.  All  Possible  Frames  (_|  y 


subject  to  i  ^  0 
11/11-7-^0 

0  /  — 1,2,  •  ■  •  ,q 
0  i-l,2,  •  ■  •  ,<7 

Before  discussing  the  complexity  (and  mathematical  intractability)  of  solving  this  prob¬ 
lem,  let  us  first  establish  that  the  feasibility  region  is  convex.  Suppose  that  j_  and  £  are  two 
feasible  time  vectors.  To  show  convexity  we  need  to  show  that  pt  +  (\—p)£  where  0</><l  is 
also  feasible.  This  is  obvious  however,  since,  using  the  first  constraint: 

The  same  argument  holds  for  the  second  constraint. 

The  complexity  of  solving  this  capacity  assignment  problem  arises  from  the  fact  that  the 
optimization  occurs  over  all  feasible  time  vectors,  /,  and  all  possible  ways  of  allocating  periods 
in  frames  satisfying^.  Even  if  we  were  given  (by  some  divine  mathematician!)  the  time  vector 
which  was  optimal,  finding  the  frame  (the  ordering  of  all  T  slots)  that  achieves  minimal  delay 
is  in  itself  mathematically  intractable.  In  the  following  sections  we  will  present  an  approach  that 
finds  optimal  time  vectors  over  the  class  of  frames  that  have  been  randomly  selected,  but  first 
we  will  address  the  question  of  feasibility. 


4.5.2  Feasibility 


i 


In  this  section  we  establish  a  procedure  for  determining  if  a  given  clique  cover,  for  a 
given  set  of  queues  and  corresponding  flows  o,,  permits  a  feasible  time  vector.  It  is  easy  to 
create  cases  where,  no  matter  how  one  adjusts  the  components  of  a  time  vector,  some  queues 
will  have  more  flow  into  them  than  they  can  accommodate.  We  have  already  seen  that  the 
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feasible  region  is  convex.  This  allows  us  then  to  formulate  the  feasibility  problem  as  a  linear 
program: 

Minimize  1 1  _r  1 1 

l 

subject  to  i  >  0 
Min-1- L'm>  0  /  =  1,2,  •  •  • ,? 

Mil-  L\n  -  Ux>  0  /-1,2,  •  •  •  ,q 

If  the  solution  to  this  linear  program  has  IU II  <  T  than  any  vector  with  _r  +  e  where  e>0  and 
Nell  —  T—  1 1  1 1  is  a  feasible  time  vector,  otherwise  no  feasible  time  vector  exists.  We  might 
add  that  phase  1  of  solving  this  linear  program  can  easily  be  solved  by  a  time  vector  having  all 
components  equal  to  r;— Max  ( L  +  L^x) ,  and  that  phase  2  of  the  program  can  halt  as  soon  as 

i 

the  objective  function  1 1  HI  becomes  less  than  T.  Thus  this  feasibility  program  can  be 
efficiently  solved. 


4.5.3  Average  of  All  Randomly  Generated  Frames 

Due  to  the  mathematical  intractabilities  of  the  frame  assignment  problem,  we  seek  in 
this  section  to  create  an  approximation  to  the  average  delay  of  messages  in  passing  through  a 
queue  in  which  the  alterations  of  frame  periods  have  been  randomly  selected.  This  will 
represent  an  upper  bound  to  the  mean  delay  obtained  with  an  optimal  ordering  and,  as  we  will 
see,  will  lead  to  a  tractable  program.  Let  us  therefore  intuitively  motivate  our  approximation  to 
this  delay.  Recall  that  we  have  normalized  the  time  axis  into  time  slots  equal  to  the  transmis¬ 
sion  of  one  message  unit.  Suppose  we  are  given  T,n,  and  T,  and  have  selected  one  of  the 
possible  random  frames  having  these  time  values.  Select  a  point  in  this  frame,  T  where  T<  T 
and  let  Tex,  Tm ,  and  T,  be  the  number  of  external,  internal,  and  service  slots  used  in  generating 
the  frame  up  to  time  T.  Since  the  frame  was  randomly  generated,  the  next  slot.  T  + 1,  will  be 
drawn  with  a  uniform  probability  from  the  remaining  population  of  slots.  For  the  purpose  of 
clarity,  we  will  focus  our  discussion  in  terms  of  service  slots  (analogous  statements  can  be  made 
for  external  and  internal  arrival  slots).  At  the  F  +  l  step,  the  probability  that  the  next  slot 
selected  is  a  service  period  is  given  by  Pf  -  (f,— f,)/(r-  f,-  T„-  T,„ ).  The  original  proba¬ 
bility  of  selecting  a  service  interval  was  given  by  P,  «  TJT.  and  thus  we  see  this  probability 
changes  as  the  cycle  is  generated.  If  T  is  large  howevf,  over  most  of  the  frame,  one  would 
expect  that  the  probability  of  selecting  a  service  period  at  each  step  does  not  differ  much  from 
/*,.  Intuitively,  T ,,  Tm.  and  Tex  decrease  at  a  rate  in  proportion  to  their  sizes,  and  since  their 
sizes  are  assumed  to  be  large,  their  relative  proportions  do  not  alter  significantly  over  most  of 
the  frame.  Only  when  most  of  the  slots  have  been  given  out.  when  T  is  close  to  T.  will  we 
expect  a  wide  variation  in  P,.  The  main  point  is  that  until  this  time,  Ps  has  remained  close  to 


P,  and  the  major  component  of  the  deiay,  as  determined  by  the  backlog  curve  of  the  fluid 
approximation,  has  already  been  determined.  We  thus  conclude  that  the  delay  obtaineu  from 
generating  frames  using  the  initial  probabilities,  P5.  P„.  and  Pex  throughout  the  frame  genera¬ 
tion  will  not  differ  much  from  those  generated  by  allowing  Ps  to  change  throughout  frame  gen¬ 
eration. 

Now  that  we  have  motivated  the  approximation,  let  us  describe  it  more  precisely.  Sup¬ 
pose  we  are  at  a  height  k  on  the  backlog  curve.  From  the  previous  discussion,  the  probability 
that  the  next  slot  is  an  external,  internal,  or  service  slot  is  given  by  Pex=T„!T.  P,„*=TJT ,  and 
Ps=*T$/T  respectively.  Recall  that  !m=LjT  is  the  probability  that  a  randomly  selected  slot 
contains  an  internal  arrival.  Define  u,  /~1,2  to  be  the  probability  that  /  arrivals  come  into  the 
system  in  the  next  slot,  and  let  d  be  the  probability  that  a  message  leaves  the  system  in  the 
next  slot.  Using  the  above  definitions  we  can  then  write: 

«1  -  Pex  lex  +  Pm  I  4*0  ~  O  +40“  4r>  1 

M  2  ™  P nt  l/n  lex 

d  =  PM- U 

and  the  evolution  of  the  height  of  the  backlog  curve  is  given  by  a  random  walk  with  the  above 
probabilities.  The  state  transition  graph  for  this  random  walk  is  given  in  figure  4.11. 


FIGURE  4.11 

State  Transitions  For  the  Random  Walk. 

The  state  transition  equations  are  given  by: 

(u]  +  u2)P0-  dP}  (4.1) 

(u\  +  u2  +  d)P\  »  U\Po  +■  dP2 
(«i  +  u2  +  d)Pk  -  u2Pk- 2  +  u\Pk.\  +  dPk~\  02 

where  Pk  is  the  probability  that  the  system  is  in  state  k.  Before  solving  these  equations,  let  us 
determine  how  we  will  use  the  solution  to  generate  our  approximation  to  the  average  delay. 
We  will  proceed  in  the  same  way  as  in  the  fluid  approximation,  namely  finding  the  area  under 
the  backlog  curve  and  then  dividing  by  the  average  number  of  customers  that  entered  the  sys¬ 
tem.  Suppose  we  knew  the  steady  state  probabilities  Pk  of  being  at  step  k  in  our  approxima¬ 
tion.  Points  at  a  height  k  that  move  to  £  + 1  incur  an  increase  in  the  area  that  is  equal  to 
k  +  . 5  where  the  k  arises  from  the  rectangle  of  height  k  and  width  1,  and  the  .5.  from  the 
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triangle  having  unit  width  and  height.  Similar  calculations  can  be  made  for  other  steps  and  we 
can  write  the  average  area  for  each  step  as: 

A  —  Pq(ut  +  .5 m i )  +  £  P, {u\(i  +  .5)  +  u2(i  + 1  )+</(/  —  .5)1  (4.2) 

i-i 

—  Po(u2  +  .5  U\)  +  P(u\  +  u2  +  d)  +  (1  —  P0)(u2  +  -5u\  —  .5  d) 


In  these  equations  P  is  the  average  value  of  k.  We  can  calculate  both  P  and  Pn  by  using  stan¬ 
dard  transform  methods  on  equations  (4.1)  [Klei76aJ.  Defining  the  z-transform  of  P  to  be 

P(z)-£  P, r'.  we  have: 

-o 

/>(-)  .  2u2  +  ui-d 

u2:2  +  (u\  +  u2):~  d 

from  which,  using  standard  theory  we  can  derive: 

-  _  Su2  +  ut 
d  —  2u2- u\ 
n  d  —  2u2—  u\ 
p0 - j - 

and  using  these  in  equation  (4.2)  yields: 

A  — P(u\  +  u2  +  d )  (4.3) 

Thus  we  can  calculate  the  average  delay  by  dividing  by  the  average  number  of  messages  that 
arrive  in  a  randomly  selected  slot  to  determine  the  following  delay  equation: 


D(T„.  T,„,  r,) 


T  (3mi  +  wi) 

(Lex  +  L,J  (d  —  2u2  —  u\) 


(U\  +  Mi  +  d) 


We  have  denoted  D(  )  as  being  a  function  of  three  variables  for  the  purpose  of  clarity.  Since 
r-r,B  +  T„  +  T,  only  two  of  the  variables  are  independent. 


4.5.5  Comparing  the  Approximation  with  Randomly  Generated  Frames 

In  this  section  we  want  to  show  the  results  when  equation  (4.4)  is  compared  with  data 
obtained  when  we  generate  random  frames  and  calculate  their  mean  delay  using  the  fluid 
approximation.  Our  procedure  was.  for  a  given  set  of  parameter  values  L„.  L  ,r  T.x.  T and 
7*,,  to  generate  1000  random  frames  and  determine  the  minimum,  maximum,  mean,  and  vari¬ 
ance  of  the  delays  for  these  samples.  We  then  compared  the  average  delay  to  equation  M.4>  to 
determine  if  the  approximation  was  close.  We  performed  this  procedure  for  many  different 
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selections  of  parameter  values  and  all  showed  similar  behavior.  To  demonstrate  this  behavior 
we  have  selected  two  sets  to  plot.  Each  of  these  sets  had  T-1000.  Set  1  has  Tm- 400. 
r«-100,  7>500,  I„-100,  and  we  varied  L,„  over  the  range  2(KZ.,„<380.  Set  2  had 
r„-300,  7^-200,  7>500,  Im-200  and  L„  varied  over  20<IM.<280.  In  figure  4.12  we  show 
how  the  approximation  fared  in  relationship  to  the  mean  of  the  generated  frames. 


From  GENERATED  CYCLES 
From  APPROXIMATION 


FIGURE  4.12 

Validating  the  Approximation. 

We  see  in  this  figure  that  the  approximation  is  very  close  to  that  of  the  sample  mean  and  that 
only  for  very  large  values  of  p  does  it  break  away  from  the  generated  frames.  This  is  explained 
by  the  fact  that  the  approximation  is  "more  stochastic"  than  the  randomly  generated  frames 
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since  it  is  given  by  a  random  walk,  and  that  for  high  p  values  the  random  walk  tends  to  spend  a 
lot  of  time  "attempting  to  escape  to  infinity"  in  contrast  to  the  generated  frames  that  have  an 
upper  bound  to  their  maximum  height. 


We  can  extract  more  information  about  the  randomly  generated  frames  by  looking  at 
Table  4.1. 


p 

Mean 

Variance 

Min  Delay 

Max  Delay 

0.240000 

1.019933 

0.059683 

0.882035 

1.292519 

0.280000 

1.173994 

0.076922 

1.000477 

1.418889 

0.320000 

1.345148 

0.097195 

1.122492 

1.737658 

0.360000 

1.519549 

0.123856 

1.228147 

2.035303 

0.400000 

1.723441 

0.145813 

1.376908 

2.425304 

0.440000 

1.919901 

0.193245 

1.435101 

2.563087 

0.480000 

2.158198 

0.218624 

1.594046 

3.138185 

0.520000 

2.434490 

0.241395 

1.902604 

3.434316 

0.560000 

2.769812 

0.325750 

1.993773 

^  291 554 

0.600000 

3.172192 

0.414819 

2.390554 

4./3C672 

0.640000 

3.621928 

0.522270 

2.618372 

5.937060 

0.680000 

4.222784 

0.670217 

2.910517 

7.552518 

0.720000 

4.881158 

0.869889 

3.210058 

9.927224 

0.760000 

5.741120 

1.051270 

3.566030 

9.822580 

0.800000 

7.005569 

1.534562 

3.751551 

15.805325 

0.840000 

8.766911 

2.257944 

5.079778 

19.209290 

0.880000 

11.501342 

3.460731 

6.268044 

32.473289 

0.920000 

15.094280 

.  4.126865 

7.613554 

34.197002 

0.960000 

22.544538 

7.200002 

9.852835 

65.090919 

TABLE  4.1. 

In  this  table  we  list,  for  given  values  of  p,  for  set  1,  the  mean  delay,  variance,  minimum  delay, 
and  maximum  delay  that  was  found  over  the  generated  frames.  We  see  that  although  the 
difference  between  the  minimum  and  maximum  delays  are  often  quite  substantial  (especially 
for  large  p)  the  variance  is  usually  very  small.  This  implies  that  the  mean  of  all  the  generated 
frames  is  not  much  different  from  the  mean  of  a  particular  generated  frame.  Another  way  to 
see  this  is  to  plot  the  coefficient  of  variation  (defined  to  be  the  variance  divided  by  the  mean) 
as  a  function  of  p.  In  figure  4.13  we  plot  this,  as  well  as  the  variance,  and  see  that  the 
coefficient  of  variation  is  quite  small  throughout  all  ranges  of  p.  Using  this  approximation  we 
are  now  in  a  position  to  formulate  the  capacity  assignment  problem  for  networks  of  this  type. 


FIGURE  4.13 

The  Coefficient  of  Variation  as  a  Function  of  p. 


4.5.S  The  Capacity  Assignment  Problem  for  Random  Frames 

In  this  section  we  will  create  the  capacity  assignment  problem  for  spatial-TDMA  net¬ 
works  in  which  the  time  slots  to  cliques  have  been  performed  in  a  random  manner.  We  first 


suppose  that  there  exists  a  feasible  time  vector.  Knowing  this  we  can  then  formulate  the  capa¬ 
city  assignment  problem  as: 

Minimize  f  —  Dj(M'-t) 
i  Iti  y 

subject  to  i  >  0 

llitl-r<o 

Mln-i-L'm>0  /-1,2,  ••  •  ,q 

0  /-l, 2.  •••,<? 

Once  again  the  feasibility  region  is  convex  and  it  can  be  also  shown,  by  differentiating  D,(  ), 
that  the  Hessian  matrix  of  second  partials  is  positive  semi-definite  and  thus  the  objective  func¬ 
tion  is  also  convex.  We  thus  have  a  convex  programming  problem  which  can  be  solved  using 
any  number  of  well-known  solution  techniques  (Avri76al. 

4.6  Conclusions 

The  flexibility  offered  by  broadcast  radio  presents  clear  advantages  over  the  use,  in  con¬ 
ventional  networks,  of  point  to  point  cable.  Besides  being  more  adaptable  to  changing  topolo¬ 
gies  (  and  as  demonstrated  in  this  paper,  changing  traffic  characteristics)  radio  does  not  require 
the  construction  of  elaborate  supporting  mechanisms  for  their  implementation  (this  is  not  to 
imply  that  packet  radio  networks  do  not  have  their  technical  problems!).  The  design  of  net¬ 
works  that  use  spatial-TDMA  is  facilitated  since  precise  specifications  of  topologies  and  traffic 
characteristics  do  not  initially  need  to  be  known.  Very  often  the  true  flow  of  traffic  between 
nodes  of  a  network  are  only  known  after  the  network  has  been  used  for  a  period  of  time.  As 
shown  before,  incorrect  estimations  of  the  traffic  flows  can  be  disastrous  for  a  wire  network. 
With  radio  however,  network  designers  have  the  ability  to  simply  alter  capacities  to  match  the 
actual  traffic  flows.  A  straightforward  approach  based  on  the  results  of  this  paper,  would  be  to 
start  with  an  initial  guess  of  the  traffic  flows.  These  would  be  used  to  determine  the  initial 
capacities,  and  the  system  would  be  allowed  to  operate  for  a  period  of  time  during  which  statis¬ 
tics  of  the  actual  traffic  flows  would  be  gathered.  These  more  precise  estimations  would  then  be 
used  to  again  assign  capacities.  This  process  could  be  refined  further,  where  perhaps  several 
capacity  assignment  plans  could  be  incorporated  to  match  the  flows  during  different  parts  of  the 
day,  of  the  week,  of  the  year  etc...  since  traffic  patterns  frequently  change  according  to  a 
predictable  schedule  (for  example  there  is  very  little  traffic  flowing  from  the  east  coast  after  5 
p.m.).  Scheduling  and  synchronizing  changes  between  different  plans,  would  be  a  non-trivial 
implementation  problem. 

The  above  advantages  of  radio  based  networks  can  also  be  combined  with  the  results  of 
the  extensive  research  into  the  behavior  of  wire  networks.  Since  with  spatial-TDMA.  radio  net¬ 
works  simulate  wire  networks,  many  of  these  results  can  be  readily  adapted  to  the  radio 
environment.  For  instance,  routing  and  flow  control  techniques  designed  for  wire  networks  can 


be  run  without  modifications  on  the  radio  net.  In  particular,  the  flow-deviation  algorithm 
[Frat73a]  could  be  used  to  calculate  an  optimal  (minimal  average  delay)  routing  policy.  In  this 
procedure,  flows  on  the  arcs  are  altered  (preserving  the  flow  required  between  source- 
destination  nodes  of  the  net)  in  a  fashion  that  lowers  the  average  delay  until  a  minimum  is 
found.  Since  the  function  for  the  average  delay  is  different  in  our  case  than  that  of  the  stan¬ 
dard  flow-deviation  algorithm,  we  would  have  to  modify  the  algorithm  slightly  but  the  overall 
strategy  would  be  similar  since  the  equation  for  the  mean  delay  is  convex.  After  executing  this 
algorithm,  the  flows  on  the  arcs  would  be  changed,  and  thus  we  could  run  the  capacity  assign¬ 
ment  algorithm  specified  in  this  paper  to  re-assign  capacities  to  the  arcs.  After  this  new  capacity 
was  assigned,  it  would  be  followed  by  the  flow-deviation  algorithm  again  and  this  alteration 
would  continue  until  the  change  in  the  average  delay  became  smaller  than  some  specified  value. 
This  Capacity  and  Flow  Assignment  problem  has  been  developed  for  wire  networks  [Klei76a]  and 
one  could  adapt  similar  procedures  to  be  used  in  the  radio  environment. 

Thus  it  would  appear  that  radio  networks  using  spatiai-TDMA  have  both  the  advantages 
of  a  flexible  topology  and  that  of  being  able  to  draw  from  the  wealth  of  results  created  for  con¬ 
ventional  wire  networks.  It  also  however  suffers  from  some  of  the  disadvantages  of  wire  based 
networks.  In  both  spatial-TDMA  and  wire  networks,  nodes  are  assumed  to  be  stationary  and 
besides  help  from  adaptive  routing  procedures,  there  is  no  allowance  for  dynamic  changes  to 
the  capacities  of  the  arcs  of  the  network.  Indeed  it  would  be  somewhat  tragic  if  the  most 
salient  feature  of  the  radio  media,  its  ability  to  adapt  to  topological  changes,  was  not  taken  into 
account  in  the  design  of  the  channel  access  scheme  of  a  packet  radio  network.  The  channel 
access  scheme  for  such  a  network  would  need  to  account  for  changing  topologies  and  traffic 
demands.  Certainly  the  access  scheme  as  presented  in  this  paper  does  not  have  these  charac¬ 
teristics.  However  by  modifying  the  scheme  in  a  manner  which  allows  dynamic  slot  adapta¬ 
tions,  in  a  local  fashion,  to  topological  and  traffic  demands,  spatial-TDMA  can  be  made  to 
operate  in  such  a  mobile  environment. 

Another  drawback  of  spatial-TDMA  concerns  the  nature  of  the  traffic  flowing  in  the 
net.  Since  arcs  are  enabled  for  specific  periods  of  time  during  the  cycle,  efficient  use  of  the 
channel  requires  that  messages  be  available  for  transmission  during  these  times.  To  elucidate 
the  traffic  characteristics  most  suitable  for  spatial-TDMA,  imagine  that  traffic  arrives  to  the  arcs 
in  the  network  in  a  steady  flow.  The  queue  sizes  at  each  arc  then  would  never  grow  large 
because  the  queue  could  be  steadily  decreased  during  the  transmit  periods  of  the  cycle  for  that 
arc,  and  we  would  expect  the  message  delay  to  be  small.  In  contrast  to  this,  if  we  imagine  the 
same  average  flow  arriving  to  the  arcs  in  large  bulks,  with  long  interarrival  idle  periods,  large 
queue  sizes  are  formed  which,  because  the  arc  only  has  fixed  time  from  the  cycle  in  which  to 
transmit,  would  take  a  long  time  to  empty.  This  would  then  result  in  large  message  delays.  We 
can  thus  conclude  that  spatial-TDMA  is  most  suited  for  steady  streams  of  traffic.  We  might 
note  here  that  random  channel  access  protocols  are  most  suited  for  the  case  of  bursty  traffic. 
As  an  example  of  a  network  where  spatial-TDMA  would  be  suitable,  suppose  that  terminals 
were  connected  to  nodes  of  a  spatial-TDMA  network.  If  there  was  a  sufficient  number  of 


CHAPTER  5 
RUDE-CSMA 


So  far  we  have  considered  the  interference  problem  under  the  conditions  of  using 
minimal  information  (in  Chapter  2)  and  of  using  global  information  (Chapters  3  and  4).  In 
this  section  we  consider  an  intermediate  approach  to  this  problem  that  uses  local  information 
which  nodes  obtain  when  they  sense  the  channel.  This  investigation  leads  us  to  define  a  two- 
parameter  family  of  protocols  for  the  multi-hop  environment  called  rude-CSMA. 


5.1  Introduction 

Let  us  first  review  the  CSMA  protocol  as  defined  for  the  single-hop  environment  that 
was  discussed  in  Chapter  1.  It  is  well  known  that  CSMA  [Klei75a]  is  an  efficient  channel  access 
protocol  for  single-hop  packet  radio  environments.  In  this  protocol,  nodes  of  the  network  sense 
the  channel  prior  to  transmitting  packets.  If  the  channel  is  sensed  busy,  the  sensing  node 
refrains  from  transmitting  (to  avoid  a  collision)  and  re-schedules  its  transmission  until  a  future 
time.  In  a  multi-hop  environment  sensing  the  channel  provides  information  only  about  a  sub¬ 
set  of  the  nodes  in  the  network  and  thus  a  transmission  by  one  node  that  is  surrounded  by 
silent  nodes  is  not  guaranteed  to  be  successful.  An  example  of  this  is  shown  in  figure  5.1.  In 
this  figure  suppose  node  1,  wishing  to  send  a  message  to  a  neighboring  node,  senses  the  chan¬ 
nel  idle  and  transmits  its  packet.  It  is  clear  in  this  network,  l's  packet  will  be  successfully 
received  only  if  node  4,  as  well  as  its  destination,  are  not  transmitting.  Thus  we  can  conclude, 
that  in  an  optimized  network  (where  the  objective  is  to  maximize  the  number  of  successful 
transmissions  in  the  network  over  time)  there  is  some  relationship  between  the  rates  at  which 
nodes  in  the  network  present  packets  to  an  idle  channel.  In  particular,  for  our  example,  the 
rate  that  node  1  presents  packets  to  an  idle  channel  in  an  optimized  network,  is  functionally 
dependent  upon  the  rate  that  node  4  presents  packets  to  the  channel. 

Another  relationship  can  be  hypothesized  by  supposing  that  node  1 ,  wishing  to  send  a 
packet,  senses  the  channel  busy.  In  single-hop  networks  this  implies  that  node  1  will  refrain 
from  transmitting  until  a  future  time.  In  our  example,  however,  if  node  1  senses  the  channel 
busy  due  to  a  transmission  from  node  3.  this  implies  (assuming  that  node  1  is  not  sensing  the 
channel  if  it  is  receiving  a  message)  that  node  3  is  transmitting  a  message  to  node  4.  Thus  if 
node  1  transmits  to  node  2,  its  reception  will  be  successfully  received  without  interrupting  node 


FIGURE  5.1 
A  sample  network. 

3’s  transmission.  Thus  it  appears  that  if  node  I  senses  the  channel  busy  it  should  not  neces¬ 
sarily  refrain  from  presenting  packets  to  the  channel.  Clearly  the  rate  at  which  it  offers  such 
packets  will  depend  upon  the  topology  and  traffic  generation  rates  of  the  nodes  in  the  network 
but  we  can  conclude  from  these  motivating  arguments  that,  in  some  cases,  throughput  would 
be  increased  if  nodes  transmitted  even  when  the  the  channel  was  sensed  busy.  This  charac¬ 
teristic,  namely  transmitting  even  though  other  nodes  are  using  the  channel,  explains  the  appel¬ 
lation  “rude-CSMA”  for  the  protocols  we  define  in  the  following  section.  We  should  mention 
that  these  relationships  were  also  discussed  in  Chapter  1  when  we  showed  in  Figure  1.1  that 
there  were  cases  when  it  was  advantageous  for  transmitters  r  \  and  r2  to  transmit  even  if  they 
sensed  a  busy  channel. 

This  hypothetical  reasoning  leads  us  to  question  how  the  rates  of  offered  packets  to  the 
channel  should  be  adjusted  to  achieve  maximal  throughput.  In  this  chapter  we  investigate  a 
family  of  protocols  which  attempt  to  use  the  local  information  contained  in  the  state  of  the 
channel  to  determine  the  rate  at  which  nodes  transmit  packets  on  the  channel.  The  question 
we  hope  to  shed  light  upon  in  this  chapter  is:  How  can  the  “free”  information  obtained  from 
sensing  the  channel  be  best  used  to  maximize  the  throughput  of  multi-hop  networks.  We  will 
see  that  the  family  of  protocols  we  create  to  answer  this  question  include  both  ALOHA  and 
CSMA  as  special  cases,  and,  being  generalizations  of  them,  should  have  an  optimized  perfor¬ 
mance  that  is  at  least  as  good  as  either  of  them.  To  further  motivate  the  protocols  that  will  be 
precisely  defined  in  the  next  section,  suppose  node  /  is  ready  to  transmit  a  packet  to  one  of  its 
neighbors.  After  sensing  the  channel,  suppose  /  finds  the  channel  busy.  Besides  the  binary 
information  contained  in  the  channel  status,  i  also  obtains  the  power  of  the  received  signal. 
Suppose  then  (Figure  5.2)  that  the  received  power  is  very  small.  This  implies  that  the  number 
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FIGURE  5.2 
Small  Received  Power 


FIGURE  5.3 
Large  Received  Power. 
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of  transmitters  within  range  of  /  is  small  and  furthermore  that  they  are  not  very  close  to  i.  If 
we  assume  that  all  nodes  transmit  with  the  same  power,  this  implies  that  the  immediate  area 
surrounding  i  is  free  of  transmitters  and  thus  that  there  a  possibility  that  if  i  transmits,  the  sig¬ 
nal  will  be  successfully  received  by  its  destination.  The  percentage  of  area  of  the  shaded  region 
of  the  circle  shown  in  the  figure  is  in  fact  the  probability  that  /  will  be  successful.  On  the  other 
hand,  suppose  the  received  power  is  large  (Figure  5.3).  This  implies  that  either  there  are  a  few 
transmitters  very  close  to  i  or  many  far  away.  In  the  first  case  the  environment  covered  by  a 
transmitter  very  close  to  /  would  greatly  overlap  that  covered  by  /  and  thus  a  transmission  by  / 
would  most  likely  cause  a  collision.  In  the  second  case.  Figure  5.3,  since  there  are  a  large 
number  of  transmitters  around  /,  the  probability  that  one  of  /’ s  neighbors  hears  a  silent  channel 
is  very  small  and  thus  if  /  transmits,  its  signal  would  have  a  high  probability  of  causing  a  colli¬ 
sion.  We  can  conclude  from  this  hypothetical  scenario  that  the  decision  to  transmit  depends 
upon  the  number  of  transmitters  and  their  locations  in  /’ s  immediate  environment.  In  the  pro¬ 
tocols  defined  in  this  chapter  we  will  assume  that  given  the  power  of  the  received  signal  on  the 
channel,  the  number  of  transmitting  nodes  within  the  hearing  distance  of  any  given  node  can 
be  estimated.  Later  we  will  show  that  this  assumption,  although  being  unrealistic,  is  not  too 
restrictive. 


5.2  The  Protocol 


We  will  define  the  state  S  of  the  network  to  be  a  binary  vector 
S  —  (sj(5),  s2(S),  •  •  •  ,  s„(S) )  where  n  is  the  number  of  nodes  in  the  network  and  s,(S)  «  1 
if  node  /  is  transmitting  and  0  otherwise.  For  a  given  node,  say  node  /,  there  exists  a  subset  of 
the  other  nodes  of  the  network  which  are  its  neighbors.  This  set  is  denoted  by 

A,  —  [j  |  node  j  is  within  hearing  distance  of  node  / } 

In  general  for  a  multi-hop  network.  A,  is  a  proper  subset  of  all  the  nodes  in  the  network.  Thus 
i's  decision  to  transmit  will  depend  only  upon  the  states  of  the  nodes  contained  in  its  hearing 
region  A,.  For  a  given  state  S,  suppose  node  /  has  a  packet  which  is  ready  to  be  transmitted. 
Node  i  senses  the  channel  and  based  upon  the  power  of  the  received  signal,  estimates  the 


number  of  transmitters  within  its  local  hearing  distance.  Denote  this  number  by: 


-Vi  (S)  -  £  s;(5) 

a  4 


Likewise  the  number  of  neighbors  not  transmitting  in  /' s  hearing  range  can  be  estimated  as: 

Nb(S)  -\A,  I  -  .VI  (S) 


We  assume  in  our  mathematical  model  that  these  estimates  are  exact.  According  to  these 
values  node  /  adjusts  the  rate  at  which  it  presents  packets  to  the  channel,  and  in  particular,  for 
state  S,  the  rate  at  which  node  i  transmits  (assuming  that  s,(S)  -  0)  is  given  by: 
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(5.1) 


where  x  and  y  are  given  parameters  of  the  protocol  and  yo  is  a  given  arrival  rate  of  packets  to 
each  node  of  the  network  (from  the  attached  host  computer).  In  the  case  that  node  /  is  a 
transmitter,  (s,(S)  —  1),  we  assume  the  packet  transmission  time  is  exponentially  distributed 
with  an  average  length  of  1/p,  time  units  and  thus  the  rate  at  which  i's  transmission  stops  is 
given  by: 


rj(S)«M  (5.2) 

For  any  given  state  S,  with  these  rate  definitions,  one  can  determine  the  rate  at  which  the  sys¬ 
tem  changes  into  any  other  state  S\  and  it  is  clear  that  the  state  vector  and  rate  definitions 
define  a  continuous  time,  finite  state,  Markov  process.  A  generalized  form  of  this  process  can 
be  found  in  [Kell79a].  Observe  that  this  two  parameter  family  of  protocols  contain  ALOHA 
(x-1,  y-1)  and  CSMA  Cx-l,  y«0)  as  special  cases.  To  optimize  the  performance  of  the  net¬ 
work  we  need  to  determine  the  ( x,y )  values  that  maximize  the  expected  number  of  successful 
transmissions  ovdr  a  unit  time  interval.  For  a  given  state  S,  define  E(S)  to  be  the  expected 
number  of  successful  transmissions  for  that  state.  It  is  clear  that  E(S)  can  be  calculated  if  the 
traffic  matrix,  topology,  and  capture  ratio  for  the  system  are  known.  For  example,  suppose  in 
Figure  5.4  that  nodes  are  equally  likely  to  transmit  to  any  of  their  neighbors.  Nodes  in  Figure 
5.4  are  labeled  1  if  they  are  transmitters  and  0  otherwise.  We  can  determine  E(S)  by  calculat¬ 
ing  the  probability  that  silent  nodes  in  the  network  successfully  receive  a  packet  destined 
towards  them.  Thus  node  2  in  the  figure  is  adjacent  to  only  one  transmitting  node  (node  1) 
which  transmits  in  that  direction  with  probability  1/4.  The  expected  success  then  for  node  2  is 
1/4  as  it  is  for  nodes  3  and  4.  In  the  same  manner  nodes  6  and  7  have  a  1/4  probability  of 
receiving  node  5’s  transmission  successfully,  whereas  node  8  has  no  chance  of  receiving  either 
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FIGURE  5.4 
Calculating  E. 

node  5’s  or  node  9's  reception  since  they  will  collide.  For  the  state  shown  then,  the  expected 
success  is: 

E(S)  —  3  (-7)  +2<t>  -  t 
4  4  4 

Similar  procedures  can  be  generated  that  handle  cases  with  different  assumptions  of  the  traffic 
matrix  and  capture  parameter  than  those  of  the  previous  example.  Observe  that  these  assump¬ 
tions  correspond  to  a  packet  radio  network  in  which  terminals  have  only  a  transmitter  buffer, 
and  radio  units  are  assumed  to  capture  the  first  transmitted  signal  they  receive.  This  perfect 
time  capture  is  an  ideal  assumption  and  could  alternatively  be  modeled  as  an  upper  bound  on 
throughput  to  a  packet  radio  network  without  time  capture  in  which  the  exponential  message 
length  corresponds  to  a  bulk  of  very  small  packets.  To  explain  this  suppose  a  radio  unit  i  is 
successfully  receiving  a  transmission  from  another  unit  j  over  the  time  period  [0,r].  It  is  possi¬ 
ble  that  another  node  k,  within  i's  hearing  radius  but  not  within  s  range,  starts  transmitting  at 
time  t.  Clearly,  assuming  non-capture,  this  implies  that  /  will  hear  a  collision  at  r,  and  will  not 
successfully  receive  the  rest  of  j's  message  if  j  was  sending  a  single  packet  message.  However 
if  we  assume  that  j  was  sending  a  series  of  very  short  packets,  then  only  those  sent  after  time  t 
will  be  lost.  In  the  limit  as  the  packet  size  goes  to  zero,  E(S)  then  will  measure  the  fraction  of 
time  the  channel  that  is  being  used  successfully  when  the  system  is  in  state  5. 


Let  n(S,x,^)  be  the  steady  state  probability  of  state  S  for  a  given  (xyy).  In  Appendix 
C  we  show  that  with  rate  equations  (5.1)  and  (5.2),  11(5, x,.y)  is  given  by: 

n/c  \  r  W(S>  -®0<sl  BAS) 

II(S,jr,^)  -  C  p  x  u  y  1 


where  A/(5)  -  s,(5) 

C  -  II(0,*,.y)  x~B°l0) 

P  “  Wm 

n(0,jc,>-)  -  [x_8°(0>  £  pw<5)  JC-flo<5,/,(*']-1 

In  this  equation  C  is  a  normalization  constant,  and  for  state  5,  M(S)  is  the  number  of 
transmitters,  B0(S)  is  the  number  of  adjacent  nodes  that  are  not  transmitting,  and  B{(S)  is  the 
number  of  adjacent  transmitters.  Using  this  we  can  write  the  expected  number  of  successes  in 
the  network  over  a  unit  time  interval  as: 

E(x,y)  n(S,x,.y)  E(S) 
s 

Naturally  we  would  like  to  maximize  this  function  over  feasible  (x,y)  values.  Besides  the 
non-negativity  of  x  and  y  however,  there  is  a  constraint  concerning  the  average  rate  at  which  a 
node  presents  packets  to  the  channel.  This  actual  rate  must  not  be  greater  than  the  arrival  rate 
of  packets  to  that  node,  which  is  defined  to  be  y0-  Rate  equation  (5.1)  shows  that  this  actual 
transmission  rate  is  a  function  of  5,  and  this,  averaged  over  all  states,  must  therefore  be  less 
than  -yo-  thus  have,  letting  s,(5)  be  the  complement  of  the  i'h  component  of  5: 

£  yo  x^(S)  y"'[ (5>  1,(5)  n(5,x,y)  ^  y0  /  *  1,2,  •  •  •  ,n 
s 

For  a  given  topology  and  traffic  matrix,  we  can  thus  formalize  the  mathematical  program  P  as: 


Program  P 

max  E  4  Y  11(5, x,y)  E(S) 

<•'•>’>  s 

subject  to: 

_  .vj(S)  .Vj<s>  _ _,c  ,  ^  . 

i  -  £  x  y  $(5)  n(s,x,y)  >  o 

s 

x  >0 


1,2, 


•  ,« 


y  >0 


This  is  the  problem  that  we  will  address  in  this  chapter.  We  should  note  however  that  optimal 


(x,y)  values  are  functions  of  the  topology  and  traffic  matrix  of  the  network.  We  will  call  the 
first  constraint  equation  in  P  the  flow  constraint  and  will  say  that  it  is  saturated  if  the  equation  is 
an  equality  for  some  equation. 


5.3  Discussion  of  Results 

In  this  section  we  will  discuss  the  results  we  obtained  when  P  is  optimized  over 
different  topologies.  We  wrote  computer  programs  which,  for  a  given  topology,  calculated 
E(S)  and  n(S,x,y)  and  then  optimized  P  over  all  feasible  (x,y)  pairs.  Besides  special  topolo¬ 
gies  which  we  studied,  we  also  ran  P  for  connected  networks  which  were  randomly  generated 
for  varying  mean  densities.  Unless  stated  otherwise,  in  all  networks  we  have  assumed  a  uni¬ 
form  local  traffic  matrix.  As  typical  of  most  Markovian  models,  the  state  space  of  the  system 
grows  exponentially  with  the  number  of  nodes  of  the  network,  and  thus  to  keep  the  program 
computationally  tractable,  we  restricted  our  optimization  to  relatively  small  networks.  We  do 
not  believe  however,  the  applicability  of  our  results  are  dependent  upon  the  size  of  the  net¬ 
work. 


5.3.1  Special  Topologies 

In  figure  5.6  we  plot  the  values  of  x  and  y  that  achieve  optimal  performance  for  the 
network  of  figure  5.5,  which  is  a  generalized  version  of  the  network  used  to  motivate  this  work. 
The  curves  for  the  graph  in  figure  5.1  show  similar  behavior.  These  curves  have  many  interest¬ 
ing  properties.  First  we  observe  that  over  the  range  of  p  shown  i.<  *e  are  three  distinct  types  of 
behavior.  For  very  small  p  values,  p<.15,  we  see  that  x  and  E  increase  rapidly  while  y 
remains  very  small.  The  increase  of  £  over  this  range  is  explained  by  the  fact  that  for  small  p 
values,  there  are  very  few  transmissions,  and  thus  very  few  which  cause  collisions  (or  inter¬ 
preted  in  the  case  of  perfect  time  capture,  few  transmissions  that  aren't  first)  and  thus  increas¬ 
ing  p  tends  to  increase  £  in  a  linear  manner.  We  see  this  behavior  clearly  as  p  goes  from  0.05 
to  0.1  during  which  £  doubles  from  about  0.35  to  about  0.6.  As  p  continues  to  grow  beyond 
0.1  the  increase  is  less  than  linear  due  to  some  collisions  in  the  network.  It  is  clear  that  over 
this  range,  since  there  are  so  few  transmissions  in  the  network,  the  y  parameter  of  rate  equa¬ 
tion  15.1)  is  non-critical.  Using  an  interactive  program  we  wrote  to  determine  how  £  varied  as 
a  function  of  y,  we  found  very  little  change  over  all  values  of  y,  thus  the  sudden  increase  in  y 
at  ps. 2  should  not  be  interpreted  as  demonstrating  singular  behavior. 


FIGURE  5.5 
A  Lattice  Network. 


The  increase  in  x  over  the  range  0.05  <0.2  can  be  explained  by  first  observing  that 

for  this  range  the  chance  of  collisions  is  small  and  thus  nodes  present  packets  to  the  channel  as 
fast  as  they  can  generate  them,  which  in  our  case  is  approximately  equal  to  yr>x2.  We  should 
mention  that  for  these  small  values  of  p,  the  optimal  x  value  lies  on  the  boundary  of  the  feasi¬ 
ble  region,  hence  the  flow  constraint  of  P  is  saturated.  We  can  explain  the  linear  behavior  of 
the  x  curve  over  this  region  with  an  intuitive  argument.  From  our  previous  discussion  we 
know  that  £  is  linear  in  p  for  small  values,  and  it  can  be  approximated  closely  by  E -4.6p  +  .1. 
For  the  most  part,  successful  states  for  these  small  p  values  correspond  to  one  transmitter  in 
the  network.  Thus  using  our  equation  for  n(S,x,.y),  we  can  write  E  -  kpx~2  33,  where  k  is  a 
constant  of  proportionality.  We  can  then  write: 

4.6p  +  .1  -  kpx~233  (5.3) 

Using  the  initial  values  from  figure  5.6  we  can  determine  A:  —  10.09.  Solving  equation  (5.3)  for 
x  therefore  yields: 

|  10,09  p  ) 1/2  33 
l  4.6  p  +  .1 J 

which  for  small  values  of  p  is  very  closely  approximated  by: 

x  —  2p  +  1.1 

This  then  explains  the  linear  shape  and  slope  of  x  over  this  region.  In  summary  then,  over  the 
first  region  of  the  curve  in  figure  5.6,  £  and  x  are  approximately  linear  with  slopes  4.6  and  2 
respectively,  and  y  is  not  a  critical  parameter  of  the  system. 


The  second  region  of  the  graph,  from  p- 0.2  to  p-0.35,  will  now  be  explained.  First 
we  observe  that  the  expected  success  rises  only  slightly  over  this  range.  Although  there  are 
more  transmissions  in  this  range,  there  are  also  more  collisions  which  limit  the  number  that  are 
successful.  The  sudden  increase  of  y  at  p=0.2,  as  mentioned  previously,  should  not  be  looked 
on  as  being  as  singular,  but  does  demonstrate  the  increased  importance  of  y's  effect  on  the 
throughput  of  the  system.  Observe  that  over  this  range  y  <  1,  and  thus  from  equation  (5.1), 
this  parameter  acts  to  inhibit  transmissions  when  there  is  a  neighboring  transmitter  (although 
not  to  the  extent  of  excluding  such  transmissions).  This  supports  our  conclusion  from  the 
scenario  in  the  beginning  that  sometimes  in  a  multi-hop  network  it  is  beneficial  to  transmit 
even  if  the  channel  is  sensed  busy.  The  decrease  of  x  over  this  range  arises  because  it  is  no 
longer  on  the  boundary  of  the  feasible  region.  However,  since  x  is  greater  than  1  in  this 
region,  we  conclude  that  this  parameter  tends  to  increase  the  rate  of  packets  offered  to  the 
channel  during  idle  channel  periods. 

Over  the  region  p>0.35,  E  is  again  not  altered  to  any  great  measure.  Both  x  and  y 
decrease  now  in  an  effort  to  prevent  collisions  on  the  channel.  As  p  becomes  much  larger,  we 
see  that  x<l,  which  indicates  that  even  the  x  parameter  tends  to  act  to  inhibit  the  number  of 
packets  on  the  channel.  Over  this  region,  the  flow  constraint  equation  is  not  saturated  for  any 
node  in  the  network  indicating  that  the  effective  transmission  rate  of  packets  on  the  channel  is 
less  than  y0-  Hence  optimizing  the  throughput  of  the  network  has  the  effect  of  limiting  the 
flow  of  packets  in  the  network.  We  will  explain  the  general  shape  of  the  x  curve  when  discuss¬ 
ing  random  topologies  in  the  next  section. 

The  next  special  topology  we  will  study  is  that  of  ring  networks  as  shown  in  figures  5.7 
and  5.10.  In  each  of  these  rings  we  have  assumed  that  messages  always  travel  clockwise  (thus 
the  arrows  in  the  figures).  We  observe  that  y— 0  for  all  p  values  and  thus  nodes  in  such  net¬ 
works  should  refrain  from  transmitting  if  they  sense  a  busy  channel.  This  is  actually  obvious 
from  the  defined  topologies  since  the  sensed  transmission  would  either  be  directed  towards  the 
sensing  node  or  be  the  intended  recipient  of  the  sensing  nodes  transmission.  In  each  case  then 
transmitting  would  cause  a  certain  collision.  Knowing  that  the  nodes  will  use  the  CSMA  proto¬ 
col  allows  us  to  solve  for  the  steady  state  throughput  using  an  alternative  method,  thus  giving 
intuitive  insight  into  the  solution  and  corroborating  our  equation  for  n(5,.v,y).  To  generate 
the  states  for  this  process,  we  observe  there  are  five  states  which  have  one  transmitter,  and  five 
states  that  have  two  transmitters  (we  do  not  allow  two  transmitters  to  be  adjacent).  We  can 
thus  form  the  Markov  chain  shown  in  figure  5.9.  In  this  figure  a  is  the  arrival  rate  of  messages 


FIGURE  5.7 
A  Ring  of  Five  nodes. 

into  each  node  and  corresponds  to  yox2  in  our  model.  Solving  this  chain  yields: 

Po  ™  1/(1  +5(a/(i)  +5(a/fi)2) 

Pi  -  S(a//x)Po 

P2  -  5(a/p)2/>0 

where  P,  is  the  equilibrium  state  probability  of  state  /'.  It  can  easily  be  checked  that  using  our 
equation  for  n(S,x,y)  yields  the  same  set  of  equations  with  a-yo*2.  In  each  of  states  1  and  2 
there  is  an  expected  success  of  1  packet,  and  thus  maximizing  the  throughput  of  the  network 
corresponds  to  maximizing  the  probability  that  the  system  is  in  states  1  or  2.  This  can  most 
easily  be  performed  by  minimizing  the  probability  the  system  is  in  state  0.  It  is  clear  that  this 
occurs  for  (a/p.  )-*<»,  and  in  figure  5.8  we  see  the  monotonic  rise  of  throughput  with  p,  and 
the  increase  in  x  which  corresponds  to  lying  on  the  boundary  of  the  feasibility  region. 

Quite  a  different  phenomena  is  observed  for  the  case  of  a  ring  with  an  even  number  of 
nodes  as  shown  in  figure  5.10.  To  intuitively  motivate  this  we  draw  the  corresponding  Markov 


chain  in  figure  5.1 1.  which  can  be  solved  to  yield: 

Pet  —  l/(  1  +2(a//u) '  +9(«//i)’  +6(a//x)  ) 

/>,  -  6Wm)/>0 

P,  -  6(a/»)2Po 

P:  -  UaU)2Po 

Pi  -  2(a/M)3/>o 

where  a—yojr  and  state  2  corresponds  to  a  configuration  congruent  to  1  and  3  being 
transmitters  and  state  2'  to  1  and  4  being  transmitters. 


FIGURE  5.12 

Curves  for  the  Six  Node  Ring  (y— 0). 


The  expected  success  for  each  state  is:  0  for  states  0  and  3.  and  1  for  state  2.  and  2  for  state  2'. 
Thus  we  wish  to  maximize  P2  +  2P?  over  a.  This  situation  is  much  different  than  the  previous 
odd  number  ring  because  increasing  a  has  the  effect  of  increasing  Pj  which  has  an  expected 
throughput  of  0.  The  objective  then  is: 


Max 


_ 6 (g/p)  -r  12<«<p)~ _ 

1  ■+■  2 ( <xi p )  ^  +  9(  g/p  ) "  6( a '  p  ) 

This  maximum  occurs  for  (u/pl”I  at  which  point  the  objective  function  is  equal  to  1.  In 
figure  5.12  we  see  this  behavior  in  the  curves  for  x  and  E.  For  p<.3  the  flow  constraint  in  P 
is  satisfied  and  the  («/'p)  values  are  much  less  than  the  optimal  value  of  1  because  there  is  not 
enough  input  to  satisfy  this  condition.  For  p>0.3  however,  the  input  into  the  system  is 
sufficient  to  reach  optimal  throughput  and  the  x  values  shown  act  to  control  the  rate  of  this 
input  to  achieve  optimal  throughput.  For  example  for  p= 0.5,  v=1.4  yields  «-=p.v’=0.98  (where 
we  have  normalized  p  to  be  equal  to  1).  We  notice  over  this  range  that  the  optimal  throughput 
of  1  packet  per  time  unit  is  obtained. 

We  observe  with  rings  a  striking  difference  in  the  behavior  of  the  system  (in  terms  of 
the  a— ■yo-’f"  parameter)  which  depends  on  the  number  of  nodes  in  the  network.  In  particular 
for  an  odd  number  of  nodes,  optimal  performance  was  obtained  for  a— *°°.  whereas  for  an  even 
number  of  nodes  we  saw  that  a<<»  achieved  maximal  throughput.  This  dependency  results 
from  combinatorial  considerations  on  how  transmitters  could  be  distributed  in  the  nodes  such 
that  no  two  adjacent  nodes  were  transmitters.  In  odd  node  rings,  when  the  maximum  number 
of  transmitters  is  generated,  the  expected  throughput  for  the  configuration  is  greater  than  zero 
(in  fact  it  is  always  equal  to  1),  whereas  for  even  node  rings,  this  configuration  yields  a  zero 
throughput.  This  tends  to  make  optimal  «  values  for  odd  node  rings  greater,  but  not  neces¬ 
sarily  infinite,  than  those  for  even  node  rings.  We  might  question  if  such  behavior  is  obtained 
when  the  rings  are  cut  at  one  spot  and  made  into  tandems. 


FIGURE  5.13 

A  Five  Node  Tandem  Network. 

In  figures  5.13  and  5.15  we  show  a  five  and  six  node  tandem  with  their  associated  curves  in 
figures  5.14  and  5.16.  Both  graphs  show  similar  behavior.  Once  again  for  p<0.4  the  flow  con¬ 
straint  equation  in  P  is  saturated  and  after  that  the  x  parameter  acts  to  meter  the  actual  flow  of 
packets  into  the  network. 
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FIGURE  5.16 

Curves  for  the  Six  Node  Tandem  (y-0). 

5.3.2  Random  Topologies 

In  this  section  we  report  on  the  results  which  were  obtained  when  P  was  run  on  graphs 
which  were  randomly  generated.  Many  such  graphs  were  created,  with  varying  mean  densities 
and  all  such  graphs  exhibited  similar  behavior.  This  characteristic  allows  us  to  make  the 
blanket  statement  that  in  a  random  multi-hop  network,  over  the  continuum  of  protocols 
defined  by  all  possible  (.x.y)  pairs,  optimum  performance  was  obtained  for  those  pairs  that  had 
y— 0.  As  stated  in  the  introduction  this  corresponds  to  a  CSMA  type  protocol.  From  the  many 
graphs  which  were  analyzed  we  have  picked  two  to  show  here.  This  selection  was  based  more 
on  the  aesthetic  qualities  of  the  graphs  rather  than  any  particular  characteristic  of  the  curves 
associated  with  solving  P.  In  figure  5.17  we  show  a  7  node  network  with  a  mean  number  of 
neighbors.  <V  equal  to  2.28.  The  corresponding  set  of  curves  are  shown  in  figure  5.18.  We  see 
that  for  low  values  of  p,  both  E  and  .x  are  linear  indicating  that  there  are  very  few  collisions  in 
the  network.  The  slope  of  E  becomes  nearly  zero  for  p>0.2  implying  that  after  this  point  the 
amount  of  new  successful  transmissions  in  the  network  is  small  in  comparison  to  the  number  of 


N  -  2.285 


FIGURE  5.17 

A  random  graph  with  .V  =  2.285  . 


collisions  that  arise  from  the  increased  traffic  load.  The  shape  of  the  x  curve  for  p>0 .3  can  be 
explained  by  first  noting  that  over  this  range  the  flow  constraint  of  P  is  not  saturated.  Next  we 
make  a  definition.  Suppose  for  each  node  /,  we  investigate  all  state  vectors  having  the  property 
that  i  and  all  its  neighbors  A  are  silent.  Under  these  conditions,  if  /  senses  the  channel,  it  will 
detect  that  it  is  idle  and  offer  packets  to  the  channel  at  a  rate  defined  by  the  protocol.  Let  /.  be 
the  set  of  all  state  vectors  satisfying  this  condition.  We  can  then  define  the  effective  number  of 
idle  neighbors,  denoted  by  N' ,  by  finding  the  average  number  of  silent  neighbors  a  randomly 
selected  idle  node  has  when  it  senses  the  channel  prior  to  a  transmission.  This  can  be  written 
as: 

V*  \A, I  £  n(S,.x,y) 

,  sa, 

Because  of  topological  variations,  some  nodes  are  more  likely  to  be  both  idle  and  surrounded 
by  idle  neighbors  than  others  (in  fact  the  probability  is  decreasing  with  the  number  of  neigh¬ 
bors)  and  thus  in  general  j VViV.  Define  /  to  be  the  average  time  a  random  node  in  the  net¬ 
work  waits  before  transmitting  when  it  hears  an  idle  channel.  Since  E  does  not  varv  much  over 


FIGURE  5.18 

Curves  for  graph  in  figure  5.17  (y^O). 


p>0.3  it  is  not  surprising  to  find  that  .v  varies  as  a  function  of  p  to  preserve  the  fraction  of 
time  nodes  transmit  when  sensing  an  idle  channel.  This  has  been  verified  by  observing  (bv 
using  an  interactive  program)  that  the  probability  distribution  for  state  vectors  does  alter 
significantly  as  p  ranges  greater  than  0.3  for  optimized  x  values.  We  thus  can  equate  rates  to 
obtain: 

pxs  -  1  // 

We  can  solve  for  /  and  N'  by  picking  two  points  from  figure  5.18.  At  p-0.72  we  have  v=l 


and  thus  can  write  0.72—1/  /'  or  /'-l  .36.  At  p=0.32  we  have  .32 <  18) v  -.72  which  implies  that 
V— 1.4.  Thus  we  can  write  v  as  a  function  of  p  as: 


v  -  (l  '/p)1  v’  -  <-72/p) p>0.3  (5.4) 

Equation  (5.4)  was  used  to  generate  the  second  curve  in  figure  5.18  and  we  see  a  relatively 
close  match.  In  general  then,  for  a  CSMA  environment,  we  can  write  a  general  equation  relat¬ 
ing  v  with  p  as: 

.v  =  <!//>>'  v' 

Intuitively,  for  the  graph  of  figure  5.17.  we  would  expect  .V*<  V  since  nodes  having  the  greater 
number  of  neighbors,  node  4  for  instance,  have  a  much  smaller  probability  of  being  both  idle 
and  surrounded  by  idle  neighbors,  than  a  node,  say  node  1 .  that  has  far  fewer  neighbors. 


7 


FIGURE  5.19 

A  random  graph  with  .V— 3.14. 


In  the  second  graph  of  this  section,  figure  5.19.  we  again  see  in  its  associated  set  of 
curves,  figure  5.20,  similar  behavior.  As  stated  in  the  beginning  of  this  section,  this  behavior 
was  typical  of  all  the  curves  for  the  random  graphs  we  generated. 
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FIGURE  5.20 

Curves  for  figure  5.19  h  «0>. 
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5.4  Conclusions 


In  this  chapter  we  have  defined  rude-CSMA  which  is  a  two  parameter  continuum  of 
protocols  for  multi-hop  channel  access.  This  family  of  protocols  has  been  shown  to  include 
both  ALOHA  and  CSMA  as  special  cases.  The  optimum  performance  of  rude-CSMA  for  spe¬ 
cial  topologies  was  investigated,  and  we  have  seen  that  for  lattice  type  networks,  under  certain 
conditions,  transmitting  even  after  sensing  a  busy  channel  tends  to  increase  the  throughput  of 
the  system.  In  all  other  topologies  studied  including  rings,  tandems,  and  random  graphs, 
optimal  performance  was  obtained  when  the  protocol  used  was  CSMA  with  optimized  channel 
input  rates.  Thus  we  must  come  to  the  interesting  conclusion  that  for  practical  networks. 
rude-CSMA  is  not  that  “rude"  afterall.  This  justifies  the  statement  in  the  introduction  that 
assuming  nodes  could  determine  the  number  of  transmitters  from  the  received  power  of  sensed 
signals  was  not  too  restrictive  an  assumption.  The  binary  information  contained  obtained  from 
sensing  the  channel  appears  to  be  sufficient  to  determine  the  optimal  transmission  policy.  The 
curves  generated  for  these  protocols  under  varying  assumptions  have  proved  to  be  intuitively 
insightful  and  their  interesting  properties  have  been  explained. 


CHAPTER  6 
CONCLUSIONS 


In  the  previous  chapters  we  analyzed  the  interference  problem  under  various  sets  of 
policy  state  information.  We  have  also  seen  that  the  difficulty  of  the  problem  concerned  the 
fact  that  all  nodes  of  a  packet  radio  network  are  interrelated,  and  that  changing  the  behavior  of 
one  node  in  such  a  network,  causes  global  effects.  Coupled,  dependent  systems  like  these  are 
very  difficult  to  analyze  and  our  approach  in  each  of  the  models  we  analyzed,  undertook  in 
some  manner  a  means  to  control  the  dependency.  In  Chapter  2  we  essentially  ignored  the 
dependency  and  determined  the  best  performance  obtainable  for  slotted-ALOHA  networks  in 
which  only  the  mean  density  of  terminals  on  the  plane  was  known.  This  model  corresponded 
to  solving  the  interference  problem  using  a  minimal  set  of  policy  state  information,  and  in 
Chapter  3  we  established  an  upper  bound  on  the  performance  for  all  protocols  which  use  such 
information.  In  Chapter  4  we  strictly  controlled  the  dependency  between  nodes  of  the  network 
by  using  the  global  information  of  their  location  and  traffic  characteristics.  We  saw  the  strong 
inter-relationship  between  all  of  the  nodes  when  we  attempted  to  determine  optimal  time 
frames  for  networks  using  spatial-TDMA,  and  needed  to  resort  to  solving  the  capacity  assign¬ 
ment  problem  over  the  class  of  frames  that  were  randomly  generated.  The  coupling  between 
the  nodes  of  the  network  was  again  encountered  in  the  fifth  chapter  where  we  sought  to  control 
the  inter-dependencies  by  determining  a  protocol,  rude-CSMA.  which  used  only  local  informa¬ 
tion  to  determine  when  to  transmit  packets  on  the  channel.  The  coupling  in  this  model  was 
restricted,  for  each  node,  to  the  subset  of  nodes  in  its  local  environment.  The  particular  rate 
definitions  given  for  the  protocol  allowed  such  a  de-coupling  and  also  allowed  us  to  solve  for 
the  steady  state  probabilities  of  the  system  state  in  closed  form. 

We  may  pose  the  following  question:  Where  does  the  subject  go  from  here?  It  is  the 
author’s  opinion  that  probably  the  most  important  area  of  research,  that  would  shed  light  on 
the  interference  problem  as  well  as  numerous  other  such  such  problems,  concerns  the  develop¬ 
ment  of  more  sophisticated  mathematical  tools  which  are  especially  designed  to  work  with 
dependent  systems  of  interacting  processes.  Most  of  the  relevant  mathematics  used  to  analyze 
such  systems,  such  as  stochastic  processes  and  probability  theory,  are  only  tractable  for  systems 
in  which  the  random  variables  are  independent.  Short  of  having  any  other  mathematical  tools, 
the  network  analyst  often  makes  the  rather  blatant  assumption  that  the  random  variables  associ¬ 
ated  with  the  network  are  independent.  This  assumption  is  certainly  false.  Although  the 
approximations  that  arise  from  such  assumptions  often  lead  to  models  that  ire  good  predictors 
of  network  performance,  without  proper  validation  procedures  there  seems  tc  be  no  mathemati¬ 
cal  way  to  determine  the  effect  that  certain  assumptions  have  on  the  accuracy  of  the  resultant 


network  model.  This  puts  network  designers  in  the  rather  nebulous  position  of  not  knowing 
what  to  believe,  and  forces  them  to  simulate,  or  in  some  other  way.  to  validate  the  model  and 
thus  nullify  the  "precision  and  cost  effectiveness”  of  analysis.  A  current  trend  in  solving  cou¬ 
pled  systems  of  this  type  consists  in  establishing  sets  of  equations  for  the  state  of  the  system  in 
terms  of  unknown  parameters.  These  parameters  form  the  coupling  between  state  variables. 
The  solution  of  these  equations  then  involves  guessing  initial  values  for  these  parameters  and 
then  iterating  through  the  equations  until  they  converge  to  a  final  set  of  values.  These  final 
parametric  values  form  a  consistent  set  of  equations  for  the  system  but  do  not  necessarily 
correspond  to  an  optimal  solution  of  the  system.  Other  methods  for  bypassing  the  coupling  of 
the  random  variables,  such  as  assuming  average  behavior  in  all  relationships,  suffer  from  the 
same  flaws,  namely,  optimality  is  not  easily  shown  and  procedures  are  needed  to  validate  the 
model.  There  is  thus  a  great  need  for  more  sophisticated  mathematical  tools. 

In  pertaining  to  extensions  and  generalizations  of  this  dissertation,  we  can  make  the  fol¬ 
lowing  comments.  In  Chapter  2  we  assumed  that  terminals  were  uniformly  distributed  on  the 
plane  in  an  effort  to  model  a  mobile  packet  radio  network.  In  many  networks,  however,  there 
is  a  tendency  for  nodes  to  form  clusters  rather  than  to  be  uniformly  spread.  For  example,  taxis 
tend  to  be  more  densely  populated  in  areas  of  a  city  in  which  the  demand  for  transportation  is 
greatest.  To  model  such  a  network,  one  would  have  to  change  the  model  of  the  second  chapter 
to  include  the  case  where  terminals  were  distributed  on  the  plane  according  to  a  non-uniform 
density.  This  would  imply  that  the  range  and  probability  of  transmission  would  depend  upon 
the  local  topology  of  the  network.  In  Chapter  4  we  solved  the  capacity  assignment  problem  for 
spatial-TDMA  networks  over  the  class  of  frames  that  were  randomly  selected  and  claimed  that 
this  formed  an  upper  bound  to  the  mean  delay  obtained  by  an  optimal  frame.  An  important 
contribution  would  be  an  algorithm  that  efficiently  produced  an  optimal  assignment  of  time 
slots  to  the  cliques  of  the  network.  The  author  of  this  dissertation  worked  intensely  on  this 
problem  but  was  only  successful  in  developing  heuristic  solutions  that  required  substantial  com¬ 
putation  for  their  execution.  Since  spatial-TDMA  could  be  a  very  useful  protocol  for  local  area 
radio  networks,  optimal  performance  is  an  important  consideration.  We  commented  in  the  con¬ 
clusions  to  Chapter  4  that  a  form  of  spatial-TDMA  that  worked  in  a  mobile  environment  would 
be  desirable.  Although  the  author  has  developed  a  formalization  of  such  a  protocol,  more  work 
needs  to  be  done  to  make  the  protocol  both  efficient  and  practical.  We  can  make  similar  com¬ 
ments  about  possible  generalizations  to  rude-CSMA  of  the  fifth  chapter.  In  that  model,  solu¬ 
tions  were  obtained  only  for  the  case  where  the  average  rate  of  incoming  packets  to  nodes  of 
the  network  was  the  same.  Once  again  this  uniform  assumption  is  often  not  found  in  practice 
(in  the  above  example,  taxis  in  the  central  part  of  the  city  may  be  using  their  transmitters  more 
often  than  those  on  the  outskirts).  Unfortunately  assuming  different  rates  changes  the  tracta- 
bility  of  finding  a  closed  form  solution  to  the  Markov  process  defined  by  the  state  transition 
rates,  and  thus  increases  the  complexity  of  finding  optimal  parametric  values.  If  an  efficient 
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APPENDIX  A 

An  Expansion  of  the  Integral  of  Chapter  2 
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We  then  wish  to  find  F(l).  To  do  this  we  first  set  up  a  recurrence  relationship  between  Fin) 
and  /■(/!+ 1).  Integrating  by  parts  we  find 
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Using  this  relationship,  we  claim 
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where  fly  —  1-3-5 . (2./+1)  -  (2y+l)!!.  We  prove  the  assertion  by  induction.  It  is  true 

for  /  -  2  from  the  recurrence  relation  (A.2).  Therefore  assume  it  is  true  for  /  «  L,  then  we 


can  show  it  is  true  for  /  -  Z.+1  as  follows:  By  assumption 
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Using  (A. 2),  the  last  term  can  be  written  as: 
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We  can  put  the  first  term  of  (A.4)  into  the  sum  of  (A.3)  to  obtain 
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and  thus  we  prove  that  the  expression  is  true  for  the  Z.+1'st  term.  If  we  continue  in  this 
manner  we  finally  obtain 
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Simplifying  a,  we  have 
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and  using  this  in  (A.5)  we  finally  obtain  (A.l). 


APPENDIX  B 

Calculating  the  densities  for  one.  two.  and  three  hops 


We  will  first  review  the  assumptions  and  definitions  discussed  in  Chapter  3.  Let  points 
be  distributed  according  to  a  Poisson  process  on  the  plane.  Create  a  non-directed  graph  by  join¬ 
ing  points  with  an  edge  if  the  euclidean  distance  between  them  is  less  than  or  equal  to  one  unit. 
For  a  randomly  selected  point.  po,  another  point  p  is  /  hops  away,  if  there  exists  a  path  from  po 
to  p  containing  /  —  I  other  points  and  no  other  path  exists  with  fewer  points.  Assuming  there 
exists  a  path  between  any  two  points  of  the  graph,  we  derive  in  this  appendix,  the  probability 
density  functions  for  the  euclidean  distance  between  points  that  are  separated  by  one,  two,  and 
three  hops.  In  particular  looking  at  figure  B.I,  we  will  derive  the  density  functions  for  the  vari¬ 
ables  ri,  ri,  /-j,  and  also  9. 


FIGURE  B.I 

A  typical  example  of  a  three  hop  configuration. 

We  will  use  the  following  notations  throughout  the  derivations  (these  are  shown  in  figure  B.I  >: 
p,  will  be  the  i'h  point  from  po,  r.  will  be  the  distance  from  p,  and  po.  r  will  be  the  distance 
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from  p\  and  p2.  «  is  the  angle  between  the  line  segments  (/>o.  p\)  and  ( P\.p->),  and  C  is  the  cir¬ 
cle  associated  with  p,.  We  will  assume  that  each  circle  has  a  radius  of  1  unit.  The  random  vari¬ 
ables  associated  with  the  above  quantities  will  be  denoted  as  R  .  <-),  and  Z.  The  derivations 
that  are  to  follow  require  extensive  use  of  the  Law  of  Cosines  and  the  Law  of  Sines  for 
numerous  defined  triangles.  To  facilitate  reading  these  derivations,  we  have  followed  the  for¬ 
mat  of  writing  the  initial  equation  that  follows  from  these  laws  on  the  left  side  of  the  page  and 
then  showing  the  equation  for  the  unknown  quantity  along  the  right  side.  It  was  felt  in  this 
way  the  triangle  that  was  used  in  the  derivation  could  more  easily  be  identified  by  the  reader. 
We  will  first  derive  the  densities  for  R\,  R 2,  and  ©  producing  conditional  densities  that  will  be 
necessary  for  the  much  harder  derivation  of  the  density  for  R  3. 


B.l  The  Density  for  R  { 

Since  points  are  randomly  distributed  according  to  a  Poisson  distribution,  and  since  we 
are  assuming  that  the  resultant  graph  is  connected,  we  see  in  figure  B.2  that  the  probability  dis¬ 
tribution  for  the  first  hop  can  be  formed  by  the  ratio  of  the  two  areas: 

P[R\  <  ril  =  irr^/TT  =  rf 

Differentiating  this  with  respect  to  r\  will  give  the  probability  density  function: 

/r  1  ( r  1 )  ■>  2  r\  (B.l) 


B.2  The  Density  for  R2 

In  figure  B.3  we  see  that  the  conditional  probability  distribution  can  be  written  as: 
P[R2<  r2l/?,  =  ri]  -  B(rur,)fA(rxA) 

where  B(rx,  r2)  is  the  shaded  area  in  figure  B.3  and  A  ( r | .  r2)  is  the  area  contained  in  C|  that  is 
not  contained  in  Co  or  B(r\,  r2).  It  is  clear  that  B(r\ .  r2 )  =  .4  (r(,  l)  -  A  (/•■,.  r2).  and  thus  we 
need  to  derive  an  expression  for  A  (•).  We  will  derive  this  function  under  more  general  condi¬ 
tions  where  the  radii  of  the  circles  are  not  both  equal  to  1.  We  will  use  figure  B.4  to  facilitate 
this  derivation. 
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FIGURE  B.2 

Calculating  the  density  for  the  first  hop. 


In  this  figure,  T(rur2)  is  the  area  formed  by  the  triangle  of  vertices  pn.  p t,  and  P. 
W(r |,  r2)  is  the  area  which  added  to  T(r\,  r2)  would  form  the  sector  of  radius  r2  and  angle  6. 
and  A(rhr2)/ 2  is  the  shaded  area  in  C\.  Since  A(rur2)/2  is  contained  in  a  sector  of  radius  1 
and  angle  »/>,  we  can  write: 

A(r,.r,)/2  -  <1,12  -  W(run) 


but  (V(r  |,  r2)  can  easily  be  seen  to  be: 

W(rur2)  =  &  ri'/2  -  r(r,.r:) 
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We  can  thus  derive  .-J(/-|.a2)  if  we  can  derive  equations  for  <U.  6,  and  T(r,.r2).  By  following 
the  format  discussed  in  the  introductory  remarks,  we  can  use  the  Law  of  Cosines  repeatedly  to 
obtain: 


1  =  r:-  +  r]1—  2r\r2cos(6)  —  6  =  cos  1 


a22  +  a,:  —  1 

2  r |  /-l 


a22»  1  |  -  —  2riCOS<-— >b)  — •  ii/  =  cos* 


-i I  r::~  1  ~  ri: 
2  a, 


We  can  thus  write  T(r\,r2)  in  terms  of  known  quantities  to  get  T(ru  r2)  =  /•!r2sin(f6)/2.  Com¬ 
bining  this  all  together  we  have: 
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We  can  thus  write: 


P[R2^r2\  /?,=a,]  =  1  - 
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The  conditional  density  h2(r2 1  /-t)  can  thus  be  found  as: 


.  ,  —  A'(r\,ri) 

/;2(a2|a,)  =  — — - —  Kr^l+r, 
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where  after  tedious  calculation  we  have: 


,  r2>V}  ,  a,K5  f  ,  ,  a,L4I'2  , 
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where  we  have  defined: 
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Using  this  conditional  density  and  the  density  for  R\  just  found,  we  can  find  the  unconditional 
density  for  R  2  as: 

i 

h2(r2)  =  J  h2(r2\ r\)  2r\dr\  (B.4) 


B.3  The  Density  for  0 

We  will  derive  the  density  for  0  for  «  <  tt  since  the  case  for  tf  >  tt  is  completely  analo¬ 
gous.  For  this  derivation  we  will  make  use  of  figure  B.5. 


IT  29  m 


FIGURE  B.5 

Calculating  the  density  for  0. 

In  this  figure.  r,)  is  the  shaded  area.  D  is  the  triangle  formed  by  points  />,.  and  P  and 
angle  a.  and  area  C  plus  D  make  up  the  sector  of  unit  radius  and  angle  Aiso  in  the 

figure.  F  is  the  triangle  of  vertices  pn.  p\.  and  angle  and  E  is  the  triangle  made  up  of  ver- 


The  conditional  density,  kift !  r , )  then  is: 

,  dB  ( ft  |  /■  | ) 

— 7 - —  tf>cos  ( r  ,/2)  (B.)l 

and  the  unconditional  density  can  be  found  as: 

i 

k (ft)  =  J*  I  ri)  2r\dr\ 

max(2cos<«>.  0) 

While  we  are  at  this  point,  we  can  easily  derive  /'( r  !  ft.  rt ) .  the  conditional  density  of 
the  euclidean  distance  from  point  p\  to  p2  given  r x  and  ft.  In  figure  B.5  we  see  that  under  these 
conditions,  Z  must  lie  on  the  line  from  P  to  Q.  Because  points  are  assumed  to  be  distributed 
according  to  the  Poisson  distribution,  the  conditional  density  will  be  uniform  over  the  length  of 
the  line  segment  from  P  to  Q.  and  thus  using  the  derived  expression  for  y  above,  we  can 
write: 

j(:\a.rx)  -\/(\-y)  (B.6) 


B.4  The  Density  of  Ri 


The  calculation  becomes  considerably  more  complex  and  tedious  as  we  progress  to  the 
third  hop.  As  in  the  derivation  for  R2.  we  will  first  calculate  the  distribution,  conditioned 
appropriately,  by  forming  a  ratio  of  two  areas.  Figure  B.6  shows  a  3  hop  configuration  for  a 
given  r\,  ft.  and  z.  Any  node  in  the  third  hop  must  lie  within  that  portion  of  C2  that  lies  out¬ 
side  C|  and  C2.  In  the  figure  we  have  defined  the  angles  and  lengths  that  are  needed  to  calcu¬ 
late  the  area  contained  in  the  two  wedges  A„(y  I  ru  ft.  r)  and  A,(y  |  r  |,  ft.:)  that  are  highlighted 
in  the  figure.  If  the  total  area  lying  only  in  C3  is  T(r\,ft,z).  and  if  ^l+r|  (bounds  to  be 
derived  later)  then: 


/*{/?  3^.v  |  R  i“r |,  0“O,  Z—r] 


A, Ay  I  f|,  ft .  -) +  AAy  \  r<.  ft. :) 

T(r\.  ft.  :) 


The  upper  bound.  y<H-ri.  arises  because  at  v-l+ri  the  circle  of  radius  y  from  pn  becomes 
tangent  to  C\.  For  y’>l-*-r|,  as  shown  in  the  figure,  the  area  is  increased  by  the  circular 
annulus  labeled  fl(y'|/-|,0.r).  The  area  for  this  band  was  derived  in  equation  (B.2),  and  thus 
we  have: 


5(.v’|rt,0,r)  **  A  (l+rt,  l+r()  -  A  (y’.l+rj) 


I 

♦ 
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Au(y  r,.H. Z) 


T(r  ,,w.z) 


Area  only  in  C2 


A, (y  |  r,.«,z) 


Radius  Y'  from  p0 


FIGURE  B.6 

Deriving  the  third  hop  density  function. 


We  can  thus  write: 

I  - — — -  - - -  (B.8) 

/(rt, «, ;) 

where  l+ri^.v<l+/-2<3,  /!/+„(•) +/!„(•),  and  the  upper  bound  is  at  the  point  where 
the  arc  of  radius  l+r2  from  pq  becomes  tangent  to  C2.  and  the  expression  for  r2  can  be  found 
in  the  previous  section.  It  is  clear  that  T(ru  #.  :)  can  then  be  expressed  as: 

T(r |,  W,  r)  »  /</+l((l  +rl  \r  1  +  .4  (l+r|.  l+y|) 

The  only  unknown  quantities  needed  to  calculate  the  distribution  then  are  equations  for 
Au(y  \  rutt.;)  and  A,(  y  I  r]y  9,  z).  This  can  be  done  using  the  defined  angles  and  lengths 


shown  in  figure  B.6.  First  we  will  deri\e  y  «.  r). 


The  objective  of  this  derivation  will  be  to  find  the  length  U(y),  of  the  arc  subtended  by 
angle  6.  Once  this  is  found,  the  integration  of  this  from  yli,.  the  minimum  allowed  value  for 
r j.  to  v.  will  be  equal  to  the  area  of  the  shaded  wedge.  The  derivation  proceeds  by  a  rather 
long  deduction  of  angles  and  lengths: 

r-r “ r\~ +  2r(rcos(«)  — *  (/-f-er'  — 2/-iCCOs(«))_l  2 


r,-=  y2  +  l  —  2vcos(«) 


2  i  i  2 

v  +  1  -  ri 


—  n  —  COS 


y2  +  l  —  2vcos(j3) 


—  (i  =  cos  1 


v-  +  1  —  r 


I  -  r j- 


sin(y+5)  sin  (0) 


sin(<£  +  y)  sin(a) 


,  .  _i  sin(a) 

<A-sm - y 


Thus  we  have  U(y)  -yd>  expressed  in  terms  of  0,  and  r,  and  we  can  write: 


v 

M ( v  I  :)  «■  f  U(r)dr  y%,^y^l+r\ 


We  can  determine  y„  from  figure  B.7,  where  we  see  that: 


“  |ri2  + 1  —  2  T]  cos(0  -  cos  1  (r/2) )  12 


and  thus  yli,  —  max(r2,>„). 


FIGURE  B.7 

Deriving  the  upper  and  lower  bounds. 

The  derivation  of  A,(y  |  r(,  9 ,  r)  proceeds  along  similar  lines,  and 


1  -r)2+>'2-2ri>'cos(c) 

—  c  -  cos  1 

|  r,2+.v2-l ' 

l  2r\v 

r22  =  1  +y2  —  2ycos(d) 

—  cos-1 

•>  ,  ,  ■> 
(,v*  +  l  -r2- 

1  2  y 

n  1 

— *  b“  sin-i 

|  sin(<f)  1 

sin(rf)  sin(S+c+6) 

1  r,  1  C 

—  c 


and  thus  Hy)  —  by  is  also  expressed  in  terms  ot  «,  and  We  then  have: 


<>  f|.  «.  —  J*  Hr)  dr  r,„^  r  ^  1 


where  again  from  figure  B  7  we  have: 


r,  —  /•  | 3  1  —  2  r |  cos(  2tt  —  «  —  cos- 1  <  r/2 )  >  '' 


and  thus  >’4  -  max(/-2<>-,). 


We  can  now  express  the  conditional  distribution  //(/•*  I  rt.  «.  r>  in  terms  of  r n.  and  r. 
The  conditional  density  for  Rj  can  then  be  found  by  differentiating: 

.  dH(r}  I  r]%  «.  :) 

hjir)  |  /•(. «,  r)  » - - -  f B.9) 

dr  3 


Using  equations  (B.l),  (B.5).  (B.6).  and  <B.9)  we  can  find  the  unconditional  density  for  R 3  as: 

/f3(r3)  -  ns  h}(r}  \  :)  ji:  \  H,  r\) k(d\  r\) h\(r\)  ci:  ilH dr\  (B.  10) 


cos  Hr i/2)  <  H  <  77 

0  <  r,  <  1 
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APPENDIX  C 
Derivation  of  IKS.v.y) 


In  this  appendix  we  will  determine  a  formula  for  n(5..v,>’)  and  this  will  permit  us  to 
explain  the  particular  definition  of  the  rates  that  were  defined  in  Chapter  5  (Equations  (5.1)  and 
(5.2)).  It  turns  out  that  with  these  rate  definitions,  determining  the  steady  state  probability  dis¬ 
tribution  for  any  given  topology  is  mathematically  tractable  because  the  Markovian  process 
defined  by  these  rates  define  a  reversible  process  [King69a],  Intuitively  a  reversible  process  X(t) 
is  one  in  which  the  direction  of  time  has  no  effect  on  the  statistics  of  the  process  and  thus  X(r) 
and  X(— t)  have  the  same  probability  distribution.  Finding  closed  form  solutions  for  reversible 
processes  is  simplified  because  that  reversible  processes  satisfy  detailed  balanced  equations.  Let 
IKS)  be  the  steady  state  probability  of  state  S  and  let  q(S.S')  be  the  rate  at  which  transitions 
between  states  5  and  S'  occur.  The  detailed  balanced  equations  state: 

n(S)  q(S,S')  =  IKS')  q(S\S)  (C.l) 

In  general  Markov  processes  satisfy  global  balance  equations  that  equate  the  total  probability 
flux  entering  a  state  to  that  leaving  the  state.  These  equations  state  that: 

£  ri(S')  q(S\S)  -  £  H(S)  «(S.S')  (C.l) 

S'  S' 

We  can  depict  equation  (C.l)  in  Figure  C-l  where  the  equation  states  that  the  net  probability 
flux  across  the  cut  is  zero.  For  reversible  systems,  equation  (C.2)  shows  that  the  net  flow 
across  arcs  connecting  any  two  states  (the  cut  in  Figure  C.2)  is  equal  to  zero.  Equation  (C.2) 
can  be  used  to  simplify  finding  closed  form  solutions  if  reversibility  can  be  proven,  and  a  useful 
tool  to  do  this  is  Kolmogorov's  criteria. 


Theorem:  (Kolmogorov’s  criteria) 

Let  S',  S’,  •  •  ■  ,  S*‘,  S**1 —  S'  be  any  sequence  of  states.  Then  a  Markov  process  defined  over 
these  states  is  reversible  if  and  only  if  the  rate  transitions  satisfy: 

q(S\  S2)  q(S2,  S3)  -  q(Sk~KSk)  q(Sk .  Bk*')  - 

q(Sk*1,  Sk)q(Sk .  Sk~])  -  q(S\  S 2)  q(S\  S]) 


To  determine  n(S,.x,.y)  for  our  system,  we  will  first  need  to  make  some  preliminary 
definitions.  We  can  imagine  that  the  network  for  any  particular  state  S,  is  a  labeled  graph 
G(S)  in  which  the  edges  are  determined  by  the  adjacency  sets  A,  and  the  state  vector  S  deter¬ 
mines  the  labeling.  We  will  define  an  i-arc  (/— 0  or  /—l)  to  be  an  arc  in  G(S)  which  has  both 
nodes  labeled  /.  The  number  of  i-arcs  in  G(S)  will  be  denoted  as  B,(S)  and  will  be  the 
number  of  neighbors  of  node  j  that  are  labeled  i.  Let  T,(S)  be  an  operator  acting  on  state 
5  -  (s,(S),  si(S),  •  •  •  ,5,(5) )  which  complements  the  i'h  component,  thus: 

T,(5)  -  (si(5),s2(S),  •  •  •  ,s,_|(S).  s,(5),  s,+i(S),  •  •  •  ,  s„(S)) 

and  we  will  say  a  transition  is  a  0—1  transition  if  s,(S)  —0  and  for  some  /',  we  have 
Ti (5)  =  S'  (1—0  transitions  are  define  in  the  analogous  way).  We  can  now  prove  a  lemma 
that  will  be  used  to  show  that  the  process  defined  by  rate  equations  (5.1)  and  (5.2)  is  reversi¬ 
ble. 


Lemma:  For  any  state  S,  N,'(S)  -  I  B  (S)  ~  BA  T,(S ))  I 

Proof:  Let  B',(S)  equal  the  number  of  i-arcs  in  the  sub-graph  consisting  of  all  nodes  not  con¬ 
tained  in  A,  |J  I./}.  It  is  clear  that  B',(S)  -  fl',(r,(S)).  Let  k  be  the  number  of  i-arcs  (and 
thus  the  number  of  neighbors  labeled  i)  of  which  node  j  is  a  member  in  state  S  and  k'  be  the 
number  in  state  T,AS).  Since  T,(S)  and  5  differ  in  the  /"'  component  exactly  one  of  k  and  k' 


will  be  non-zero.  We  can  write: 


BAS)  -  k  +  B'  AS) 
BATAS))  «  k'  +  B'  AS) 

and  thus 

I  BAS)  -  BATAS))  I  =  !  k  -  k'\ 


Observe  that  for  a  known  0—1  transition  involving  node  j  we  can  eliminate  the  absolute  value 
sign  in  the  above  to  get  Ni)(S)  —  50(S)  -  Bn(TAS)).  In  a  like  manner,  if  node  /  is  involved 
in  a  l— 0  transition  we  have  N{  (S)  —  BATAS))  -  BAS).  We  are  now  in  a  position  to  prove 
that  the  Markov  process  described  before  is  reversible. 


Theorem:  The  Markov  process  described  by  rate  definitions  (5.1)  and  (5.2)  defines  a  reversible 
Markov  process. 

Proof:  We  will  show  that  Kolmogorov's  criteria  is  satisfied.  For  any  closed  sequence  of  states 

S].  S2,  ■  ■  ■  .  S‘ .  Kolmogorov’s  criteria  is  trivially  satisfied  if  S'*1  **  7", (S')  for  some  /  since 

the  probability  of  two  or  more  components  of  S'  and  S'+l  differing  is  zero.  Thus  assume  that 
the  sequence  of  states  consists  of  single  component  transitions.  Since  there  are  only  two  types 
of  transitions,  all  closed  cycles  must  contain  an  even  number  of  states.  Let 
S’,  S2,  •  •  •  ,  S2m,  S2m+1  -  S’  be  such  a  sequence.  Let 

-  <7 ( S’.  S2)  q(S2,  S3)  •••  q(S2'n,  S2"1^1) 

-  q(S2m*\  S2m)  ••• q(S\  S2)  q(S2,  S') 

We  must  show  that  -  ^2  and  will  refer  to  their  corresponding  sequences  as  the  forward  and 
backward  sequence  respectively.  We  first  note  that  since  there  are  an  equal  number  of  0—1 
and  1—0  transitions  in  'I't  and  'Ri  we  can  write  —  m’"  yo'” where  <t>,  contains  all  the  fac- 


* 

« 

t 


tors  of  .v  and  y.  We  thus  must  show  that  <J>|  =  4>2.  Concentrating  first  on  the  exponent  of  x  in 
these  expressions,  define  the  two  sets: 

C,  =  { <S  .S'*1)  I  3 /  T,(S‘)  =  S'*1.  s.(S  )  =01 

C2  =  !(S'+I,S') !  3/  TAS'*')  =  S'.  s,(5'*')  =  0) ! 

In  words,  C\  contains  the  0—1  transitions  for  the  forward  sequence  and  C2  for  the  backward 
sequence.  Since  a  0—1  transition  in  the  forward  sequence  is  a  1—0  transition  in  the  backward 
sequence,  we  know  that  C\  P|  C2  =  0.  Since  C\  only  contains  0—1  transitions  we  can  write, 
using  the  lemma,  the  exponent  of  x  in  the  forward  E\,  and  backward  £2.  sequences  as: 

E,  =  X  Bo (S)  -  B»(S') 

( s .  S')  =  c, 

Suppose  now  that  5  ,  S'*1,  •  •  •  S'*k  is  a  sub-sequence  of  states  from  the  forward 
sequence  such  that  (S',  S'4'1)  €  C|.  Since  the  portion  of  E\  for  this  sequence  alternates  sign, 
the  sum  telescopes  and  we  can  write  for  this  sub-section: 

1  flo(S')  -  fl0<S'+l)  -  B0(SA  -  BAS'^) 

r—i 

and  since  for  this  sequence  of  states,  (S'4-1,  S')  £  C2  j  =  /,  /+1,  •  •  •  .  i+k—l.  we  can  conclude 
that  £2  is  not  affected  by  this  telescoping.  It  is  thus  clear  that  in  calculating  £|  and  £2  we  only 
have  to  look  at  sections  of  either  sequence  where  a  change  in  the  type  of  state  transition  occurs. 
Denoting  these  places  of  state  transition  as  S"  we  have: 

£1  =  [£o(5',)-fi0(5’2)l  +  [£o(S':,)-5o(5'4)]+  •  •  ■  +[S0(S'2*-|)-S0(S'2A)] 

E1-lB9(S,x)-BdS,2k)]  +  [B0(S'lk-')-B0(S'ik-l)]+  ■  ■  ■  +  [fl„(S'3)-flo<S’:)] 

and  thus  see  that  E\  —  £2.  An  analogous  argument  can  be  made  for  the  exponent  of  v  con¬ 
cluding  the  proof. 


Knowing  that  the  Markov  process  is  reversible  allows  us  to  use  the  detailed  balance 
equations  (5.1)  to  prove  the  following  theorem. 


Theorem  C. /.The  equilibrium  probability  for  state  5  with  transitions  rates  defined  by  equations 
(5.1)  and  (5.2)  is  given  by: 

n(5..v,y)  -  C  pu'S)  x~Bl)iS)  yB',S) 

where: 


A  /(5>  =  £  s,(5> 

,-i 

_  ,  -8nlOI 

C  =  n(0,.v,y)  x  0 
p  =  Wm 

t-w  / r\  \  f  »/#c\  fli  <  S)  1  ”1 

n(o,.v,y)  -[x  Xp  x  y  j 


Proof:  To  avoid  cumbersome  notation  let  n(5,.v,y)  =  fl(5).  Since  the  process  is  reversible  we 
can  use  the  detailed  balance  equations  to  state: 

n  (5)  q(US),S) 
n (T,(S))  "  q(S,  T,(S )) 

If  5,(5)  -  1  we  can  write  equation  (C.3)  as: 

n(S)  vm>  v((S) 

YUTjS))  ~PX  y 


which  by  using  the  lemma  can  be  re-written  (for  the  case  of  a  0— I  transition;  as: 


n(S)  80<  r  <s>>  -  Bn(S)  8,  (S>  -  8,1  r, (51) 

m-isiy-v* 


5,(5)  =  1 


(C.4) 


We  can  use  this  relationship  to  write  fl(5)  in  terms  of  fl(0)  by  telescoping  a  product  of  rate 


ratios.  Suppose  m.  /3,  •  ■  ■  ,  /w,  s,  are  the  indices  of  S  which  are  equal  to  1.  Define  the  fol¬ 
lowing  operator: 

F0-\ 

F,  =  T.'Fi- 1  y  -  1.2.  •  ■  ■  ,  .V/(S) 

Observe  that  F0(S)  =  S  and  FU(S)(S)  =  0  (a  vector  of  all  zeros).  We  can  then  write  equation 
(C.3)  as: 

n(S)  MF0(S))  n (F,(S))  n(Fv/(51_i(S)) 

n(0)  “  n(F,(5))  n(F,(s)i  n (f„,S)(s)) 

q(Fl(S),F0(S))  q(F2(S),Ft(S))  q(FutSi(S).F„,s>-i(S)) 

"  ?(F0(S).F|(S))  ?(F,(5).F2(S))  ”'?(f|/(Si-|(S).Flf,.s,(5)) 

Which  can  be  simplified  by  using  equation  (C.4)  to: 

n,,,  V  80<f,_,ts>>- Z 

lH-S)  U(S)  .-0  -0 

INO)  p 

Once  again  these  sums  telescope  and  we  are  left  with: 

n(S)  1/(5,  B0lFit(s)lS))-B0iFf)iS))  81(F0(Sn-a,(fw(s)(S)) 

n(0)  p  y 


BrtIO)  -  80fS>  8,15)-  8,10) 


We  can  thus  write: 

n(S)  -  n(0)  p 1/(51  xtf°"  -  y 
and  using  the  fact  that  0)  —  0  and  defining  C  —  11(0)  x  Bf)i0>  we  finally  obtain: 

n (S) 


(C.5) 


We  can  determine  FI (0)  by  using  the  normalization  constraint  £  IKS)  =  1  to  determine  the 

s 

equation  for  n(0)  given  in  the  statement  of  Theorem  C.l. 


131 


References 


(Abra70a)  N.  Abramson,  “The  ALOHA  System  —  Another  Alternative  for  Computer 
Communications,”  AFIPS  Conference  Proceedings,  FJCC  37,  pp.281-285 
(1970). 

[Abra73a]  N.  Abramson,  “Packet  Switching  with  Satellites,”  AFIPS  Conference  Proceed¬ 

ings,  NCC 42,  pp.695-703  (June  1973). 

[Abra77a]  N.  Abramson,  “The  Throughput  of  Packet  Broadcasting  Channels,”  IEEE 
Transactions  on  Communications  COM-250),  pp.  117-128  (January  1977). 

[Akav78a]  G.  Y.  Akavia,  Hierarchical  Organizations  of  Distributed  Packet-Switching  Com¬ 
munication  Systems,  Computer  Science  Department.  University  of  California, 
Los  Angeles  (March  1978).  Ph.D.  Dissertation. 

[Avri76a]  M.  Avriel,  Nonlinear  Programming:  Analysis  and  Methods,  Prentice-Hall,  Inc., 

Englewood  Cliffs,  New  Jersey  (1976). 

[Cant74a]  D.  G.  Cantor  and  M.  Gerla,  “Capacity  Allocation  in  Distributed  Computer 

Networks,”  Proceedings  of  the  Seventh  Hawaii  International  Conference  on  Sys¬ 
tems  Sciences,  pp.l  15-1 1 7,  University  of  Hawaii  (January  8-10  1974). 

[Davi80a]  D.  H.  Davis  and  S.  A.  Gronemeyer,  “Performance  of  Slotted  Aloha  Random 

Access  with  Delay  Capture  and  Randomized  Time  of  Arrival,”  IEEE  Transac¬ 
tions  on  Communications  COM-28 (5),  pp.703-710  (May  1980). 

[Fran71a]  H.  Frank,  I.  T.  Frisch,  W.  Chou,  and  R.  Van  Slyke,  “Optimal  Design  of  Cen¬ 

tralized  Computer  Networks,”  Networks  1(1),  pp.43-57,  Wiley  (1971). 

[Frat73a]  L.  Fratta,  M.  Gerla,  and  L.  Kleinrock,  “The  Flow  Deviation  Method:  An 

Approach  to  Store-and-Forward  Communication  Network  Design,”  Networks 
3,  pp.97-133  (1973). 

[Frat80a]  L.  Fratta  and  D.  Sant,  “Some  Models  of  Packet  Radio  Networks  with  Cap¬ 

ture,”  Proceedings  of  the  Fifth  International  Conference  on  Computer  Communi¬ 
cation,  pp. 155-161  (October  27-30,  1980). 


132 


[Howa60a]  R.  A.  Howard,  Dynamic  Programming  and  Markov  Processes,  M.I.T.  Press. 
Cambridge.  Mass.  (1960). 

lKahn78a]  R.  E.  Kahn,  S.  A.  Gronemeyer,  J.  Burchfiel,  and  R.  C.  Kunzeiman, 
“Advances  in  Packet  Radio  Technology,”  Proceedings  of  the  IEEE  66. 
pp.  1468-1496  (November  1978). 

[Kell79a]  F.  P.  Kelly,  Reversibility  and  Stochastic  Networks,  John  Wiley  and  Sons  (1979). 

[Kend63a]  M.  Kendall  and  P.  Moran,  Geometrical  Probability,  Griffin,  London  (1963). 

[King69a]  J.  F.  C.  Kingman,  “Markov  Population  Processes,”  J.  Appl.  Prob.  6  (1969). 

[Klei64a]  L.  Kleinrock,  Communication  Nets:  Stochastic  Message  Flow  and  Delay, 

McGraw-Hill,  New  York  (1964).  Out  of  print.  Reprinted  by  Dover  Publica¬ 
tions,  1972. 

[Klei75a]  L.  Kleinrock  and  F.  A.  Tobagi,  “Packet  Switching  in  Radio  Channels:  Part  I 

—  Carrier  Sense  Multiple-Access  Modes  and  Their  Throughput-Delay  Charac¬ 
teristics,”  IEEE  Transactions  on  Communications  COM-23 (12).  pp.1400-1416 
(December  1975). 

[Klei76a]  L.  Kleinrock,  Queueing  Systems.  Vol  II.,  Computer  Applications,  Wiley- 

Interscience,  New  York  (1976). 

[Klei78a]  L.  Kleinrock  and  J.  Silvester,  “Optimum  Transmission  Radii  for  Packet  Radio 

Networks  or  Why  Six  is  a  Magic  Number,”  Conference  Record.  National 
Telecommunications  Conference,  pp.4.3.1-4.3.5  (December  1978). 

[Lam75a]  S.  S.  Lam  and  L.  Kleinrock,  “Dynamic  Control  Schemes  for  a  Packet 

Switched  Multi-Access  Broadcast  Channel,”  A  FI  PS  Conference  Proceedings, 
NCC 44,  pp.143-154  (May  1975). 

[Lick78a]  J.  C.  R.  Licklider  and  A.  Vezza,  "Applications  of  Information  Networks.” 

Proceedings  of  the  IEEE  66(11)  (November  1978). 

[Litt61a]  J.  Little,  “A  Proof  of  the  Queueing  Formula  Z.-AWV  Operations  Research 

9(2),  pp.383-387  (March  1961). 

[Liu80a]  J.  Liu,  “A  Distributed  Routing  Algorithm  in  Mobile  Packet  Radio  Net¬ 

works,”  Technical  Report,  University  of  Illinois.  Urbana,  Ill.  (1980). 


[Meis7 la]  B.  Meister,  H.  Muller,  and  H.  Rudin,  “New  Optimization  Criteria  for 

Message-Switching  Networks,”  IEEE  Transactions  on  Communication  Technol¬ 
ogy  COM-19,  pp.256-260  (June  1971). 

[Mine70a]  H.  Mine  and  S.  Osaki.  Markovian  Decision  Processes.  American  Elsevier 

(1970). 

[Mo!181a]  M.  L.  Molle,  “Unifications  and  Extensions  of  the  Multiple  Access  Communi¬ 

cations  Problem.”  CSD  Report  No.  810730  (UCLA-ENG-81 18).  Computer 
Science  Department.  University  of  California,  Los  Angeles  (July  1981) 
Ph.D.  Dissertation. 

[Robe72b]  L.  G.  Roberts.  “ALOHA  Packet  System  With  and  Without  Slots  and  Cap¬ 
ture.”  ASS  Note  8  (NIC  11290).  ARP  A  Network  Information  Center.  Stan¬ 
ford  Res.  Inst..  Menlo  Park.  Ca.  (June  1972)  reprinted  in  Computer  Com¬ 
munication  Review,  Vol.  5.  pp.  28-42  (April  1975) 

[Robe72a]  L.  G.  Roberts,  “Extensions  of  Packet  Communication  ~  chnology  to  a  Hand 

Held  Personal  Terminal,”  AFIPS  Conference  Proce  ogs.  SJCC  49.  pp  295-298 
(1972). 

[Silv80a]  J.  Silvester,  “On  Spatial  Capacity  of  Packet  Radio  Networks.  '  UCLA-ENG- 

8021.  Computer  Science  Department.  University  of  California.  Los  Angeles 
(May  1980).  Ph.D.  Dissertation. 

[Toba74a]  F.  A.  Tobagi,  “Random  Access  Techniques  for  Data  Transmission  over 
Packet  Switched  Radio  Networks,”  UCLA-ENG-7499.  Computer  Science 
Department,  University  of  California,  Los  Angeles  (December  1974).  Ph.D. 
Dissertation. 

(Toba75a]  F.  A.  Tobagi  and  L.  Kleinrock,  "Packet  Switching  in  Radio  Channels:  Part  II. 

—  The  Hidden  Terminal  Problem  in  Carrier  Sense  Multiple-Access  and  the 
Busy  Tone  Solution,”  IEEE  Transactions  on  Communications  COM-23U2). 
pp.1417-1433  (December  1975). 

[Toba80a]  F.  A.  Tobagi,  “Analysis  of  a  Two-Hop  Centralized  Packet  Radio  Network- 
Part  I:  Slotted  Aloha,”  IEEE  Transactions  on  communications  Com-28(2). 
pp.196-207  (February  1980). 

[Toba80b]  F.  A.  Tobagi,  "Analysis  of  a  Two-Hop  Centralized  Packet  Radio  Network- 
Part  II:  Carrier  Sense  Multiple  Access,”  IEEE  Transactions  on  Communications 
Com*28(2),  pp.208-216  (February  1980). 


